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Das Problem der Ordnungs—Unordnungs- 


übergänge in Kristallen 
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H. JAGODZINSKI 
Marburg 
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I. Einleitung 


Das Problem der Ordnungs—Unordnungsübergänge gehört zu den 
interessantesten der Kristallstrukturforschung und ist im Verlauf der Ent- 
wicklung besonders der vergangenen zwanzig Jahre von vielen Autoren 
sowohl experimentell als auch theoretisch untersucht worden. Viele physi- 
kalische und physikalisch-chemische Eigenschaften der Kristalle sind weit- 
gehend vom Ordnungszustand abhängig, und da die damit verbundenen 
Fragen nicht nur vom wissenschaftlichen, sondern auch vom technischen 
Standpunkt sehr wichtig sind, ist es kein Wunder, daß auch in den Kriegs- 
jahren eingehende Untersuchungen ausgeführt wurden. Insgesamt gesehen 
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bilden aber alle theoretischen und experimentellen Ergebnisse ein sehr 
inhomogenes Gebiet, und es ist deshalb nicht ganz leicht, das ganze vor- 
liegende Material einheitlich zusammengefaßt darzustellen. Hinzu kommt 
noch, daß ein großer Teil der vorliegenden Arbeiten in oder unmittelbar 
nach dem Kriege im Ausland erschienen und somit dem deutschen Leser 
nicht leicht zugänglich sind. Es scheint deshalb notwendig, den ganzen 
Fragenkomplex einmal möglichst unter einheitlichen Gesichtspunkten dar- 
zustellen. Ich will aber dabei nicht das Ziel verfolgen, einen vollständigen 
Bericht über alle in dieses Gebiet fallenden Arbeiten zu geben. Das ist allein 
aus Raumgründen nicht möglich. Es kommt mir vielmehr darauf an, die 
bedeutendsten Fortschritte der theoretischen Untersuchungen kritisch zu 
beleuchten und vor allen Dingen die Grenzen ihrer Anwendbarkeit klar 
hervorzuheben. Die zahllosen experimentellen Ergebnisse werden hier nur 
gelegentlich herangezogen, um sie mit dem zu vergleichen, was die Theorie 
über den entsprechenden Fall aussagt. 

Leider stößt die theoretische Bearbeitung des Gebietes auf große 
Schwierigkeiten nicht nur rein formal mathematischer, sondern auch begriff- 
licher Art. Ich habe deshalb angestrebt, überall den physikalischen Gedan- 
kengang so klar wie möglich herauszuheben und die Berechnungen nur 
so weit in den Text einzuarbeiten, wie sie mir zum Verständnis der großen 
Linie der vorliegenden Theorien notwendig erscheinen. Daß dies an einigen 
Stellen auf Kosten der mathematischen Exaktheit der Darstellung ging, 
bitte ich den mathematisch interessierten Leser mir zu verzeihen. Die 
Arbeit wird sehr durch einen von NIx und SHOCKLEY (1) veröffentlichten 
Bericht aller bis zum Jahre 1938 erschienenen Arbeiten über Ordnungs— 
Unordnungsübergänge in Legierungen erleichtert. Die Spezialisierung auf 
Legierungen ist kaum ein Nachteil, da sich fast alle bisherigen Unter- 
suchungen auf dieses Gebiet beschränkten. Ich werde jedoch für den fol- 
genden Bericht nicht die Kenntnis der Arbeit von Nıx und SHOCKLEY 
voraussetzen, sondern die wichtigsten theoretischen Teile kurz zur Ein- 
führung darstellen. Das wird im Kapitel II, $ 1, 3 und 4 getan werden; 
der Leser, dem diese Arbeiten bereits bekannt sind, kann diesen Teil der 
Arbeit ruhig überspringen. 

Während sich Kapitel II im wesentlichen mit Berichten über bereits 
erschienene Arbeiten befassen wird, soll im Kapitel III bei der Behandlung 
des röntgenographischen Beugungsproblems versucht werden, eine Brücke 
zu den Ergebnissen der statistischen Thermodynamik zu schlagen. In diesem 
Teil beschränkt sich diese Arbeit nicht auf reine Berichte, sondern es werden 
auch neue Gedankengänge und Berechnungen diskutiert, um abschätzen zu 
können, wie weit die Theorie mit dem Experiment im Einklang steht. Das 
Gebiet der sogenannten Rotationsumwandlungen wird hier nicht näher 
besprochen werden, weil die wichtigsten Fortschritte bereits vor dem Kriege 
erzielt wurden; vom theoretischen Standpunkt stehen hier die Verhältnisse 
in voller Analogie zu den in dieser Arbeit behandelten Problemen. 

Schon bald nach der Eröffnung der fruchtbaren Periode der Erforschung 
der Struktur der Materie mittels Röntgenstrahlenbeugung ergab sich, daß 
der Bau der Kristalle bei weitem nicht so ideal sein konnte, wie man es 
sich ursprünglich vorgestellt hatte. Man führte den Begriff des Mosaik- 
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kristalls ein, entwarf Vorstellungen über randliche Störungen und mußte 
zur Klärung vieler physikalischer Eigenschaften auch Störstellen im Gitter 
annehmen. Erinnert sei hier nur an die schönen Untersuchungen der licht- 
elektrischen Leitfähigkeit, der Ionenleitfähigkeit, der Mischkristalle und 
ihrer Überstrukturen, des spezifischen Widerstandes, der Lichtabsorption 
usw. Daß solche Fehlstellen vorhanden sind, bedarf keiner Frage, lediglich 
die Art der in einem Problem vorliegenden Fehlstellen ist unter Umständen 
eine strittige Frage. Zur Klärung der experimentellen Ergebnisse haben sich 
bestimmte Vorstellungen entwickelt, die schließlich von WAGNER und 
SCHOTTKY (2) (3) (4) in ihren grundlegenden Arbeiten über die Theorie der 
geordneten Mischphasen ihren Niederschlag fanden. Sie gliederten die 
Fehlstellen in 


a) Fehlbesetzungen (Substitutionstypus) 
b) Leerstellen (Subtraktionstypus) 
c) Zwischengitterstellen (Additionstypus) 


ein, wobei wohl die dazu entworfenen geometrischen Vorstellungen natürlich 
nur Modelle sind und in ihrer strengen Geometrie nicht genau genommen 
werden dürfen. Es soll hier nicht im einzelnen diskutiert werden, wie weit 
die einzelnen Vorstellungen im Kristall realisiert sind oder nicht. Sicher ist 
nur, daß es verschiedene Arten von Fehlstellen gibt, die sich in wesentlichen 
Punkten in ihrer Geometrie unterscheiden. — Erwähnt seien in diesem 
Zusammenhang noch die für die plastische Verformung so wichtigen Ver- 
setzungen der Atome in Metallen, die von BraGG und Nye (5) mittels des 
bekannten Seifenblasenmodells auf elegante Weise evident gemacht wurden. 


Es ist einleuchtend, daß eine Fehlstelle im Gitter ein energetisch un- 
günstiger Zustand ist, dessen Existenz im allgemeinen zwei Gründe haben 
kann: 

1. Die Fehlstelle ist eine Folge der Temperaturbewegung der Atome 
oder Moleküle. (Die Wärmeschwingungen bringen in Abhängigkeit 
von ihrer Amplitude mehr oder weniger häufig die Energiedifferenz 
zur Bildung der Fehlstelle auf.) 

2. Die Fehlstelle ist aus Gründen der geometrischen Schwierigkeiten 
nicht ohne weiteres aus dem Gitter zu entfernen, sei es, daß die 
Schwingungsenergie zu gering ist, um durch Platzwechsel der Atome 
einen Abbau der Fehlstelle zu bewirken, oder daß es Atome oder 
Moleküle gibt, die sich selbst oder in bezug auf ihre Umgebung mit 
der Geometrie des Gitters schlecht oder gar nicht vereinbaren lassen. 
Hierunter fallen also insbesondere auch die 2 an Korn- 
grenzen. 

Ist die Anzahl der Fehlstellen gering zur Anzahl der möglichen Fehl- 
besetzungen, so kann man sie in erster Näherung als unabhängig voneinander 
ansetzen, und man gelangt mittels des BoLtzmann’schen Gesetzes der Ver- 
knüpfung von Entropie und Wahrscheinlichkeit leicht zu einer einfachen 
Beziehung für den wahrscheinlichsten Zustand des Kristalls als Funktion 
der Temperatur. Ist das aber nicht der Fall, so müssen die Wechselwirkungen 
der Fehlbesetzungen untereinander zur Berechnung herangezogen werden. 
Damit führt die Fragestellung unmittelbar auf das Gebiet der sogenannten 
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kooperativen Erscheinungen, deren mathematische Behandlung nicht sehr 
einfach durchzuführen ist. 

Die erste Entwicklung der Theorie der Fehlordnung, die von Autoren 
wie GORSKY (6), DEHLINGER und GRAF (7), WAGNER und SCHOTTKY (2) (3) 
(4), BoRELIUS (8) eingeleitet wurden, ging zunächst von der klassischen 
Thermodynamik aus, in die das statistische Problem mehr oder weniger 
eingearbeitet wurde. Die Entwicklung, die schließlich in den klassischen 
Arbeiten von BRAGG und WirLıams (9) (10) (11) ihren Abschluß fand, kann 
aber heute weitgehend als überholt betrachtet werden; ihre Arbeitsweise, 
die Energie einer Fehlbesetzung unabhängig von den Schwankungen im 
Kristallinneren durch eine gemeinsame Funktion des Ordnungszustandes 
anzusetzen, dürfte dem tatsächlichen Zustand nur in den Grenzfällen der 
völligen Ordnung bzw. Unordnung nahekommen. Eine wirklich befriedi- 
gende Lösung kann man nur von der Statistik her erwarten. Solche Lösungen 
sind heute, wenn auch in etwas vereinfachter Form, herbeigeführt worden. 
Diese Methoden haben sich den klassischen thermodynamischen Verfahren 
als überlegen erwiesen, wenn auch die mathematische Behandlung der frü- 
heren Methoden sowie die Erstreckung der Berechnungen auf Nichtgleich- 
gewichtszustände bedeutend einfacher sind. Da hier die Arbeiten von BRAGG 
und WILLIAMS zusammen mit den modernen Arbeiten diskutiert werden, 
kann an dieser Stelle verzichtet werden, auf die Zusammenhänge und die 
Unterschiede zwischen den beiden Vorstellungen näher einzugehen. 

Ein strenges Einteilungsprinzip der verschiedenen Fehlordnungstypen 
ist eigentlich nicht zu geben; aber immerhin kann man auf Grund der vor- 
liegenden experimentellen Ergebnisse gewisse idealisierte Typen annehmen; 
man muß sich jedoch immer klar darüber sein, daß sie in der Idealform 
kaum oder nur in seltenen Fällen anzutreffen sind. Tabelle 1 soll einen 
Überblick über die in Frage kommenden Arten geben: 


Tab. 1: Einteilungsschema für Fehlordnungstypen 


(eindimensionale) 
Fehlordnung (zweidimensionale) 
es a 
I 
Erhaltung der Gittergeometrie Veränderung der Gittergeometrie 
Substitutions- Leer- Überstruktur- Entmischung Zwischen- 1. Mehrere Be- 


mischkristall stelle bildung und und Übersruk- gitterplätze setzungsmög- 
Entmischung turbildung mit Versetzun- lichkeiten un- 


ohne Änderung Endgliedern gen gleichwertiger 
der geometri- stark abwei- Gitterplätze. 
schen Atom-  chender Gitter- 2. Völlige Ände- 
lagen geometrie rung der Git- 
tergeometrie 
durch Gitter- 
typen ver- 
schiedener 


Symmetrie. 
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Fehlordnung kann im dreidimensionalen Kristall in einer, zwei oder 
drei Translationsrichtungen auftreten. Entsprechend wollen wir die Kristalle 
als ein-, zwei- oder dreidimensional fehlgeordnet ansehen. Eindimensional 
fehlgeordnet würde dann bedeuten, daß in zwei Translationsrichtungen 
strenge Ordnung (Netzebenen) herrscht, die dritte aber Fehlordnung auf- 
weist. Wir kennen sehr viele praktische Beispiele dieser Art (Zinkblende— 
Wurtzit, SiC, Maucherit, organische Schichtstrukturen, Graphit usw.). Ein 
zweidimensional fehlgeordneter Kristall würde aus streng geordneten Ketten 
bestehen, deren Anordnung zueinander nach mehr oder weniger statistischen 
Prinzipien erfolgt. (Hierzu gehören Faserstrukturen, z. B. Chrysotil.) Wir 
unterscheiden in diesen drei Hauptgruppen nun zwischen 


1. Fehlordnung unter strenger oder annähernder Erhaltung der Gitter- 
geometrie und 
2. Fehlordnung unter Störung der Gittergeometrie. 


Diese Unterscheidung ist insofern besonders wichtig, weil jeder Unordnungs- 
vorgang entscheidend von der Geometrie des Gitters bestimmt wird. Neigt 
z. B. ein bei hohen Temperaturen stabilisierter Mischkristall bei tieferen 
Temperaturen zur Entmischung oder Überstrukturbildung, so ist es von 
sehr einschneidender Bedeutung, ob die Überstruktur oder das Auszu- 
scheidende mit der Geometrie des Gitters einigermaßen vergleichbar ist. 
Ist das nicht der Fall, so interessiert hier besonders die Kinetik des Vor- 
ganges, da die hohen zu überwindenden Potentialschwellen oft die Er- 
reichung des stabilen Endzustandes verhindern. Diese Art der Entmischung 
hat eine große Bedeutung bei der Herstellung aushärtbarer Legierungen 
gewonnen, bei denen immer solche Verhältnisse vorliegen. Die Härte- 
steigerung wird dabei durch die Behinderung der Gleitfähigkeit der Kristalle 
infolge der Störstellen und den damit verknüpften Spannungen im Gitter 
hervorgerufen. . 

Im folgenden soll nun zunächst ein Überblick über die vorhandenen 
theoretischen Ansätze der Fehlordnung gegeben werden. Es kann kein 
Zweifel darüber bestehen, daß alle Ansätze nur dazu beitragen können, das 
charakteristische Verhalten des Übergangs zu erläutern; von einer quanti- 
tativen Anwendbarkeit kann nur in einzelnen Idealfällen die Rede sein. 
Trotzdem kommt diesen Theorien eine sehr hohe Bedeutung zu, weil sie 
das grundsätzliche Verhalten der physikalischen Eigenschaften richtig 
wiedergeben. Ich werde mich in der Darstellung nicht an die historische 
Entwicklung halten, sondern die verschiedenen theoretischen Methoden 
zeitlich unabhängig, dafür aber in sich geschlossen entwerfen. 


II. Die Theorien der Ordnungs—Unordnungsübergänge 


1. Die Fernordnungstheorie von BRAGG und WILLIAMS 
BraGG und WirLıams (9) (10) (11) haben in sehr einfacher Weise eine 
Theorie entwickelt, die bereits weitgehende Aussagen zu machen gestattet. 
Sie beschränken ihre Untersuchungen in ihrer Anwendbarkeit auf Metall- 
legierungen. 
Gegeben sei ein Gitter, das mit zwei Atomarten (A, B) besetzt ist, die 
zueinander in einem bestimmten stöchiometrischen Verhältnis, z. B. AB 
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oder AB,, stehen. Es wird angenommen, daß bei hohen Temperaturen der 
Kristall dem idealen Mischkristall sehr nahe kommt, während bei tiefen 
Temperaturen eine Neigung zur Überstrukturbildung besteht. Beim abso- 
luten Nullpunkt ist der streng geordnete Kristall stabil; die A-Atome 
besetzen bestimmte Gitterplätze (a-Gitterplätze), desgleichen die B-Atome 
(b-Gitterplätze); (z. B. Atombesetzungen im NaCl- oder CsCl-Gittertyp). 
BrAGG und WILLIAMS setzten nun voraus, daß die A- und B-Atome ‚‚wissen‘“, 
welches für sie die ‚richtigen‘ (a-, b-Gitterplätze) oder die „falschen“ 
(b-, a-Gitterplätze) sind. Weiterhin nehmen BrAGG und WILLIAMS an, daß 
die Gittergeometrie für alle möglichen Besetzungen nicht verändert wird. 
Nennen wir nun r, den Bruchteil der A-Atome auf a-Plätzen und ent- 
sprechend r, denjenigen der B-Atome auf b-Plätzen, so ist 1—r, und 
1 — r» der Bruchteil der A- bzw. B-Atome auf falschen (b- bzw. a-) Gitter- 
plätzen. Mit F, und F', bezeichnen wir den Bruchteil der a- bzw. der b-Gitter- 
plätze (F, + F» = 1). Haben wir rein zufällige Verteilung der A- und 
B-Atome, so befinden sich offenbar von den A-Atomen der Bruchteil F, 
auf a- und F, auf b-Gitterplätzen. Das gleiche gilt auch für die B-Atome. 
BRAGG und WiırLıams definieren nun den Fernordnungsparameter S wie 
folgt: 
N — 15 ig F) 
) Se Seen 


Man sieht, daß unter Berücksichtigung der vorhin angestellten Be- 
trachtungen für die zufällige Verteilung S = 0 und für die geordnete Ver- 
teilung S = 1 wird. Der Ordnungsgrad wird hier also durch die Anzahl der 
„falsch‘‘-sitzenden A-Atome definiert, damit ist natürlich auch die Anzahl 
der Fehlbesetzungen für die B-Atome festgelegt. 


Die weitere Aufgabe ist nun, die Kristallenergie E als Funktion des 
Ordnungsstandes S abzuleiten (E=E (S)). Diese Aufgabe wird von BRAGG 
und WILLIAMS nicht exakt gelöst, sondern durch einen Lösungsansatz er- 
setzt, den sie wie folgt begründen: Im völlig ungeordneten Zustand, bei dem 
also die Unterscheidung zwischen richtigen und falschen Gitterplätzen 
keinen Sinn mehr hat, wird durch die Umbesetzung eines A-Atoms von einem 
a- auf einen b-Platz kein Energieaufwand benötigt (abgesehen von der 
Energie zur Überwindung der Potentialschwelle eines Platzwechsels, die ja 
nach Beendigung des Platzwechsels wieder gewonnen wird). Im geordneten 
Zustand dagegen muß die Arbeit V, geleistet werden, um eine Fehlbesetzung 
durchzuführen. Der Zusammenhang der aufzuwendenden Arbeit V für einen 
Zwischenzustand der Ordnung S wird nun von BRAGG und WILLIAMS durch 
den Ansatz 


(IL, 2) Wev.s 


definiert. Es soll hier davon Abstand genommen werden, über die Zulässig- 
keit dieses Verfahrens zu diskutieren. Sicher ist, daß V, ausgehend von O,V, 
auf irgendeinem Wege in Abhängigkeit von S erreicht, der sicher keinen 
sehr komplizierten Funktionszusammenhang zeigen wird. Da in dieser 


Theorie nur die qualitativen Ergebnisse interessieren, ist die genaue Kenntnis 


dieses Zusammenhanges ziemlich belanglos. 
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Die Berechnung der Energie kann nun in einfacher Weise erfolgen, in 
dem man sich überlegt, wie viele neue Fehlbesetzungen vorgenommen werden 
müssen, um den Ordnungszustand S um dS zu verändern. Haben wir ins- 
gesamt N A- und B-Atome (also auch Gitterplätze) im Gitter, so ist, weil 
mit einer Fehlbesetzung immer ein Paar verändert werden muß, die Anzahl 
der Fehlbesetzungen gleich NF, (1 —r,). Ausgedrückt durch S ergibt das 
gemäß (II, 1) 

Ber DeNDgTe m (1 Pn)SIeNMIR,d> 8: 
Nach S differenziert, erhalten wir 

da[NF, (1 —r,)]=—NF,(1—F,)dS=-— NF,FydS. 
Anwendung von Gleichung (II, 2) ergibt für dE 
dE = —NV,F,F,SdS. 


Daraus gewinnt man durch Integration vonS—=1bisS 


15:3) E(S)=4NV,F,F, (1 — 8), 
mit E(O)=E,=4NV,F,F, wird Gleichung (II, 3) 
(II, 4) E(S)=E,(1— 8). 


Damit ist das erste Problem, die Auffindung der Kristallenergie als 
Funktion des Ordnungszustandes, als gelöst zu betrachten. Die zweite 
Aufgabe besteht nun darin, den für eine bestimmte Temperatur T wahr- 
scheinlichsten Ordnungszustand des Kristalls zu ermitteln. Man erreicht 
die hierzu notwendige Bedingungsgleichung am einfachsten über das Mini- 
mum der freien Energie F, das ja die Bedingung für das thermodynamische 
Gleichgewicht darstellt. F ist definiert durch 
(IL, 5) F(S,T)=E($)—To($), 
wobei ® (S) die Entropie des Ordnungszustandes darstellt. © (S) ist nun 
mit der Wahrscheinlichkeit W (S) des Ordnungszustandes durch die BoLTz- 
MANN’sche Beziehung 

8(S)=klnW (S) 
verknüpft. W (S) läßt sich aber leicht mit Hilfe der bekannten kombi- 
natorischen Formen ermitteln, wenn man die A- und B-Atome alle möglichen 
Permutationen auf den richtigen und falschen Plätzen durchlaufen läßt und 
bedenkt, daß die Permutationen der ‚richtigen‘ und ‚falschen‘ A-, B- 
Atome unter sich keine neue Anordnung erzeugen. 


1,6) WS) In} Be) N 
1,9 DER ENT m) EN HANIA—n) Fon] 
Anwendung der STIRLIng’schen Formeln und Ersetzen von r, und ry 
durch S gemäß Gleichung (II, 1) ergibt 
(EL, 7) In W (5) =—N {Fr (F, + F,S) In (F, u F,S) 7 F,F, (1 —S) 
[In F, (1 Er S) + In Fy(l a“ S)] ir Fy (F, im F,S) In (F1 Sr F,S)}. 
Für das Minimum der freien Energie gilt bei einer bestimmten Tem- 
peratur 


(I, 8) 


—=.0, 


oF dE d [InW(S)] 
—— m kT — = 
08 ds ds 
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Wir nehmen also an, daß die Energie E temperaturunabhängig ist. 
Differentiation von Gleichung (IL, 3) und (II, 7) nach S und Einsetzen der 
Differantialquotienten in (II, 8) ergibt 
1—F,(1—S$) 

F,(1—S) 


BR dESpen 
FB, (1 —8) 7 


—NV.SF,F, + KENEAR,| In san 


Daraus erhält man die Beziehung 


(II, 9) I 


F, (1—8$) TuS KESaET 


oder in Worten 


Anzahl der A-Atome auf B-Blätzen x Anzahl der B-Atome auf A-Plätzen E* 
Anzahl der A-Atome auf A-Plätzen x Anzahl der B-Atome auf B-Plätzen 
exp — V/KT. 


Mit Hilfe der Gleichung (II, 9) ist es nun nicht mehr so schwer, den 
Fehlordnungsgrad als Funktion der Temperatur auszureehnen. Immerhin 
ist der funktionelle Zusammenhang nicht in eine übersichtliche Form zu 
bringen. Im allgemeinen löst man deshalb Gleichung (II, 9) graphisch, in 
dem man einerseits 


V. 
= mS 
en 
u oz Fr, (8% 


als Funktionen von S aufträgt. Der Schnittpunkt der beiden Kurven ergibt 
dann den Gleichgewichtszustand für eine bestimmte Temperatur (s. Fig. 1). 


S 


fı (5) 


k tıls 
P 


0 #1, #2 


Fig. 1. Graphische Ermittlung des Ordnungszustands als Funktion der Temperatur. 


Die Frage, ob ein Minimum, Maximum oder Wendepunkt der freien Energie 
vorliegt, kann dann mit der üblichen Methode der Betra chtung der Differen- 
tialquotienten entschieden werden. 
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In Fig. 2 ist das prinzipielle Verhalten des Ordnungszustandes als 
Funktion der Temperatur für den AB-Fall graphisch dargestellt. Man sieht 


Fig. 2. Fernordnung S als Funktion der Temperatur. 
BrAGG und WILLIAMS ..... 
BETHEE OD.‘ uno. 
KIRKWOOD 


‘daraus, daß unterhalb einer kritischen Temperatur T, jegliche Fernordnung 
‘ aufhört zu existieren, das Verhalten der spezifischen Wärme (Konfigu- 
rationsanteil) für den Fall AB, findet der Leser in Fig. 3; diese zeigt eine 


Fig. 3. Spezifische Wärme als Funktion der Temperatur. 
BrAGG und WILLIAMS ..... 
BE 
KIRKWoOOoD Tr 
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Unstetigkeitsstelle für die kritische Temperatur T... Es sei hier noch bemerkt, 
daß durch diese Kurven nur qualitativ das wirkliche Verhalten der spezifi- 
schen Wärme, wie sie etwa bei der Überstrukturbildung des AuCu, oder 
des -Messing gemessen werden, wiedergegeben wird. Experimentell wird 
immer noch über die kritische Temperatur hinaus ein anomaler Kon- 
figurationsanteil zu der spezifischen Wärme gefunden. Dieses Verhalten ist 
nicht verwunderlich, wenn man bedenkt, daß nach dem Zusammenbruch 
der Fernordnung im Kristallgitter sicher noch eine gewisse Nahordnung 
vorhanden sein wird. Das Verhalten der spezifischen Wärme für die BRAGG- 
Wirrıams-Theorie beruht eben auf der Tatsache, daß man ohne Begründung 
die Gitter in falsche oder richtige Plätze für die einzelnen Atomarten auf- 
teilt, ohne sich zu überzeugen, ob eine solche Einteilung einen Sinn hat. 
Ein weiterer Nachteil ist der Ansatz Gleichung (II, 2), der den statistischen 
Schwankungen der Gitterplatzbesetzung in keiner Weise Rechnung trägt. 
Aus diesen Gründen wurde bald nach den von BrAGG und WILLIAMS ver- 
öffentlichten Arbeiten eine neue Theorie von BETHE (12) entwickelt, welche 
die Hauptschwäche des BRAGG-WILLIAMS-Modells beseitigt, dafür aber neue 
in Kauf nehmen muß. 


II, 2. Dieexakte Lösung des eindimensionalen Modells mit Hilfe 
der Wechselwirkungsenergien der nächsten Nachbarn. 


Ehe wir jedoch dazu übergehen, die BETHE’schen Näherungsverfahren 
zu diskutieren, wollen wir am eindimensionalen Modell die exakte Lösungs- 
methode durchführen. Der Grundgedanke ist hier der, daß man die Wechsel- 
wirkungsenergien der nächsten, übernächsten usw. Nachbarn als bekannt 
voraussetzt und aus diesen Werten die Zustandssumme des Systems zu 
finden sucht. In vielen Fällen, z. B. bei metallischer Bindung, wird man nicht 
allzu große Fehler machen, wenn man nur die Wechselwirkungsenergien 
der nächsten Nachbarn berücksichtigt und die höheren Wechselwirkungen 
als vernachlässigbar betrachtet. Für dieses eindimensionale Modell ist der 
exakte Lösungsweg bereits von IsınG (13) auf kombinatorischem Wege für 
den (A, B)-Fall mit beliebigem Mengenverhältnis A : B angegeben worden. 
Seine Berechnungen beziehen sich zwar auf ein eindimensionales ferro- 
magnetisches Modell, aber im Prinzip gilt die Rechnung in der gleichen 
Weise für den Ordnungs-Unordnungsübergang. Der kombinatorische Weg 
Isıng’s ist folgender: Gehen wir von der geordneten Kette ABABAB...... 
aus; man kann nun in diese Kette aus N-Gliedern 1,2...... N—1 Fehler 
(zwei gleiche Atome nebeneinander) einbauen. Dabei ist zu beachten, daß 
sowohl O0 und N—1 Fehler je einen Ordnungszustand darstellen. Man be- 
rechnet nun mit Hilfe der Kombinatorik, auf wieviele Weisen man x Fehler 
auf die N-Glieder der Kette verteilen kann; jedem Fehler kommt ein un- 


günstiger BoLTzMmAnN-Faktor exp — V/kT, x Fehlern also exp — 27 zu. 


Die Summierung über 


Das Problem der Ordnungs—Unordnungsübergänge in Kristallen 105 


ergibt dann die Zustandssumme Z. Aus der Zustandssumme kann man nun 
alle wichtigen thermodynamischen Eigenschaften ableiten. Gemäß bekann- 
ten Sätzen aus der statistischen Thermodynamik ist die freie Energie 

(IL, 10) F=—kTlnZ. 

Wir werden diese Beziehung später noch mehrfach benutzen. 

Für das eindimensionale Modell hat sich die Differenzengleichungs- 
methode als besonders machtvoll erwiesen. Diese wurde von A.‘J. C. WILSON 
(14) erstmalig für das Röntgenbeugungsproblem angewendet und läßt sich, 
wie im folgenden gezeigt wird, sehr leicht zur Ermittlung der Zustands- 
summe verwenden. Zu diesem Zweck führen wir die Wechselwirkungs- 
energien zwischen gleichen und ungleichen Teilen ein. Vs (W. E. zwischen 
zwei A-Atomen), Vpg (W.E. zwischen zwei B-Atomen), Va = Vpa (W.E. 
zwischen A- und B-Atomen). Der Einfachheit wegen setzen wir Vaa = VB»: 
Den noch frei zu wählenden Nullpunkt der potentiellen Energie legen wir so, 
daß VB = — V4p: Es ist also VaB = (&s VpB — — ei) 

Schreiben wir noch n — — H, so können wir nun leicht 


kT 


nach folgendem Schema die Zustandssumme ermitteln. 
Ei B 


exp H \  exp—H exp He exp H 
A NB A7 B 
u. re a en are 
| 
A B AB u ls Aue 


Tabelle 2. Schema für die Aufstellung der Differenzgleichung zur Ermittlung der 
Zustandssumme. 


Wenn wir nun das Schema bis zum N-ten Glied ergänzt denken, so 
sieht man ohne weiteres, daß die insgesamt 2N—1-Endpunkte der beiden 
Schemata, die mit A oder B bezeichnet sind, die gesamten Konfigurationen 
enthalten; wenn man weiterhin die an den Verbindungsstrichen eingesetzten 
Multiplikationen mit den Exponentialfaktoren durchführt, erhält man die 
zu jeder Anordnung gehörigen BoLTzMmAnn-Faktoren. Wir nehmen nun an, 
uns sei das Ergebnis der (m — 1)-ten Zeile bekannt; Be repräsentiere die 
auf A endenden Glieder und DI die auf B endenden Glieder der Zeile. Wir 
können dann sofort daraus das Ergebnis der m-ten Zeile ableiten. Es ist 
nämlich (s. Tab. 2) 

Ps = Pe espH + Pl, exp— H, 
(IT, 11) Poren eHh.t Pb, .,eospH. 


Faßt man Ba Pe und pt Be als die Komponenten des zwei- 
dimensionalen Vektors pm und pm_ı auf, so sieht man unmittelbar, daß die 


*) Es entsteht nur eine geringe Komplizierung der Rechnung, wenn Va # Vp5 
gesetzt wird. 
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Gleichungen (ILL, 11), weiter nichts aind als eine Transformationsgleichung, 
mit deren Hilfe man die Komponenten des Vektors pm aus den Kom- 
ponenten des Vektors pm berechnen kann. Die Koeffizienten dieses Glei- 
chungssystems bilden eine Matrix H vom Rang 2. 
H— Be ann! 
exp— H expH 
Man kann deshalb für die Gleichungen (II, 11) auch wie folgt schreiben 
(IL, 12) Pm —= H Pm—ı 
und diese Gleichung nach der Matritzenmethode ausrechnen. Wir werden 
auf dieses Verfahren an späterer Stelle zurückkommen und uns hier dem 
mathematisch einfacheren Verfahren der Rekursionsgleichungen zuwenden, 
die sich lösen lassen. Man kann nämlich durch Bildung von Pen PS die 
beiden Gleichungen zunächst auf zwei nur P, bzw. P, enthaltende Glieder 
reduzieren. Sie lauten 
(1,13) Pa — Pu, 2expH—P, s(esp  2H- ep2H) 0 
pP, Pro 2pH pP (exep— 2H—exp2H) —=0. 
Man sieht, daß die Lösung beider auf die gleiche charakteristische Gleichung 
zurückgeführt werden kann. Man erhält diese durch den Ansatz Pı=P2 
— x®, Das ergibt die Gleichung 
x®22[x?— (2expH)x +exp2H—exp—2H] =0. 

Damit ist die Lösung auf die in der eckigen Klammer stehende cha- 
rakteristische quadratische Gleichung zurückgeführt worden. Die beiden 
Lösungen dieser Gleichung lauten 


(IL 14) x, =2c0shH, 
x; =i2enhH. 
Die allgemeine Lösung unserer Gleichung (II, 11) ist dann 
PB, Aox ar Bir,s: Be Az u Bern 
Die Konstanten lassen sich aus den Randbedingungen ermitteln: 
Schema 1 Schema 2 Schema 1 Schema 2 
Perl By 0 Por =0 | 
Pa expH PA =exp—H Pr =ıp—H P>eop 


Die Angabe zweier Konstanten entspricht den beiden Schemata unserer 
Tabelle 2. Man sieht aus diesen Bedingungen, daß das Ergebnis für die 
beiden Schemata der Tabelle 2 identisch ist (nur Vertauschung von A und B). 
Wir berechnen daher die Konstanten nur für Schema 1 und erhalten 


AA: Bi =}: Au: BB — — 


wm 


Daraus ergibt sich die Zustandssumme: 
(II, 15) Z=2(Ps_ı + Pi) = 2x1 =2(2 cosh H)N—ı. 

Mit Hilfe der Gleichung (II, 15) kann man alle interessierenden Daten 
errechnen. Wir wollen das Ergebnis hier nicht näher erläutern. Es sollte 


damit nur die Methode am eindimensionalen Modell dargestellt werden. 
Die Beschränkung auf die Wechselwirkung nur der nächsten Nachbbarn 
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bietet hier keine Schwierigkeit. Der Verfasser hat an einem speziellen Bei- 
spiel (15) die Berechnung bis zu Wechselwirkungen der drittnächsten Nach- 
barn erweitert. Allerdings wurde dabei nicht wie hier mit den relativen, 


- sondern mit den absoluten Wahrscheinlichkeiten zur Berechnung des Rönt- 


genbeugungsproblems gearbeitet. Grundsätzlich lassen sich die Berech- 
nungen in analoger Weise auch auf die Zustandssumme erweitern. 


Die Anwendung der Differenzenmethode auf das zwei- und dreidimen- 
sionale Fehlordnungsproblem bringt große Schwierigkeiten mit sich, die 
nicht in der Aufstellung der Gleichungen, sondern in ihrem hohen Grad 


‚ liegen. Diese führen dann auf Rekursionsgleichungen sehr hohen Grades, 


deren Lösung praktisch unmöglich ist. Zwar benötigt man nur die größte 
Wurzel zur Berechnung der Zustandssumme; denn da das Endergebnis in der 
Form Z= X A,x‘ gegeben wird, ist für große N für die Zustandssumme 
nur das Glied mit xmax von Bedeutung. Eine Näherungslösung nach der 
Differenzenmethode wurde von F. ZERNIKE (16) für ein dreidimensionales 
Modell durchgeführt, die er dadurch erreichte, daß er die relativen Be- 
setzungswahrscheinlichkeiten der nächsten Nachbarn unabhängig vonein- 
ander ansetzte, eine Verfahren, das vor allen Dingen für das Übergangs- 


‚ gebiet zwischen Ordnung und Unordnung bedenklich erscheint. 


II, 3. Das BETHE’sche Näherungsverfahren. 


BETHE (12) hat nun, analog zu den untenstehenden Berechnungen am 
eindimensionalen Modell, unter alleiniger Berücksichtigung der Wechsel- 
wirkungsenergien der nächsten Nachbarn einen Ansatz durchgeführt, mit 
dessen Hilfe er über zwei Näherungswege versucht, der exakten Lösung 
nahe zu kommen. Der grundlegende Unterschied gegenüber der Methode 


' von BrAacG und WiLrıams liegt darin, daß man die a-priori-Aufteilung in 


„falsche“ und „richtige“ Gitterplätze für die einzelnen Atomarten vermeidet. 
Hier wird also nicht allein der BracG-WiıLrıams’sche Fernordnungspara- 
meter S verwendet, sondern eine weitere Definition für den Ordnungs- 
zustand eingeführt. 


Wir betrachten ein Gitter, bei dem jeder Gitterplatz von z-nächsten 
Nachbarplätzen umgeben ist. Die Gesamtzahl der nächsten Nachbar- 


r 


: IN 
bindungen beträgt, wenn insgesamt N-Gitterplätze vorhanden sind, —.-- 


Nennen wir für eine bestimmte Anordnung die Wahrscheinlichkeit, daß 
irgendein herausgegriffenes Paar nächster Nachbarn ein A—B-Paar ist, q, 
und bezeichnen wir weiterhin die gleiche Wahrscheinlichkeit für den geord- 
neten Zustand mit q, und für den völlig ungeordneten Zustand mit qu, SO 
lautet die BETHE’sche Definition des Nahordnungsparameters o 
(IT, 16) Se 
Io — u 
Man sieht daraus, daß 
0 02 füre=g = gu 
efinng=dg, ist. 
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Man erkennt den Unterschied der Definition von S und o am besten 
bei Betrachtung der Tabelle 3. 
ABABBA 
BABAAB 
ASBEASBEBER 
BEA®BEATATBSAEB 
Tabelle 3. Erläuterung der Verschiedenheit des Fernordnungsparameters S und des 
Nahordnungsparameters o an einem extremen Beispiel. 


. In Tabelle 3 ist der Fernordnungsparameter S = 0, während o = 0,8 
ist (Auswertung durch Abzählen der Bindungen). Schon rein gefühlsmäßig 
möchte man dem letzten Wert den Vorzug geben, da das angeführte Beispiel 
doch einen erheblichen Ordnungsgrad besitzt. 

Man kann nun wie vorhin die Kristallenergie als Funktion des Ord- 
nungszustandes ausrechnen. Wir wollen die Berechnung hier nur für den 
A—B-Fall (A:B=1; z. B. Cs—Cl- bzw. Na—Cl-Typ) durchführen. Für 
diesen Fall vereinfacht sich, da q,u = 1 und qu = # wird, Gleichung (II, 16) 
zu 
(11,17) == 

bzw. q=4(1+o0), 1—q=4(1— 0). 


Führen wir wieder, wie bisher, die Wechselwirkungsenergien Vaa, Vs 
und Var ein, so ist E 


Nz V V 
E= [avi + aa al. 
Einsetzen von o aus Gleichung (IL, 17) ergibt 
Nz V V 
(IT, 18) = Tut Ya+l gut tel, 
Nz & V V N 
= 12 Van zu ( _ = = Van)); 


Var + VeB 
5) 


u 


Nz : 

— const Ey V(1l—-o) mit V= — Vıp- 

Der Unterschied zwischen der BRAGG-WıLLıams’schen und der BETHE- 
schen Berechnung tritt hier klar zutage. Hier werden unmittelbar die 
Bindungen abgezählt, es wird also auch den Schwankungen im kleinen 
Rechnung getragen, während bei BRAGG-WILLIAMS die Energie durch den 
makroskopischen Ordnungszustand S mittels des stark vereinfachten An- 
satzes V = V,S gewonnen wurde. | 

Es muß nun wieder, wie es auch im BrAGG-Wirrıams-Verfahren getan 
wurde, der Nahordnungsgrad o als Funktion der Temperatur gefunden 
werden. Aber der Vergleich mit der BraGG-WiLLIams-Theorie zeigt schon, 
daß o allein zur Beschreibung des Ordnungzustandes nicht ausreicht. Man 
kann das leicht einsehen, wenn man sich zwei Fälle gleichen Nahordnungs- 
grades (gleiche Zahl von AA-, BB- und AB-Bindungen) einmal mit Fern- | 


Das Problem der Ordnungs—Unordnungsübergänge in Kristallen 109 


ordnung und ein andermal ohne Fernordnung (wie z. B. in Tab. 3) kon- 
struiert; beide werden dann durch dasselbe co gekennzeichnet. Die Lösung 


“ des Problems hat sich also in erster Linie damit zu befassen, eine Brücke 


zwischen S und o zu schlagen. Zwar bestimmt o allein, ob auch Fernordnung 
vorhanden ist oder nicht, aber man kann dem Parameter o nicht unmittel- 
bar ansehen, ob er Fernordnung bedingt oder nicht. Dieser Zusammenhang 
von S und o, wobei also o durch die Temperatur festgelegt wird, ist eines 
der Kernprobleme, mit deren Lösung wir uns im folgenden zu beschäftigen 
haben. 


Wir wollen zunächst das BETHE’sche Verfahren am zweidimensionalen 
Problem des quadratischen Netzes erläutern (A—B-Fall). In Figur 4 ist 
ein quadratisches Netz mit einigen Gitterplätzen aufgezeichnet. Ausgehend 


Fig. 4. Einteilung der Schalen um ein Zentralatom. 
Die Doppelpfeile kennzeichnen die Bindungen der nächsten Nachbarn. 
+ Zentralatom 
@ |. Schale 
© 2. Schale 
& 3. Schale 


vom Zentralatom werden die benachbarten Gitterplätze als 1. Schale, die 
übernächsten als 2. Schale usw. bezeichnet. Der Gesamtkristall besteht aus 
sehr vielen Schalen. Gefragt wird nach der Wahsrscheinlichkeit. der Be- 
setzung des zentralen Gitterplatzes (beim unendlichen Gitter ist aber jeder 
Gitterplatz Zentralplatz) in Abhängigkeit von seiner Umgebung. Diese 
Wahrscheinlichkeit ließe sich sofort angeben, wenn wir die Wahrscheinlich- 
keit der Besetzung der benachbarten Gitterplätze kennen würden. Da diese 
aber wieder von den Wahrscheinlichkeiten der zweiten, diese wieder von 
denen der dritten usw. abhängen, ist die Berechnung der Besetzungs- 
wahrscheinlichkeit des Zentralatoms sehr schwer anzugeben. Prinzipiell 
läßt sich natürlich die Aufgabe lösen, wenn die Besetzung der Randatome 
gegeben ist. Aus dieser kann man dann sukzessive die Besetzung der inneren 
Schalen berechnen. Man kann das Problem auch anders betrachten. Gegeben 
ist die Besetzung des Zentralatoms, man berechne die Wahrscheinlichkeiten 
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der Besetzung der ersten, zweiten, dritten usw. Schale, die also wieder, wie 
im eindimensionalen Fall, auf die Aufstellung von Differenzengleichungen 
zurückgehen. Die Frage, ob Fernordnung vorhanden ist oder nicht, läßt 
sich in beiden Auffassungen leicht präzisieren. Hat nämlich die Rand- 
besetzung noch einen Einfluß auf die Besetzung des Zentralplatzes bzw. hat 
der Zentralplatz noch einen Einfluß auf die Besetzung des Randes, so ist 
Fernordnung vorhanden. In diesem Fall hat dann die Aufteilung in A- und 
B-Plätze von BRAGG und WILLIAMS wieder einen Sinn, weil diese Unter- 
scheidung durch die Randbesetzung bzw. die Zentralbesetzung vorgegeben 
wird. Diesem Gedanken trägt BETHE dadurch Rechnung, daß er die Um- 
gebung des Zentralatoms in zwei Teile teilt. In der ersten Näherung 
bilden das Zentralatom und die ersten Nachbarn den Innenraum, alle 
übrigen den Außenraum. In der zweiten Näherung werden noch die 
zweiten Nachbarn in den Innenraum einbezogen. Den Außenraum kenn- 
zeichnet BETHE so, daß er in der Lage ist, die Plätze des Innenraumes in 
„richtige“ und ‚‚falsche‘‘ einzuteilen. Man sieht also, daß auch die BETHE- 
sche Theorie nicht ganz von der Hypothese der Einteilung der Gitterplätze 
loskommt, aber diese Einteilung ist nun nicht mehr an ein makroskopisches 
System von Gitterplätzen gebunden; sie werden vielmehr von den lokalen 
Schwankungen bestimmt. Wir wollen jetzt das BETHE’sche Verfahren für 
die erste Näherung diskutieren (first shell approximation). Die zweite, die 
der ersten vollkommen analog ist, wollen wir aus Raumgründen nicht mit 
heranziehen. Es sei nur bemerkt, daß die zweite Näherung keine sehr großen 
Verbesserungen mit sich bringt, so daß man wohl annehmen darf, daß schon 
die erste Näherung eine ansprechende Lösung darstellt. Allerdings muß 
hier bemerkt werden, daß die Konvergenz der beiden Näherungen noch kein 
absoluter Beweis für die Richtigkeit ist; sie gelten sicher in den beiden Grenz- 
fällen Ordnung und Unordnung, im Übergangsgebiet der Fernordnung ist 
die Frage aber sehr strittig. Mit der Erweiterung des Innenraumes werden 
die Ergebnisse wahrscheinlich genauer, weil ein größerer Prozentsatz der 
Umgebung zur exakten Berechnung herangezogen wird. Die Einteilung des 
Innenraumes in falsche und richtige Gitterplätze wird aber gefährlicher, da 
einerseits lokale Schwankungen im Innenraum stören können und der 
Außenraum selbst in seiner Oberfläche durch lokale Schwankungen eine 
solche Definition sehr schwierig macht. 

In der ersten Näherung führt BETHE für jedes ‚‚falsche‘‘ Atom in der 
ersten Schale eine mittlere ungünstige Energie U ein, die also in der Zustands- 
summe als BOLTZMANN-Faktor 

e = exp — U/kT 

in Erscheinung tritt. Diesen gemeinsamen Faktor für jedes falsche Atom 
kann BETHE durch die numerische Durchrechnung der Besetzungswahr- 
scheinlichkeit des Zentralatoms im quadratischen Gitter mit der 1. Schale 
allein und der 1. und 2. Schale plausibel machen, indem er zeigt, daß für 
jedes falsche Atom in der ersten Schale bei Berücksichtigung der 2. Schale 
ein Faktor auftritt, der also den Wechselwirkungen der ersten mit der 
zweiten Schale Rechnung trägt. 

Wir haben jetzt ein beliebiges Gitter der Koordinationszahl z vor uns. 
Das statistisch mechanische Gewicht von n falschen Atomen in der ersten 
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Schale ist gleich der Anzahl der Permutationen der z Plätze dividiert durch 
die Permutationen der n falschen Atome und (z—n) richtigen Atome unter- 


einander, also 
AN! zZ 
n!(e—n)! \ny 


Ist der Zentralplatz ‚richtig‘ besetzt, so erhält jedes falsche Atom in der 
1. Schale den BoLrzmann-Faktor x = exp — V/kT, außerdem erhält es den 
schon erwähnten Faktor e = exp — U/kT durch den Außenraum. 


Die relative Wahrscheinlichkeit, daß der Zentralplatz ‚richtig‘ besetzt 
ist, wenn n Plätze der ersten Schale ‚falsche‘ Atome haben, ist demgemäß 


Z 
= xrel 


Die relative Wahrscheinlichkeit r, daß das Zentralatom „richtig“ ist, 
beträgt dann 


(II, 19) u > nn > (4) esar 
n=0 


Auf analoger Weise kann die Wahrscheinlichkeit, daß der Zentralplatz 
„falsch‘‘ besetzt ist, ausgerechnet werden. 


(IT, 20) = 5" en (et x), 
2=0 


Die relative Wahrscheinlichkeit p„, daß in der 1. Schale n ‚falsche‘“ 
Atome sitzen, ist aber r, (f„), da ja r„ (fn) die Wahrscheinlichkeit ist, daß 
n-,‚falsche‘‘ Atome in der ersten Schale waren und zugleich der Zentralplatz 
„richtig“ (‚‚falsch‘‘) besetzt war. Damit kann man die Wahrscheinlichkeit f’ 
berechnen, daß ein Atom in der 1. Schale ein ‚‚falsches‘“ ist. 


(II, 21) 


1 - 1 EM EIX € 
I — ] = — — & —— ı 
f z > N Pn z > n (In +) rn f) lex a: 
n=0 


Da im unendlichen Kristall die physikalische Bedeutung von f und f’ 
gleich ist, setzt BETHE f = f’. Man erhält also aus (II, 20) und Il, 21) 


| ie) 


u ee I ex 
oder 1 
e+X Se erg 
1 ER 22 
2 il 9 
Setzt man 
sep 292 —|), 
so wird 
sinh (z — 2) ö 
5) Fe an, Ey: 
B1,22) z sinh z ö 
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Aus Gleichung (II 22) läßt sich e als Funktion der Temperatur ermit- 
teln. Aus den Eigenschaften der sinh-Funktion sieht m&an sofort, daß 
(x = 0/0) für alle Temperaturen eine Lösung darstellt. Ist ö eine Lösung 
von (II, 22), so ist auch — 6, (.) eine Lösung (Austausch von „richtigen“ 

1 
und ‚‚falschen‘‘ Atomen). Das Wesentliche von Gleichung (II, 22) ist aber, 
wenn überhaupt eine von 0 verschiedene Lösung für ö existiert, x niemals 
größer als 1 — 2/z (Grenzwert von x für ö = 0) werden kann. Istx > 1— 2/z, 
so ist nur noch die Lösung ö = 0 möglich, d. h. es ist keine Fernordnung 
vorhanden. 

Mit Hilfe der hier berechneten Beziehung lassen sich die anderen Zu- 
sammenhänge in verhältnismäßig einfacher Weise ableiten. Der wesentliche 
prinzipielle Unterschied zwischen der Berechnung von BRAGG und WILLIAMS 
und der von BETHE liegt im Verlauf der spezifischen Wärme (Fig. 3), die 
noch über die kritische Temperatur hinaus einen anomalen Verlauf zeigt, 
der von dem für die Nahordnung notwendigen Anteil der Energie herrührt. 
Die theoretische Deutung der experimentellen Ergebnisse ist damit verbessert 
worden. Wir wollen hier die anderen Zusammenhänge nicht explizit aus- 
rechnen, weil dabei keine grundsätzlich neuen Gedanken auftreten. Der 
näher interessierte Leser sei auf die Originalarbeit von BETHE (12) ver- 
wiesen. 

Die zweite Näherung wird in ganz analoger Weise durchgeführt, nur 
ist sie mathematisch etwas schwieriger. Die Ergebnisse weichen aber, wie 
vorhin schon erwähnt, nicht wesentlich von der 1. Näherung ab. BETHE 
hat seine Berechnungen nur auf den A—B-Fall erstreckt, für den Fall AuCu, 
sind diese Ergebnisse nicht mehr gültig, weil ja bei der Berechnung an- 
genommen wurde, daß die Bindungen der nächsten Nachbarn nur zwischen, 
aber nicht innerhalb der Schalen liegen. Für die dichtesten Kugelpackungen 
ist dieses Verfahren aber nicht anwendbar, da ja in diesem Fall die eben 
erwähnte Bedingung nicht erfüllt ist. Für diesen Fall (AuCu,) ist das 
BETHE’sche Verfahren von PEIERLS (17) erweitert worden. Die Ansätze 
folgen der gleichen Linie, wie es hier bei der ersten Näherung dargestellt 
wurde. Allerdings treten erhebliche mathematische Komplikationen dabei 
auf. — Auch die Beschränkung auf Wechselwirkungen der nächsten Nach- 
barn bedeutet für dieses Verfahren keine grundsätzliche Schwierigkeit. Das 
gleiche gilt auch für das strenge stöchiometrische Verhältnis. Über Er- 
weiterungen der BEerHe’schen Theorie in diesen Richtungen, die hier aus 
Raumgründen nicht diskutiert werden können, sei der Leser auf Arbeiten 
von EASTHOPE (18), SHOCKLEY (19) und CHangG (20) (21) verwiesen. Es sei 
nur erwähnt, daß bei diesen Erweiterungen keine wesentlich neuen Züge 
in das Verhalten der thermodynamischen Funktionen gebracht werden. Die 

Änderung der kritischen Temperatur in Abhäng igkeit vom Mischungs- 
ver. hältnis ist in Figur 5 für den AB-Fall w iedergegeben. 


LI, 4. Kırkwoop’s Semi-Invariantenmethode. 
KIRKwooD (22) hat eine elegante mathematische Methode entwickelt, 


die besser v erallgemeiner ungsfähig ist als das BETHE’sche Verfahren. In der 
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BraGG-WirLrams-Theorie wurden jedem Fernordnungsgrad S gemäß Formel 
(IL, 2) eine bestimmte Energie zugeordnet. Wie wir aber bei der Behandlung 
des Fehlordnungsproblems nach der BETHE’schen Methode gesehen haben, 
liegt in dieser Zuordnung eine große Schwäche. Jedem S sind nach dem 
- Wechselwirkungsprinzip eine Fülle von Konfigurationen mit verschiedenen 
Energiewerten zuzuschreiben. Die Energie als Funktion von S kann also 
bei BraGG-WILLIamS nur den Charakter eines Mittelwertes tragen. Die 
direkte Berechnung der Energieverteilungskurve für ein bestimmtes S wurde 
noch nicht durchgeführt. Kırkwoop’s Methode beruht auf Reihenentwick- 
lungen und bringt keine grundlegend neuen Gedanken in das Problem hinein; 
- sein Lösungsverfahren ist deshalb mehr vom mathematischen Standpunkt 
aus interessant. Wir werden den Gedankengang in verschiedenen Punkten 
nicht explizit wiedergeben, sondern nur das Grundsätzliche des Verfahrens 
so klar wie möglich herauszustellen suchen. 


Wir bezeichnen mit W die Anzahl aller möglichen Konfigurationen des 
Kristalls und mit E; die Energie des Zustandes W;. Die Zustandssumme ist 


dann 
w 


(IT, 23) ve > exp — Ey/KT. 
i=1 
Aus ihr läßt sich die freie Energie F berechnen, es ist 
ne 
(1,24) F=—kTinZ bzw exp— FkT=Z- > exp — Ey/KT, 


i=1 


| 


ö 1l 
dafür schreiben wir mit x = — KT 
IT, 25) expxF = > exp xE,;, 
| i=1 
' F ist dabei eine Funktion der Variablen x und E;,. Für die folgenden Über- 

legungen nehmen wir wieder an, daß die E; temperaturunabhängig sind, 
' so daß die E; für ein bestimmtes Problem als konstante Größen angesehen 
werden können. Man kann dann die Funktion x F in eine Reihe entwickeln. 
' In dieser Reihenentwicklung treten die an sich unbekannten ersten zweiten 
usw. Ableitungen der Funktion xF an der Stelle x = 0 auf, die wir mit A, 
‘“ (nullte Ableitung), A, (erste Ableitung), Aa, -- ----- bezeichnen; dieses sind 
_ die sogenannten THIELE’schen Semi-Invarianten 4;, die natürlich von den 
E; abhängig sind. Wir führen also für die linke Seite von Gleichung (II, 25) 
die Entwicklung durch, während wir die Exponentialglieder der rechten 
Seite in ihrer Reihendarstellung aufschreiben. Damit erhalten wir 


e x? x 
exp urahtghtadrt... 


wi BE „a w 
Wis D)E+ DB t+n D,Eit 
i=1 Il en) 


Fortschritte der Mineralogie 1949. Bd. 28. 8 


D| 4 
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Da aber 
w Ww w 
Du-WENm-WE;) Er-Wwe TEN, 
al 2 i=1 


erhalten wir 


(II, 26) exp Bra rapie ap = 


=wi4sE4, B+5, en. 


d.h. expf,(x)=1,(x); 
differenziert ergibt das 
at, X) dEK&) _ at,(2) 


(127) exp f,(x) ech: 6) 


Da die Beziehung (II, 27) für alle x gilt, müssen die Koeffizienten der 
beiden Reihenentwicklungen gleich sein. Setzen wir die durch Gleichung 


(II, 26) definierten Beziehungen für f,(x) und f,(x) in Gleichung (II, 27)' 


ein, so erhalten wir 


- — X 00 xd— 
(IT, 28) wWiE+sB+4,W+ JEi+.... 


x3 


F  < en ale x? ” 
-W|I4xE4,P+ SB + en It 


3 
A, ist der Wert für xF an der Stelle x = 0. 


Wir können ihn aus Gleichung (II, 25) ableiten. Weil alle Exponential- 
faktoren 1 werden, ist xF = In W, also A, = In W. Die anderen Werte für 
DNS bestimmen wir durch Koeffizientenvergleich der ausmulti- 
plizierten rechten Seite mit der linken Seite von Gleichung (II, 28) und er- 
halten 


el me nn +34,4,)—3 (EP — Bi = 4, — 34% 


Damit ist das Problem der Auffindung der Semi-Invarianten auf die 
Bestimmung der verschiedenen Mittelwerte der Energie zurückgeführt wor- 
den. Sind diese Mittelwerte bekannt, so kann man aus (11,29) die Semi-Inva- 
rıanten ausrechnen. Damit erhält man die freie Energie aus Gleichung (II, 26) 
A, E As Ay — 34,5? 


I, 30) Ver FR = 
ar IE 2IkKT "3IKT® AKT) 


Das Kernproblem ist, wie schon erwähnt wurde, die Berechnung der 44 


Dies ist eine Aufgabe der statistischen Mechanik und kann je nach dem vor- 


| 
| 
| 
| 
| 
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liegenden Gittermodell gelöst werden. Allerdings ist die Berechnung für 
die höheren A sehr mühsam. Von Kırkwooo ist sie für den AB-Fall bis Ay 
ausgeführt worden, eine Erweiterung auf }, und A, wurde von CHAnG (23) 
vorgenommen. Wir wollen hier die Berechnung nicht ausführen, sondern 
nur die Ergebnisse der ersten Werte angeben. 


E NzV 

A,=E=-E,(1—8), E, = r (z = Koordinationszahl), 
PU ge 

ER .NaV (1 — 82)2, en 


Durch Einsetzen in (II, 30) erhält man die freie Energie 
NzV? (1 — 82 


F=—kTnW+B,(1—8) on 


Da In W bereits in Gleichung (II, 7) gegeben ist, erhalten wir fürF, =F,—=} 
(IL, 31) e 
Er E18) NKT|n2- 3049 ml+9-3A—9ma-8) 


E2(1— 82% 
NzkT | 
Setzt man mit KıRkwooD 
e zV en Ey 
Tat Ne 
und ermittelt das thermodynamische Gleichgewicht aus 
oF 
Be) 
os : 
so erhält man schließlich 
R& 5 2 
(II, 32) S=tgh a—— (1 zu], 


Für Gleichung (II, 32) kann man nun ähnliche Überlegungen ausführen, 
wie wir es mit Gleichung (Il, 22) für das BErTHE’sche Näherungsverfahren 
getan haben. Man erhält im Prinzip die gleichen Ergebnisse wie bei BETHE, 
also S = 0 oberhalb einer kritischen Temperatur T.. Man kann nun auch 
Formeln für die spezifische Wärme und die Energie ableiten; wir wollen 
jedoch hier nicht darauf eingehen. In den Figuren 2 und 3 sind die Ergebnisse 
der KIRKwooD’schen Theorie im Vergleich zu den Ergebnissen von BRAGG- 
WiırLıams und BETHE graphisch dargestellt. Sie sind ein guter Maßstab für 
die Güte der bisher behandelten Verfahren. 


Eine sehr wichtige Darstellung für die ‚Unordnungsübergänge ist die 
Abhängigkeit der Entropie von der Energie, wie die folgende Überlegung 
beweist. Die freie Energie ist definiert durch 

F=E-To» ® — Entropie des Systems. 
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Für das thermodynamische Gleichgewicht bei konstanter Temperatur 
gilt dF=dE—-Td®=0 
dE 
do 
d. h. die Tangente mit der Neigung T an die Energie-Entropie-Kurve gibt 
das Minimum der freien Energie, also die. wahrscheinlichste Energie des 
Systems bei der Temperatur T an. Man kann diese Beziehung auch un- 


mittelbar aus der graphischen Darstellung der Energie—Entropie-Kurve 
ablesen (s. Fig. 6). Es wurde auf der Abzisse ® und auf der Ordinate E 


bzw. = |} 


Energie 


$ Entropie 


Fig. 6. Energie-Entropiedarstellung des Fehlordnungsproblems. 


aufgetragen. Im Dreieck EBE’ entspricht die Neigung E’B der Temperatur 
T des Systems. Das System befindet sich in einem Zustand, der durch die 
Werte E und ® für Energie und Entropie bestimmt ist. Da nun aber 
EE' 
gwpo=T=-— 


ist also die Länge der Strecke ER’ —= T®@. Die freie Energie F ist also offenbar 
durch E— T® —0E’ gegeben. Damit sieht man unmittelbar, daß die 
Tangente der Neigung T an die Energie-Entropie-Kurve das Minimum der 
freien Energie festlegt. 


Kırkwoon hat nun den Verlauf der Energie-Entropie-Abhängigkeit 
explizit ausgerechnet, und da mir dieser Zusammenhang sehr wichtig er- 
scheint, sei der Gedankengang der Rechnung hier kurz ausgeführt. Wir 
wählen dazu den noch freien 0-Punkt der Energie so,daß E = 0 ist. Die neue 
Energievariable U ist dann 
(IL, 33) U=E-.E. 
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Entsprechend wird die freie Energie 


(II, 34) F=F-—E. 

Wir entwickeln die Entropie in eine Reihe nach U und erhalten 
(II, 35) (U) =-&,+9 U a an 
Entsprechend entwickeln wir die Funktion y \o) =. nach 1/T, das ergibt 
au, 30) u ee 2 


Die in (II, 30) gewonnene Reihendarstellung für die freie Energie F ver- 
mittelt unter Berücksichtigung von Gleichung (LI, 34) eine einfache Be- 
ziehung zwischen den y, und den A, bzw. A, (Gleichung (II, 29)). Es ist 


Y=—khW, 
yvı=), 
(11737) Yw=— 4;lk, 
vw; — As/k’, 
—A,+ 342 
MT 5 k3 — 
Da nun 
ar 1% 


haben wir damit eine Beziehung zwischen F’ und @. 


Eine zweite Beziehung, die für unsere Zwecke brauchbar ist, gewinnt 
man aus (II, 38) durch Differentiation nach 1/T und erhält 
9 (F’/T) 
am 
Aus Gleichung (II, 38) findet man unter Berücksichtigung von (IL, 39) 
IIo(RöT 
Be uuTg Er 
Berücksichtigt man weiter die Reihendarstellungen (II, 35) und (II, 36) 
sowie die Gleichung (II, 39), so erhalten wir 
(II, 41) 


® ® 
B,+ 8, (m+7 at...) l)+3(-)+ 


1 Yy UF Yı Yı %Ys Ws 


(IL, 39) 


(IL, 40) 
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Durch Koeffizientenvergleich der beiden Seiten von Gleichung (II, 41) und 
Einsetzen der Ergebnisse für die y, aus Gleichung (II, 37) erhalten wir die 
ersten Werte für die ®, 


9, =—-Y=klaW, 
9%, =yı=0, 
1 k 
(II, 42) re 
= Yı=155, 


Damit gewinnen wir auch die ersten Glieder der Reihenentwicklung (LI, 35) 
der Entropie als Funktion der Energie unter Berücksichtigung von 
U=E-E 


1 > A — 
zn ne 
2 
Man = 
+ er 
. 2 


In ähnlicher Weise kann man auch die spezifische Wärme und die 
anderen interessierenden thermodynamischen Eigenschaften ableiten. Die 
Genauigkeit der Kırkwoop’schen Methode hängt davon ab, wie gut die 
Reihen konvergieren. Man kann natürlich die Näherungen beliebig ver- 
bessern, wenn man die Reihenentwicklungen sehr weit treibt, jedoch ist 
die Ausrechnung der höheren Glieder sehr mühsam, sie hängt auch von dem 
gewählten Gittertyp und der Reichweite der eingeführten Wechselwirkungen 
ab. Von KIRKwooD und CHANG wurden die Semiinvarianten bis A, berechnet; 
dieses Ergebnis dürfte schon eine recht gute Näherung darstellen. 


Wir wollen hier auf eine weitere ausführliche graphische Darstellung 
verzichten und nur noch darauf hinweisen, daß auch die KıIRkwooD’sche 
Theorie das eigentliche Problem des Eintritts der Fernordnung nicht direkt 
als Ergebnis der Berechnungen enthält. Diese wurde durch die Einführung 
des Fernordnungsparameters S vorgegeben, der Vorteil gegenüber der 
BrA6G-WiırLıams’schen Methode liegt hier aber in der Berechnung der 
Energie nach dem Wechselwirkungsprinzip der Nachbarn. KIRKWOOD’s 
Methode ist daher als Bindeglied zwischen BrAGG-WILLıAMS und BETHE 
zu betrachten, wobei hier der Vorteil besteht, daß die aufgestellten Formeln 
allgemein gültig sind und nur die Bestimmung der TuıELe’schen Semi- 
Invarianten eine Frage der gewählten Bedingungen ist. Die etwas bedenk- 
liche Methode der Einführung eines Einteilungsprinzipes im Innenraum in 
falsche und richtige Gitterplätze fällt hier fort. 
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II, 5. Die quasichemische Methode. 


Wir wollen uns nun mit der sogenannten quasichemischen Methode 
befassen, die sich in ihrer Anwendbarkeit als sehr bedeutungsvoll erwiesen 
hat. Der Grund dafür liegt offenbar in ihrer mathematischen Einfachheit 
der Darstellung des physikalischen Sachverhalts; deshalb scheint mir diese 
Methode besonders gut auf kompliziertere Probleme anwendbar zu sein. 
Eine ihrer Schwächen liegt darin, daß an die Spitze aller Berechnungen eine 
Hypothese gestellt werden muß, welche die Gleichgewichtsanzahlen der AA-, 
BB- und AB-Bindungen nach der Art des bekannten Massenwirkungs- 
gesetzes der Chemie festlegt, und deshalb quasichemische Methode genannt 
wird. Sie geht zurück auf eine statistisch mechanische Behandlung der 
sogenannten regulären Lösungen, die von GUGGENHEIM (24) entwickelt, von 
RUSHBROORE (25) korrigiert und verallgemeinert und schließlich von FOWLER 
und GUGGENHEIM (26) auch auf das Fernordnungsproblem in Kristallen 
anwendbar gemacht wurde. 


Um ein tieferes Verständnis für die Lösungsmethode zu erhalten, wollen 
wir zunächst ganz kurz auf die regulären Lösungen eingehen. Es handelt 
sich hierbei um Lösungen, die folgende Eigenschaften besitzen: 

1. Zwischen den Teilchen (Molekülen und Atomen) der Lösung sind die 
Wechselwirkungsenergien zwischen weit entfernten Teilchen vernach- 
lässigbar (Nächste-Nachbarnprinzip von BETHE). Die Energie wird 
also wie bisher aus den Wechselwirkungsenergien der nächsten Nach- 
barn allein berechnet. 

. Die Teilchen können als Kugeln betrachtet werden und haben an- 
nähernd die gleiche Größe (auch hier besteht ein Analogon zu unserer 
alten Forderung des geometrisch ungestörten Mischkristalls). 

3. Jedes Teilchen ist von z anderen Nachbarn umgeben; die Koordina- 
tionszahl z ist unabhängig von der Art der Teilchen. (Sehr kristalliner 
Charakter der Lösungen, aber trotzdem scheint er berechtigt zu 
sein, da sich für alle Probleme die Einführung einer mittleren 
Koordinationszahl als sehr nützlich erwiesen hat.) 

4. Jedes beliebige Mischungsverhältnis bedingt keine Änderung des 
Volumens der gesamten Lösung (siehe Bedingung 2). 

Die Verwandtschaft der oben aufgestellten vier Forderungen mit dem 
Kristallproblem ist evident und man sollte von vornherein erwarten, daß 
die Ergebnisse ebenso gut auf Kristalle anwendbar sind. 

Wie bei den meisten Methoden gehen wir hier wieder aus von der Zu- 
standssumme Z, mit deren Hilfe wir die freie Energie aus der schon oft 
benutzten Gleichung 


—D 


— F/kT=1InZ 
berechnen können. Die Anteile der kinetischen Energie (Atom- bzw. Molekül- 
schwingungen) werden, wie bisher immer stillschweigend angenommen wurde, 
als konfigurationsunabhängig betrachtet und aus der Zustandssumme 
herausgezogen. Wir müssen also wieder den Konfigurationsanteil einer 
näheren Untersuchung unterziehen. Dieser entspricht nun aber einer Inte- 
gration der mit ihrem zugeordneten BoLtzMANnN-Faktor exp — E/kT multi- 
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plizierten Volumenelemente des Konfigurationsraumes, der bei N-Teilchen 
3N-dimensional ist. Jedem Punkt im Konfigurationsraum entspricht eine 
bestimmte Energie E;. Wir nehmen an, daß N, A-Teilchen und N, B-Teil- 
chen in der Lösung vorhanden sind. Man kann dann leicht einsehen, daß 
die Energie irgendeiner Konfiguration bereits durch Angabe der Anzahl 
der A—A-Bindungen gegeben ist, wenn die geforderten Bedingungen 1—4 
erfüllt sind. Ist nämlich jedes Teilchen von z Nachbarn umgeben, so gibt es 
insgesamt 

Nı+Ng 

2 

Da nun aber durch die Angaben der Anzahl der A—A-Bindungen zwar nicht 
eine bestimmte, so doch aber eine begrenzte Anzahl von Konfigurationen 
(in bezug auf nächste Nachbarn) vorgegeben ist, können die B-Atome nur 
noch jeweils auf die freien Plätze gesetzt werden. Es läßt sich zeigen, daß 
dies für alle Anordnungen immer nur derselben Anzahl von BB- und AB- 
Bindungen führt. Vertauscht man nämlich bei einer vorgegebenen Kon- 
figuration ein A- und ein B-Teilchen, wobei das A-Atom r,A-Nachbarn und 
z—r, B-Nachbarn, und das B-Teilchen r, A-Nachbarn und z—r, B-Nach- 
barn hatte, so ist unmittelbar klar, daß nach der Vertauschung alle A—A- 
Bindungen, A—B-, und alle AB-Bindungen AA- bzw. BB-, sowie alle BB- 
Bindungen AB-Bindungen werden. Eventuell gemeinsame Bindungen der 
vertauschten Atome ändern ihren Charakter nicht; sie brauchen nicht 
berücksichtigt zu werden, weil man das gleiche Ergebnis für die Änderung 
der Bindungszahl erhält, wenn man diese ausläßt und die Koordinationszahl 
entsprechend erniedrigt. 


z Bindungen. 


Die Anzahlen der Bindungen sind: 


vor der Vertauschung |nach der Vertauschung | Differenz 
(II, 44) AA Te L, | Ir. 

AB z— nr, +5 I, T2—-n — 2 (nr, —r,) 

BB Z—T, zZ —Tj | fir 


Man sieht daraus, daß sich die Gesamtzahl der Bindungen bei einer 
Vertauschung immer nur um ein positives oder negatives Vielfaches von 
AA + BB— 2 AB ändert. 


Darin liegt bereits der Beweis der oben aufgestellten Behauptung. Wir 
können jetzt die Gleichung für die Energie ableiten, wenn wir annehmen, 
daß die Anzahl der AB-Bindungen (= zx) bekannt ist, und im Zustand 


getrennter Phasen die Anzahl der AA-Bindungen en und entsprechend — 
für die BB-Bindungen betragen. 2 


Die Energie ist dann durch na Ausdruck gegeben. 
Z 


(OEL, 45) E= > Nu) Vat- 5 N —x) VB + zxV4R 


9 


2 


== = INA VER + NpVpp) +zx (Yan Yasııh a), | 
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Das Integral der Zustandssumme ist, wenn man den kinetischen Anteil 
vernachlässigt, 


Se 


N N 
wobei dv“ dvp" das Volumenelement des 3(NA+”B) dimensionalen Kon- 
figurationsraumes darstellt. Nach dem zweiten Mittelwertsatz der Integral- 
rechnung gilt 


(II, 46) 
ex NV HN + EN | f.../ antane 
A EN aa + NB BB tx | 3 VA NUN, 


[z NaVaa+N8VpB +x(2Vap — Vaa— Vap)]| BE dv? 
A Be> 


2 kT 


Die Ausführung der Integration von Gleichung (II, 46) ist nun möglich. 
Sie ergibt nach Berücksichtigung aller möglichen Permutationen der A- und 
B-Teilchen, die zu neuen Anordnungen führen 
Seal 
Nun 
dabei sind v, und vg die mittleren freien Volumina der Teilchen A und B. 
Aus Z lassen sich nun die interessierenden thermodynamischen Größen der 
regulären Lösungen ableiten. Das Hauptproblem ist hier die Bestimmung 
von x als Funktion der Temperatur und der Wechselwirkungsenergien. 
GUGGENHEIM machte für x den Ansatz 
(IL, 47) x? —= (N, — X) (N —X)exp— 2V/kKT. 
RUSHBROOKE hat nun durch eine Äquivalenzbetrachtung mit der BETHE- 
schen Näherung gezeigt, daß der richtige Ansatz wie folgt lauten muß 
(II, 48) x2—= (N, — x) (N, —x)exp—2V/kT, 
wobei x durch 


en N N 
x f... [sp —xVikT dad: = 

35 N N 

Sfr VvIRTans dvp” 


Wir wollen hier auf den Beweis nicht eingehen, zumal wir bereits im 
nächsten Abschnitt die Bedeutung dieser Gleichung kennen lernen werden. 
Diese Gleichung in Verbindung mit der Gleichung 
0x 
01 
die man unmittelbar aus den Definitionsgleichungen für X und x ableiten 
kann, vorausgesetzt, daß wir, wie im Vorhergehenden immer angenommen 
wurde, V als temperaturunabhängig ansehen dürfen, ermöglicht die Berech- 
nung aller wichtigen thermodyynamischen Eigenschaften derregulären Lösung. 
Die Übertragung dieser Ergebnisse auf Kristalle liegt nach den einleitend 
hervorgehobenen Bemerkungen auf der Hand. GUGGENHEIM selbst (27) 
hat gezeigt, daß die Ergebnisse der Theorie auf alle sogenannten regulären 


definiert ist. 


zs-r-T 
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Ansammlungen anwendbar sind. Zu diesen gehört auch der Mischkristall 
mit Überstrukturbildung oder Entmischung mit den in allen bisher behandel- 
ten Theorien ohnehin gemachten Einschränkungen. Allerdings konnte der 
Einbau der Fernordnung in die Theorie von GUGGENHEIM noch nicht vor- 
genommen werden. FOWLER und GUGGENHEIM (26) haben bewiesen, daß 
diese Einschränkung fallen gelassen werden kann, und damit beginnt diese 
Theorie auch vom kristallographischen Standpunkt sehr interessant zu 
werden. Wir wollen hier wenigstens für den Fall AB, und zwar zunächst für 
das Mischungsverhältnis 1 : 1 die Übertragung der Methode der regulären 
Lösungen auf einen Kristall mit Fernordnung vornehmen. Wir haben in 
diesem Fall also N—A-Atome und N—B-Atome (insgesamt also 2 N-Atome). 
Die Koordinationszahl, die für beide Atomarten wieder gleich angenommen 
wird, sei z. Wir teilen das Gitter wieder in a-Plätze und b-Plätze ein, und 
nehmen an, daß a-Plätze nur von b-Plätzen und umgekehrt b-Plätze nur 
von a-Plätzen umgeben sind. Für die einzelnen Bindungspaare führen wir 
folgende Symbole ein: 


(11, 49) [0,1] = Anzahl der AA-Paare mit A-Atomen auf a-Plätzen mit 
A-Atomen auf b-Plätzen; 

[1,0] = Anzahl der BB-Paare mit B-Atomen auf a-Plätzen mit 
B-Atomen auf b-Plätzen; 

[0,0] = Anzahl der AB-Paare mit A-Atomen auf a-Plätzen mit 
B-Atomen auf b-Plätzen; 

[1,1] = Anzahl der AB-Paare mit B-Atomen auf a-Plätzen mit 
B-Atomen auf a-Plätzen. 


Unter Berücksichtigung der Definitionen (LI, 49) gelten folgende Beziehungen 
(II, 50) [0,0] + [0,1] = [0,0] + [1,0] =zNr 

warez de 
r ist hier die Anzahl der ‚richtig‘‘ besetzten A- und in diesem Fall also 
auch B-Gitterplätze. Setzen wir in Analogie zum GUGGENHEIM’schen Ansatz 
[0,1] = zNx, so lassen sich mittels Gleichung (II, 50) die Anzahl der übrigen 
Paare daraus berechnen. 


Es ist 
(II, 51) 
[00] =z2N(r—x); 1] =zNxs; [1,0] =zNx; [,J =zNIA <r 33 
Damit berechnen wir die Energie 


(IL, 52) E=zN V,g 4 zNx (Var —+ Vpp —2V ip) 


7 T 
—zN (Van +2xV); vor t Ve u, 
Die quasichemische Gleichgewichtsformel lautet 
[0, 0] : BE 1] 
II, 53 [0,0 N navıkm 
a KA SN: 0] exp V/kT 


Man beachte hier die Verwandtschaft der analogen Beziehungen, die bei 
der Behandlung der BRAGG-WILLIAmS-Theorie in Gleichung (II, 9) gefunden 
wurde. Der einzige Unterschied von (II, 9) und (II, 53) liegt in der Definition 
von V. Hier ist die Analogie zum chemischen Massenwirkungsgesetz voll- 
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kommen klar. Man behandelt die A- und B-Atome auf den einzelnen Plätzen 

wie verschiedene Teilchen eines Gases, die Moleküle (Paare) bilden können, 

und wendet auf diese die bekannte Gleichgewichtsformel an. Setzen wir die 

Ergebnisse (IT, 51) in (II, 53) ein, so erhält man die Beziehung 
(r—x)(1—r—x)=x?exp2V/KkT, 

deren Wurzel 


Vırdrl r)fesp2Vkr 1] 1, 
2(exp2 V/kT— 1) a 
Wir führen nun an Stelle von r den BraaG-WirLıams’schen Fernordnungs- 
parameter S ein (Gleichung II, 1) fürF,=F, =#; 
1 
7 n—_ RZ 
Führen wir das Ergebnis in (II, 54) ein, so erhalten wir für x 
V1+4—8% (ep2VkKT—N)—I 
2(exp2 V/kT—|]) 
Damit ergibt sich durch Einsetzen in Gleichung (II, 52) die Energie als 
Funktion der Temperatur 
E Vıp ir 
SW VE Sy eVeop2ViETn-1 
ı Die uns hier interessierende freie Energie erhalten wir durch Integration 
' der thermodynamischen Beziehung (siehe Gleichung (II, 39)). 


F 

m) 
| o (1/T) 
' Wählt man T = & als untere Grenze der Integration, so ergibt sich die 
Integrationskonstante zu —kInW (folgt aus der Zustandssumme) mit 

F N! 
en 

I  LNTn!IN A—n)]! 


‘“ Führt man an Stelle von r wiederum S ein, so erhält man nach Ausführung 
der Integration unter Anwendung der StirLin@’schen Formel für W 


st. 


(II, 54) Zi 


— 27 1\@r 


' (II, 55) Z- 


(II, 56) 


2 F 1 
=. H, also m fE d 2) —- const. 


12 


BI DVa2 te snmazstu—- ni Sn 
NkT 
H [ Ä EA a—S 9] a+1] 
+42 | - es (1 S) Ing n—, | 
mit a= VYl+ (1-8) (exp2V/kT—1). 


; 0 
Für das thermodynamische Gleichgewicht erhalten wir nun aus FR 0 


a+s u 2\] 1+S 
Ba z ze 


(II, 58) In 
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Eine nähere Untersuchung der Gleichung (II, 58) enthüllt grundsätzlich 
die gleichen Beziehungen wie die bisherigen Theorien. Insbesondere erhält 
man für die kritische Temperatur exakt den gleichen Wert, den man bei 
der Durchführung der BETHE’schen 1. Näherung für den gleichen Fall aus 
Gleichung (22) erhält. Dieser Wert war ja dort durch 

x—=1— 2/z (mit x = exp — V/kT) 
V 2 V z 
gegeben, also — ıT In ( — -) oder IT In = 
die gleiche Beziehung erhalten wir auch aus Gleichung (II, 58). 


Aber nicht nur die kritische Temperatur, sondern auch alle übrigen 
Ergebnisse sind mit dem BErHe’schen Ergebnis der 1. Näherung völlig 
äquivalent, der Beweis dieser Tatsache wurde von FOWLER und GUGGEN- 
HEIM gegeben und soll hier nicht durchgeführt werden. Diese Äquivalenz 
liegt nicht nur für das hier geführte Beispiel vor, sondern gilt ganz allgemein. 
Damit gewinnt die quasichemische Methode ihre besondere Bedeutung, sie 
ist in allgemeinen Fällen sehr viel einfacher durchzuführen als die BETHE’sche 
Berechnung. Ihre Hauptschwierigkeit liegt oft in der Integration der Energie- 
gleichung zur Gewinnung der freien Energie. Schon die von FOWLER und 
GUGGENHEIM mitgeteilten Berechnungen eines A—B-Gitters mit variablem 
Mischungsverhältnis von A und B läßt die Vorzüge der quasichemischen 
Methode deutlich werden; denn hier liegen für das BETHE’sche Verfahren 
Berechnungen von EASTHoPE (18) und CHAnG (20) (21) vor, die im Verhältnis 
zu der hier im folgenden kurz skizzierten Darstellung erheblich komplizierter 
sind. Wir wollen hier diesen allgemeineren Fall kurz streifen, um die beson- 
deren Vorzüge der quasichemischen Methode herauszustellen. 

Wir bezeichnen mit N r die Anzahl der A-Atome auf a-Gitterplätzen 
und mit Nq die Anzahl von A-Atomen auf b-Gitterplätzen. (N (r+dq) 
—= ,2N = Anzahl der A-Atome). Die entsprechenden Größen für die B- 
Atome sind dann also N (1—r) und N (1—.q) (Anzahl der B-Atome 
— N [2 — (r + q)] = 9 2N). Die entsprechende Anzahl der Paare A auf 
a- und B auf b-Plätzen bezeichnen wir wieder wie im vorigen Fall mit 
[0,0], [1,0], [0,1] und [1,1], und nehmen an, daß die Anzahl [0,1] =zNr 
sei. Dann gelten für die einzelnen Paarzahlen folgende Beziehungen 
(II, 59) [0,0])=zN(r—x); [(,1]=zNx; [1,01)=zN1—r— gq+x% 

l,1]=zN (q—x). 
Gleichung (II, 53) wird dann mit den Werten (IT, 59) zu 
(IL, 60) r—x)(g—x)=x(1—r— q+xXx)exp2V/kT. 
Die Energie E wird in diesem Fall gleich 
A161) E=zN(r+q)Vss +tzN (1—r— g)Vgg tzNx2V. | 

Wir müssen nun den Fernordnungsparameter S neu festsetzen, für den 

ja die Definition in Gleichung (II, 1) voraussetzt, daß im geordneten Zu- 


stand (S = 1) alle A-Atome sich auf a-Gitterplätzen befinden müssen. Die 
Jetzt eingeführte Definition von S lautet | 


(II, 62) r=y(l+9; q=,Ml—8. 
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Man sieht ein, daß für S=1Nr=2NO,, also alle A-Atome auf 
a-Plätzen sind und fürrS=0Nr=NO,, also die Hälfte aller A-Atome auf 
a-Plätzen sind. Die Beziehungen für I —rund 1 —qsind, da, + 9% =1 

1—r= Op — 9,Sundl — O5 + O8: 
Die mögliche Darstellung mit ©, allein wurde hier aus Symmetriegründen 
) vermieden. Man kann nun mit Hilfe der Lösung für Gleichung (II, 60), die 
/ man in Gleichung (II, 61) einsetzt, unter Berücksichtigung der Beziehungen 


J (IT, 62) die Gleichung für die Energie aufstellen und erhält 


R E—-E, V...—1-— (0.—9,) (exp2 V/kT—1) 
E63) NzV_ exp2 V/kT —1 


V1+[2(05-9,)+4922(1-8)] (exp 2V/kT-1) + (05-04)? (exp 2 V/kT- 1) 


Die freie Energie F gewinnt man dann durch die analoge Integration der 
thermodynamischen Beziehung (II, 39). Da zur Bestimmung des Gleich- 

F ß 
gewichts nur rn benötigt wird, haben FOWLER und GUGGENHEIM die Inte- 
gration nur für diesen Ausdruck vorgenommen. Wir wollen diese neue formal 
mathematische Rechnung nicht durchführen, sondern nur noch das Ergebnis 
angeben: 


(II, 64) 


I OF] A+8)(% +98) 

NK os IE (1—S) (98 — QS) 
| ER NIEREN A Ve 

| Di exp 2 V/kT(1—S) (0, — 0,8) 


nit y... = 


V1+l2(0-9,)+492? (1-8)] (exp 2V/kT-1)+ (9-04)? (esp2V/kT-1)? 


Nullsetzen dieses Ausdrucks ergibt die Gleichgewichtsbedingung. Für die 
kritische Temperatur erhält man die Bedingungsgleichung 


a2 
z Es O4 05 


(II, 65) exp2V kl. > 


DIE men 

Dieser Ausdruck wurde bereits von EASTHOPE (18) unter Verwendung 
der BerHe’schen 1. Näherung gefunden. Die graphische Darstellung der 
Gleichung (IT, 65) ist bereits in Fig. 5 wiedergegeben worden. Weitere 
Untersuchungen über den AB-Fall mit variabler Konzentration findet der 
Leser in der schon erwähnten Arbeit von CHAn&G (21). 


In diesem Zusammenhang ist es lohnenswert, eine weitere Methode kurz 
zu streifen, die eigentlich die logisch einfachste, aber mathematisch im 
allgemeinen recht kompliziert ist. Es handelt sich dabei um eine rein kom- 
binatorische Methode, die direkt von der Zustandssumme ausgeht. Wir haben 
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den Gedankengang schon bei der Besprechung des eindimensional ferro- 
magnetischen Modells von Isına (13) besprochen. Es liegt nahe, diese 
Methode auch auf drei Dimensionen zu übertragen. Wenn man in der Zu- 
standssumme 
Z= 2 exp — E,/kT 
Konf. 

alle Glieder mit gleichem Wert für E; durch einen mathematischen Ausdruck 
angeben kann, ist die Durchführung der Summation zur Bildung der Zu- 
standssumme möglich. Die Energie einer Konfiguration läßt sich nun mit 


a ). 1.0 
os), 


0. 
02 


0 
04 02 0,3 0.4 05 ®R 


Fig. 5. Kritische Temperatur als Funktion der Konzentration für 
z = 6 (NaCl-Typ) 
z=8 (OsCl-Typ) 


Hilfe der x gemäß Formel (II, 61) angeben. Die Anzahl W; der Kon- 
figurationsmöglichkeiten mit der durch E; gegebenen Energie, die für den 
eindimensionalen Fall streng berechnet werden kann, wurde von CHANG (28) 
(29) durch direkten Vergleich mit der BETHE’schen Methode bestimmt. Das 
rein formal mathematische Verfahren sei hier nicht erläutert. Diese Methode 
wird dadurch näher begründet, daß CHang für das eindimensionale Modell 
die Identität seines Ergebnisses mit der mit Hilfe der BETHRr’schen ersten 
Näherung erzielten Formel für z = 2 nachweist. FOWLER und GUGGENHEIM 
(26) können auf der anderen Seite zeigen, daß die Gleichgewichtsbedingung 
der Paare in der quasichemischen Methode unmittelbar erhalten wird, wenn 
man eine kombinatorische Näherungslösung für die Anzahl der gesuchten 
Verteilungsmöglichkeiten der Paare für ein bestimmtes x ausführt und dabei 
nicht berücksichtigt, daß im Kristallgitter die Paare eigentlich nicht un- 
abhängig voneinander betrachtet werden dürfen (non-interference of pairs). 
Die dadurch erhaltene kombinatorische Formel muß natürlich wegen der 
schon gezeigten Äquivalenz der quasichemischen Methode mit der BETHE- 
schen ersten Näherung mit der von CHANG (28) erhaltenen Formel identisch 
sein. 

Wir können im Rahmen dieser Zusammenfassung nicht auf die in dieser 
Richtung erschienenen Arbeiten eingehen, zumal die darin enthaltenen 
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Gedankengänge rein formal mathematischer Art sind und keine neuen 
physikalischen Gedankengänge entwickeln. Ich möchte jedoch auf einige 
wichtige Erweiterungen, die von CHAnG und Ho (30) durchgeführt wurden, 


) hinweisen. 


Wichtig erscheint mir im Zusammenhang mit der quasichemischen 
Methode noch eine neuerdings von Yang (31) angegebene Verallgemeine- 


! rungsmöglichkeit zu sein, die den Vorteil besitzt, daß sie erstens bessere 


Näherungen erzielt und weiterhin auf kompliziertere Systeme leicht an- 
gewendet werden kann. Der Vorteil des Verfahrens liegt darin, daß die für 
den A—B-Fall schon recht schwierige Integration der Energiegleichung 
hier durch eine LEGENDRE-Transformation ausgeführt wird. Der physikalische 
Gedankengang ist aber der gleiche wie bei FOwLER und GUGGENHEIM und 
wir können uns deshalb in der Darstellung auf äußerste Kürze beschränken. 
Wir erläutern das Verfahren am A-B-Fall mit konstantem Mischungs- 
verhältnis 1 :1. Wir ordnen den [0,0]-Paaren die Atomart x, den [1,0]- 
Paaren die Molekülart xy, den [0,1]-Paaren die Molekülart xz und den 


- [1,1]-Paaren die Molekülart xyz des dreiatomigen Gases mit den Atomarten 


x, y,z zu. Man sieht sofort, daß die Anzahl der x-Atome gleich 


(HI, 66a) [0,0] + [1,0] + [0,1] + [L1] =Nz, 


die Anzahl der y-Atome gleich 
(II, 66b) 0,1 + [LI =Nz (Ir), 
und die Anzahl der z-Atome gleich 


(LI, 66) [1,0)+[1,1]=Nz(1—.r) ist. 


Die Definitionen von r und z sind die gleichen wie in (IL, 50). Wir wollen 


| dabei in Erinnerung behalten, daß die Paare [0,0] und [1,1] A-B-Bin- 


dungen, [1,0] B-B-Bindungen und [0,1] A-A-Bindungen darstellen. Es 
sollen jetzt für die mittlere Anzahl der einzelnen Paaranzahlen die Gleich- 
gewichtsformeln des dreiatomigen Gases aufgestellt werden, die bekannten 
Gleichungen sind (siehe FOWLER, Statistical Mechanics (32): 

[0,0]=£ exp — Vır/KT, 
(LI, 67) [0,1] =,&% exp — Vaa/KT, 

[1,0] = £u exp — V»s/KT, 

[1,1] = Euv exp — Vır/KT, 


wobei die Größen £, », u aus den Bedingungsgleichungen (II, 66a—c) zu 
ermitteln sind. Wir erhalten daraus für die mittlere Energie 


(II, 68) E= (£exp — Vap/kr) Van + (Evexp — Varfır) Vaa 
+ (Eu exp — Vrp/kr Vpp + (Evu exp — Vap/er) Var: 
Führen wir die Hilfsfunktion 
9 (£&,»,u, T) = Eexp— Vap/ır + Evexp— Var/ır + Euexp — Vap/kr 
+ Eru exp — Var/kr 
ein, so sieht man unmittelbar, daß die Energie E sich durch p ausdrücken 
läßt. Es ist offenbar 


(II, 69) E=krT 


128 H. Jagodzinski, 


Es gelten folgende Beziehungen 


(II, 70) et <aN; nzaNd—n); ee 
Hier haben wir so getan, als ob die £, «u, v und T unabhängige Variable 
seien; in Wirklichkeit sind sie es aber nicht, sondern Funktionen von T undr. 
Faßt man daher die Gleichungen (II, 70) als Definitionsgleichungen auf, 
so kann die LEGENDRE-Transformation 
In&e>zN, n»>zN(l—r, nı>zN(1—r), T>T 


einführen. 


Die Integration der Energiegleichung gestaltet sich nun sehr einfach. 
Die thermodynamische Beziehung 
0 (F/KT) E °p 
ol  Km—aT 
zeigt uns unmittelbar den Zusammenhang zwischen der freien Energie und 
unserer Funktion 9. Die noch offene Frage der Integrationskonstanten wird 
dabei durch die LEGENDRE-Transformation gelöst. Für den Leser, dem die 
Technik der LEGENDRE-Transformation nicht geläufig ist, sei hier die Inte- 
gration über das totale Differential durchgeführt. Wir fassen 9 (£, »v, u, T) 
als Funktion der unabhängigen Variablen auf und bilden das totale Diffe- 
rential von 


(II, 71) 


0Y 
9u 


0R, 


IV 


de 


| 
T 


OP am 


(IL, 72) do IT ; 


de- dv +4 


ä ) i j 
Gesucht ist der Wert des ” m T, integriert von T = » als untere 
0 


Grenze und T als obere Grenze. Das totale Differential hat nun die angenehme 
Eigenschaft, daß die Integration unabhängig von Integrationsweg ist. Wir 
können also hier ruhig die Variablen als unabhängig betrachten, und müssen 
nur beachten, daß die für die einzelnen Variablen eingesetzten Grenzen 
mit den durch T vorgeschriebenen Grenzen im Einklang sind. Wir setzen 
die Beziehungen (II, 70) in (IL, 72) ein und erhalten 


d 


ze 
£ 


Ss 


9% am 


N 
a OT 


Ne Sr 
v [2 


Integrieren wir diesen Ausdruck von T = & bis T, so ergibt sich 


= ar 
2) 
(11473) [RAT p—aNin es —aN nv —zN (nina 

YroEl 

[6>) T=o 
Ah 

= Mlrsn) = % , D-Xf, 5) 

N= cs 


mit 


(1,74) Fe, D=g,W—-zNm:—-zN(I-ninv»—zN(l—r) In u. 
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Das ist aber das gleiche Ergebnis, das uns dieDurchführungder LEGENDRE- 
‚ Transformation liefern würde. Die Werte für £&, », u, lassen sich aus den 
| Gleichungen (II, 67) und (II, 66a—c) ermitteln. Für T = ist 


R Ir 
SEND we 


r 
Für T führen wir die Rechnung nicht aus. 


Damit erhalten wir 


AN 
F ö 
(EI, 75) im) Fam —eNhe zN Aria 


—zN (ln) Invep — 95 +tzNinzNr2+2zN(l—rIn en 
Nr EB? 
Sn 2 BE 
Ai TT un vn 


da oT = po =ZN ist. 


Weil wir hier an sich nur die Methode des Verfahrens zeigen wollten, 
können wir an dieser Stelle mit der Weiterführung der Rechnung ab- 
schließen. Man sieht unmittelbar, daß der quasichemische Ansatz nicht nur 
auf Bindungspaare, sondern auch auf Atomgruppen angewendet werden 
kann, ohne daß die Berechnungen irgendwie gestört werden. Dabei würde 
lediglich die Anzahl der Bedingungen (II, 66) und (II, 67) erhöht, das 
Rechenverfahren bleibt grundsätzlich das gleiche. YanG hat diese Methode 
auch auf den Fall AuCu, ausgedehnt, dabei allerdings nicht die gesamten 
Gruppen der zwölf Nachbarn, sondern nur eine tetraedrische Gruppe von 
ı Atomen erfaßt. Eine weitere Anwendung auf das System Au—Cu findet 
der Leser in einer neueren Arbeit von Lı (33). 

Damit sollen die Erörterungen über die quasichemische Methode ab- 
geschlossen werden. Die Vorzüge liegen offenbar in ihrer mathematischen 
Einfachheit und deswegen großen Verallgemeinerungsfähigkeit. Insbeson- 
dere ist das von Yana angewandte Verfahren der Einführungen der 
LEGENDRE-Transformationen sehr leicht durchführbar, da sich die Grund- 
gleichungen im Prinzip nicht ändern. Die ursprüngliche Beschränkung der 
Genauigkeit auf diejenigen der BurHe’schen 1. Näherung wurde durch das 
_ letztere Verfahren ebenfalls beseitigt. Damit scheint die quasichemische 
' Methode eines der wertvollsten Hilfsmittel zur theoretischen Bearbeitung 
komplizierter fehlgeordneter Mischkristalle zu sein. Die Beschränkungen 
ihrer Anwendbarkeit sind natürlich die gleichen wie bei den bisher be- 
handelten Problemen, da auch in dieser Theorie eventuelle Gitterstörungen 
mechanischer Art nicht erfaßt werden können. 


II, 6. Die Methode der Cluster-Integrale. 


Wir wollen jetzt eine weitere Methode der Behandlung unseres Fehl- 
ordnungsproblems diskutieren, die wesentlich allgemeiner als die voran- 
gegangenen Thecrien, aber in der numerischen Ausrechnung ziemlich müh- 
sam ist, und deshalb noch nicht sehr oft auf Kristalle angewendet wurde. 


Fortschritte der Mineralogie 1949. Bd. 28. 9 
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Die Methode der Cluster-Integrale geht zurück auf einige Arbeiten von 
MAYER und Mitarbeitern (34) (35) (36) (37), in der die statistische Behand- 
lung des Kondensationsproblems realer Gase abgeleitet wird. Eine gute 
zusammenfassende Darstellung über dieses Gebiet findet der Leser in einer 
Arbeit von Born und Fuchs (38) sowie in dem Buch ‚Statistical Mechanics“ 
von MAYER und GÖPPERT-MAYER (39). Von Fuchs (40) wurde diese Methode 
auf Kristallgitter übertragen; diese Arbeit gewinnt deshalb besonderes 
Interesse, weil in ihr gerade die andere Seite des Fehlordnungsproblems, 
nämlich die Entmischung, besonders berücksichtigt wird. 

. Zum besseren Verständnis wollen wir hier zunächst auf die Behandlung 
der Gase eingehen, um die Ergebnisse dann in analoger Weise auf Kristall- 
gitter zu übertragen. Ich will im folgenden weitgehend die Bezeichnungs- 
weise von BoRN und Fuchs (38) verwenden, da ich hier die Berechnung 
nicht explizit durchführen werde; der Leser kann dann in einfacher Weise 
die entsprechenden Beweise in der Originalarbeit nachlesen. 


Die Zustandssumme eines Gases wird allgemein dargestellt durch 
1 » » 
A N | al exp (TH U)/kTdp; ... dpardd, dam 


Hierbei ist 
T = kinetische Energie 
U = potentielle Energie 
h = Pranck’sches Wirkungsquantum. 


Die Integration erfolgt über den gesamten Phasenraum der 3N-Orts- und 
3N-Impulskoordinaten. Wenn wir annehmen, daß die kinetische Energie 
unabhängig von der jeweils eingenommenen Konfiguration ist, kann die 
Integration über den Impulsraum getrennt vollzogen werden; das ergibt in 
bekannter Weise 


2 nm kT\ 921 r 
(IT, 76) 2 5 sn f.. fer VeTda --- dam. 


Liegen nun Zentralkräfte vor, so ist die potentielle Energie nur eine 
Funktion des Abstandes und nicht der relativen Lage der Teilchen zu- 
einander. Dann kann man das Potential irgendeiner Konfiguration aus- 
rechnen. Es ist 


(IL, 77) IN 
ISIS JEN 


r;; ist dabei der Abstand des i-ten Teilchens vom j-ten Teilchen. Die Sum- 
mation hat unabhängig über alle Teilchen zu erfolgen, wobei die Ein- | 


schränkung 1 <i< j < N bedeuten soll, daß jedes Potential nur einmal 


gezählt wird. Jedem Punkt im Konfigurationsraum ist also ein bestimmter 
Wert von V zugeordnet. An Stelle der v (r;;) führen wir die Funktion f;; 
ein, die durch die Beziehung 


(II, 78) exp —v(n)/kT)—1=1f; 
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definiert ist. Dann wird das Integral (II, 76) über den Konfigurationsraum 


-/. ./\+22%+ 5 fuhr + day... dam. 


1>3 i>7 - 


Wir führen jetzt dieselbe Vernachlässigung ein, die bisher allen Metho- 
den zugrunde lag, nämlich v (r;;) = 0 für Teilchen, die nicht unmittelbare 


‘ Nachbarn sind. Das vereinfacht die Auswertung der Summen in Gleichung 

(IL, 79) sehr, da wir bei der Integration nur solche Teile des Konfigurations- 
‚ raumes zu berücksichtigen haben, für welche die in den Summen stehenden 
Produkte nicht verschwinden. Es kann also die Integration der verschiedenen 
- Glieder der Summen nur über ein Teilchen unabhängig durchgeführt werden, 


die anderen Teilchen geben ja nur einen Beitrag, wenn sie sich in der un- 
mittelbaren Nachbarschaft dieser Teilchen befinden. Damit ergibt sich ein 
Zuordnungsprinzip zu jedem Produkt der f;; zu Cluster, dessen Bedeutung 
wir an einem Beispiel erläutern wollen. 


Das Glied 


IEPE Pe PETE TE eT Ber 
in dem die einzelnen Faktoren nur dann Anteile ergeben, wenn die durch 
die Indizes gekennzeichneten Teilchen benachbart sind, zeigt, daß die 


' Teilchen 1—4, 5—7 und 8,9 Cluster bilden müssen, wenn ihre Beiträge zur 
‚ Integration über den Konfigurationsraum nicht verschwinden sollen. Wir 


haben also im obigen Beispiel eine Einteilung in 1 Cluster von 4, 1 Cluster 
von 3 und 1 Cluster von 2 Molekülen. Die Indizes geben an, welche Ver- 
knüpfungen bestehen müssen, ohne jedoch zunächst auf die geometrische 


_ Realisierbarkeit Bezug zu nehmen. Das ist eine Aufgabe, die dann bei der 


vorzunehmenden Integration über den Konfigurationsraum berücksichtigt 
werden muß. Alle übrigen Teilchen müssen in bezug auf das obige Glied als 
frei beweglich angesehen werden, da ja deren Nachbarschaftsanteile nicht 
in diesem Glied enthalten sind. Man sieht auch unmittelbar, daß die Cluster 


unter sich frei beweglich sind. Wir fassen nun alle möglichen Cluster von 
insgesamt 1 Teilchen zu einem Clusterintegral bı zusammen, das wir wie 
_ folgt definieren: 


(II, 80) 1b /... [zn f4 dgı --- dgsaı) 


Die Summierung im Integral muß über alle Möglichkeiten, auf die man 
einen Cluster von vorgegebenen 1 Indizes bilden kann, erfolgen. Die Inte- 
gration darf nur über I—1 Teilchen ausgeführt werden, da, wie vorhin 
bemerkt wurde, ein Teilchen als frei beweglich angesehen werden kann und 
nur den Beitrag V im Konfigurationsintegral liefert. 

Jeder Term unseres Integrals über den Konfigurationsraum enthält 
m, Cluster, bestehend aus 1 Teilchen, m, Cluster aus 2 Teilchen, ...... mj 
Cluster aus I Teilchen. 


132 H. Jagodzinski, 


Dabei muß natürlich die Nebenbedingung 
al mj = N 
1 


erfüllt sein. Das statistische mechanische Gewicht dieses Terms ist durch 
die Anzahl der Permutationen der N-Teilchen bestimmt, allerdings muß 
man dabei beachten, daß die I! Permutationen des Clusters 1 und die mı! 
Vertauschungen der gleichen Cluster keine neuen Anordnungen ergeben. 
Damit erhalten wir für das statistisch mechanische Gewicht des betrachteten 
Terms 


N! 1 
— N! ee 
(1m (21m... (Ihm my!... mi! 1 (Um mı! 


Das gesamte Integral über den Konfigurationsraum erhält man durch die 
Summierung über alle möglichen Kombinationen von m, Oluster der 
Grobe laser mı Cluster der Größe l mit Nebenbedingung Zlmı =N 


= (VI! Bm (Nobjm V 
(II, 81) f- ‚fe! u on -n!y Il, mito=g 
mı 1 


=oN] 


v = mitt]. freies Volumen pro Teilchen. 


Damit ist das Problem der Ermittlung der Zustandssumme auf die 
Auswertung der sogenannten Cluster-Integrale zurückgeführt worden. Für 
die Auswertung wird von MAYER und GÖPPERT-MAYER eine Zerlegung der 
Cluster-Integrale in reduzible und irreduzible Cluster vorgenommen, die 
eine große Vereinfachung bedeuten. Unter einem irreduziblem Cluster ver- 
steht man einen Cluster, dessen zugeordnete Faktoren f;; nicht in zwei Grup- 
pen aufgespalten werden können, daß beide Gruppen nur einen gemeinsamen 
Index besitzen. Ein Beispiel möge das erläutern. Der zu fj3 fa3 fs fy4 fs 
gehörige Cluster ist reduzibel, er kann nämlich in die Gruppen fj3 fa; fı3, 
fz4, 15, aufgespalten werden, so daß jede Gruppe nur einen Index gemeinsam 
hat. Der Cluster fj, fa; fj, ist dagegen irreduzibel, denn jede weitere Ab- 
spaltung bringt immer eine Gruppe mit zwei gemeinsamen Indizes zu- 
stande. 


Man kann nun auf rein kombinatorischem Wege zeigen, daß jedes 
allgemeine Cluster-Integral in eine Summe von Produkten irreduzibler 
Cluster-Integrale zerlegt werden kann, auf diese Weise ist das kombinatori- 
sche Problem der Auffindung aller möglichen Cluster der Größe 1 wesentlich 
vereinfacht. Wir wollen hier den Beweis nicht antreten, der Leser findet ihn 
ausführlich dargestellt in der bereits erwähnten Arbeit von BoRN und 
Fuchs (38). 

Man kann die irreduziblen Cluster auch graphisch darstellen, nur muß 


man dabei beachten, daß diese Darstellung nicht etwa die geometrische 
Konfiguration wiedergibt, sondern nur die im betreffenden Cluster a | | 
vorherrschende ‚Bindung‘ darstellt. Die Frage der Geometrie muß dann erst | 
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bei Ausführung der Integration über den Konfigurationsraum berücksichtigt 
werden. Wir wollen hier die einfachsten irreduziblen Cluster darstellen: 


ßı A irreduzibler Cluster re dis 


Bl ” 5 2 N 
i oe _——® fja fgs fs 


o—® 
Ba 3 > ” zo | | fj2 fy3 tz, fja 
oo ——-® 
e—-® 
6 ” ” Pz31ı DI f]2 fj3 fja tag tz, 
ee —® 
1 ER) » a2 U fj2 fj3 f4 ty; fya tz. 


Die Anzahl der irreduziblen Cluster erhält man durch die Ermittlung 
der möglichen Änderungen der Indizes der dazugehörigen Produkte, ohne 
daß die dazugehörige graphische Darstellung sich ändert. 


Man erhält die allgemeinen Cluster-Integrale aus den möglichen ir- 
reduziblen Clusterintegralen vermittels der Gleichung 


(II, 82) > It; Im: i v bezeichnet dabei die Teilchenzahl des 
irreduziblen Clusters und u, die Anzahl 


Di gr en der irreduziblen Cluster der Größe v. 
> v 


Daß hier Zvu, =1— 1 ist, liegt an der Verknüpfung der irreduziblen 
Cluster bei der Integration. Die Summierung über die u, hat also wieder 
über alle Möglichkeiten der Zusammenstellung der irreduziblen Cluster zu 
erfolgen unter Beachtung der Nebenbedingungen Zyu, —=1— 1. Wir be- 
merken hier noch, daß die Gleichungen (11,81) und (II, 82) sehr verwandt 
sind, und wir aus diesem Grunde die Auswertung dieser Gleichungen, die 
man recht gut mit funktionentheoretischen Hilfsmitteln vornehmen kann, 
mit einem einzigen Ansatz lösen können. Wir wollen die rein formal mathe- 
matische Ausrechnung des kombinatorischen Anteils des Konfigurations- 
integrals und der Cluster-Integrale nur ganz kurz in ihrem Gedankengang 
ausführen. Gemäß (II, 81) und (II, 82) ist 


BEE (I x ymı Mean 
(II, 83) => a = E (ML 2,9% „oe obe 
ee 
Zlm=M 
B, M=1—]1 
Sl > m —P.(M 1,3220) mir 
Tr a 


a — il 
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Wir führen nun eine neue Variable £&, ein und können zeigen, daß die 
Gleichungen (II, 83) mit den Koeffizienten der Entwicklung der Exponential- 
funktion exp (1 2x,£”) identisch sind. Diese lautet 


1 
eXPpl 22, 2 (De 
2 an 
Man kann durch einfache Kunstgriffe der Umnumerierung der Summierung 


und Zerlegung von m zeigen, daß die obige Gleichung durch folgenden 
Ausdruck dargestellt werden kann (siehe dazu BORN und Fuchs (38). 


1] (1x,)* 
= N eM 
xplrne)—= NE > Ei 


Zvu,=M 

Definiert man nun eine neue Funktion 
rules = 22, pi, 

welche die Eigenschaft besitzt: . 

(II, 84) m G=z En 

dann ist 


EXP L2x, EPG Er 2 (ME 81,8%...) 


Nach dem CAaucHv’schen Integraltheorem ist dann 
1 Wr) 


2ni EM+1 


(IL, 85) ENT ze 


or 
Ye 


Die in (II, 85) dargestellte analytische Funktion ist offenbar innerhalb 
des Kreises mit dem Radius des Konvergenzkreises von G, (&,Xy-..) 
regulär mit Ausnahme an der Stelle & = 0. Wählen wir zur Ausführung der 
Integration einen Kreis mit dem noch frei zu wählenden Radius € =Zexpi 
bzw. pexpig und führen wir die Lösungen des asymptotischen Wertes für 
große M und I nach der Sattelpunktmethode durch, so erhält man schließlich 
die Näherungslösungen für große N und | 


SEM nt 2 
2 exp1G, (P; = a). 
’bı ( 
al EL E| (P, 2. .) 
mit Gulz eb DD z’ A; Gy (e, Br , B 
il al 


wobei z und go durch die Gleichungen ° 
G sobre 
| definiert sind. 
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Wir wollen an dieser Stelle mit den Betrachtungen über reale Gase 
abbrechen. Die Berechnung der interessierenden thermodynamischen Größen 
ist mit der Rückführung der Integration über den Konfigurationsraum auf 
die irreduziblen Cluster-Integrale als gelöst zu betrachten. Man kann schon 
ohne Kenntnis der einzelnen Integrale sehr wichtige Aussagen machen. 
Allerdings sind sie in bezug auf die Kondensation wohl deswegen nicht gültig, 
weil die irreduziblen Cluster-Integrale als volumenunabhängig betrachtet 
wurden. Das gilt aber sicher nicht im Bereich der Kondensation. Über eine 
Erweiterung in dieser Richtung siehe die bereits erwähnte Arbeit von MAYER 
und HARRISON (35). 


Die Berechnung der irreduziblen Cluster-Integrale ist aber eine sehr 
mühsame Arbeit; sie wurde von MAYER und Mitarbeitern mit dem LOnDonN’- 


; schen Potentialansatz für van-DER-Waars-Kräfte nur für 8, durchgeführt. 


N 


Mit Exponentialansätzen, für welche die Integration einfacher auszuführen 
ist, gelang ihnen die Berechnung bis einschließlich f,. Damit scheint die 
quantitative Anwendung der Theorie sehr in Frage gestellt. Trotzdem sind 
ihre Ergebnisse, die man auch ohne Kenntnis der irreduziblen Cluster- 
Integrale erhalten kann, von weittragender Bedeutung, und sind zweifellos 
ein wesentlicher Fortschritt in der Entwicklung der Theorie realer Gase. 


Es erhebt sich nun die Frage, wie weit die hier angegebenen Berech- 
nungen auf Kristalle anwendbar sind. Wir haben es bisher ja nur mit einer 
Teilchenart zu tun gehabt. Es läßt sich aber leicht zeigen, daß diese Ergeb- 
nisse auch auf einen Mischkristall mit zwei Komponenten übertragbar sind. 
Unser Mischkristall habe insgesamt N Atome, von denen ein Bruchteil 
C,A- und ©, =1-—-C, B-Atome sind. Die Anzahl der A-Atome beträgt 
also C, N und die der B-Atome C, N = (1—C,)N. Wie wir bereits bei 
der Besprechung der quasichemischen Methode gezeigt haben, ist aber die 
Berechnung der Energie durch die Angabe der Verteilung der A-Atome 
gegeben. Bezeichnen wir die Anzahl der A-A-Bindungen einer vorgegebenen 
Verteilung mit x, so ist gemäß Gleichung (II, 44) die Energie für eine 


‘ bestimmte Konzentration durch 
«L, 87) IN x (Via + V5R —2 EN) —- const 


— xV + const 
gegeben. 


Man braucht also nur die Verteilung der A-Atome zur Berechnung 
der Energie. Damit haben wir an Stelle der Cluster im Gasraum nur die 
Cluster der A-Atome auf den N Gitterplätzen zu berücksichtigen. Die Inte- 
gration im Konfigurationsraum wird durch die Summation jedes einzelnen 
Atoms über N Gitterplätzen ersetzt. In Analogie zu Gleichung (II, 78) 
führen wir wieder die f;;, definiert durch 


(IL, 88) exp — V;/kT—1=f;; 

ein, wobei 

vj; = V für Atome, die nächste Nachbarn sind, also fj; = exp — V/kT—1 
vj; = © für identische Atome, also = —1 


vj; =0 für alle übrigen Atome, also =. 
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Damit wird die Zustandssumme 


N N 
1 

(II, 89) Zexp— const/kT = a >> u 1 ME era 
ns erlegen 


Die Äquivalenz von Gleichung (II, 89) und (II, 79) ist unmittelbar 
evident, und der nun folgende Gedankengang ist vollkommen gleichwertig 
mit den Überlegungen, die wir bei der Behandlung der realen Gase angestellt 
haben. Statt der Integrationen haben wir es hier nur mit Summationen zu 
tun. Die Cluster-Integrale werden also zu Cluster-Summen, die irreduziblen 
Cluster-Integrale zu irreduziblen Cluster-Summen. Auch die funktionen- 
theoretische Auswertung der kombinatorischen Formeln (II, 81) und (II, 82) 
ändert sich nicht, da ja die Cluster-Summen als unabhängige Größen in 
Erscheinung treten und nicht vom Auswertungsverfahren betroffen werden. 
An Stelle von N haben wir natürlich c,N einzusetzen, und statt des Volumens 


N ll 
V tritt hier nur N auf. An Stelle von “ haben wir hier also — — einzu- 
N CN Ca 


setzen. Formal gesehen benötigen wir hier nur die Existenz von N Gitter- 
plätzen, während an die Geometrie des Gitters keinerlei Forderungen zu 
stellen sind. Praktisch sind wir hier aber auf solche Fälle beschränkt, bei 
denen auch die Gittergeometrie unabhängig vom Ordnungsgrad ist, denn in 
diesem Ansatz sind nicht die Verzerrungsarbeiten des Gitters enthalten, 
die auftreten müssen, wenn die Gittergeometrie der A-Atome mit denjenigen 
der B-Atome nicht einigermaßen übereinstimmt. 


Wir wollen die einzelnen von Fuchs (40) gegebenen Anwendungen 
des Verfahrens nicht ableiten. Im Kristall können die Cluster-Summen 
natürlich ausgewertet werden, auch hier ist das Verfahren trotz des Fortfalls 
der schwierigen Integration mit komplizierten Kraftgesetzen noch sehr 
umständlich. Es ist klar, daß diese Methode besonders gut geeignet ist, den 
bisher weniger betrachteten Fall der Unordnungsübergänge zu beschreiben, 
nämlich 


vayat Ye -2Vis 


Bisher haben wir ja immer stillschweigend vorausgesetzt, daß V > 0, also 
bei tiefen Temperaturen Überstrukturbildung vorliegt. V < 0 führt aber 
bei tiefen Temperaturen zu einem Zweiphasengebiet, also zur Entmischung. 
Fuchs leitet nun alle wichtigen Beziehungen, wie die Löslichkeitsgrenze in 
Abhängigkeit von Temperatur und Konzentration, die Festlegung der kri- 
tischen Temperatur für die Entmischung, das Verhalten der spezifischen 
Wärme als Funktion der Temperatur usw. ab. Er vergleicht seine Ergebnisse 
mit den von KIRKWOOD nach der Semi-Invariantenmethode (41) und von 
RUSHBROORE (25) nach der quasichemischen Methode erhaltenen Ergebnissen 
und findet, daß die Reihenentwicklungen bis zu den dritten Potenzen über- 
einstimmen. Das von FucHs aufgestellte Zustandsdiagramm einer binären 
festen Lösung des kubisch raumzentrierten Gitters hat die bekannte Form 
und ist in Fig. 7 in bezug auf ihren qualitativen Umlauf wiedergegeben. 
Alle Ergebnisse beziehen sich aber nur auf das thermodynamische Gleich- 


| 
2 
\ 
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gewicht. Oberhalb der Löslichkeitslinie haben wir den homogenen Misch- 
kristall, während unterhalb Entmischung in zwei Phasen der durch das 
Zustandsdiagramm angegebenen Zusammensetzung eintritt. 


os) 


A 0.2 0.4 0.6 0.8 B 


Fig. 7. Zustandsdiagramm einer binären festen Lösung mit Entmischung. 


II, 7. Ein kinetischer Ansatz für das Entmischungsproblem. 


Wir haben im letzten Absatz das Gebiet der Entmischung kurz ge- 
streift. Das Problem in den bisher behandelten Theorien gewinnt aber 
insofern noch kein besonderes Interesse, weil sie nur in der Lage sind, solche 
Ausscheidungserscheinungen zu beschreiben, bei denen keine Änderungen 
der Gitterkonstanten eintreten. Diese sind aber sowohl vom technischen als 
auch theoretischen Standpunkt viel weniger interessant, als Fälle, bei denen 
die Gittergeometrie des Ausgeschiedenen wesentlich verschieden von der 
des Mischkristalls ist. Durch die bisher behandelten Probleme ohne Änderung 
der Gittergeometrie wird der Anschein erweckt, als ob für V > O beigenügend 
tiefer Temperatur immer Überstrukturbildung und für V < 0 immer Ent- 
mischung eintreten muß; das braucht aber keineswegs der Fall zu sein. 
Starke Veränderungen der Gittergeometrie können sich auch im Fall V > 0 
als „Entmischung‘ zu erkennen geben. Man kann das leicht qualitativ 
einsehen. Nehmen wir an, A- und B-Atome haben die Neigung, unter Aus- 
bildung eines wesentlich geringeren Abstands als er im Wirtgitter realisierbar 
ist, feste Bindungen einzugehen. Die Einnahme eines geringen Abstandes 
zur Gewinnung des Energieminimums von Va bedingt aber eine Ver- 
zerrungsarbeit des Wirtgitters. Für ein Paar allein ist die Verzerrungsarbeit 
zu groß gegenüber dem Energiegewinn aus der Wechselwirkungsenergie V4p; 
dagegen werden unter bestimmten Bedingungen mehrere Paare zusammen 
diese Verzerrungsarbeit leisten können, d. h. es gibt keine „‚ Überstruktur- 
bildung‘, sondern praktisch Ausscheidung von ‚„Überstruktur‘“-Gebieten, 
die besonders dann deutlich in Erscheinung treten, wenn die Konzentration 
der einen Komponente gering ist; das ist aber in solchen Fällen meist der 
Fall, weil die Mischbarkeit durch die Existenz der starken AB-Bindung 
sehr beeinträchtigt wird. Diese Fragen sind für Aushärtung technisch wich- 
tiger Legierungen von großer Bedeutung und ihre Lösungen können sehr zur 
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Klärung mancher interessanten Probleme beitragen. Es ist unmittelbar 
klar, daß für einen solchen Fall nicht Berechnungen mit der Annahme des 
thermodynamischen Gleichgewichts herangezogen werden können, weil 
dieses oft gar nicht erreicht wird. Es ist von viel weitgehenderem Interesse, 
einiges über die Kinetik solcher Vorgänge zu erfahren. Schon allein Platz- 
wechselschwierigkeiten können verhindern, daß sich das thermodynamische 
Gleichgewicht einstellt. 


In dieser Richtung hat R. BECKER (42) (43) einige Ansätze ausgeführt *). 
Sie gehen auf die EINSTEIN-SMOLUCHOWSKT’sche Theorie der Kondensation 
zurück. Wie BECKER und DörınG (44) gezeigt haben, kann man für das 
Kondensationsproblem folgenden Lösungsansatz ausführen. Nennen wir pr 
den Dampfdruck eines Tröpfchens vom Radius r und p. den Dampfdruck 
eines unendlich großen Tropfens, so gewinnt man aus dem Vergleich der 
beim Kondensationsprozeß gewonnenen Arbeit mit der gegen die Öber- 
fläche zu leistenden Arbeit bei der Entfernung der gleichen Maße aus dem 
Tröpfchen folgende thermodynamische Beziehung zwischen dem Dampf- 
druck des Tröpfchens vom Radius r und der Oberflächenspannung. 


Pr _ M 20 M = Molekulargewicht 


Sl pBoere e — Dichte der Flüssigkeit 
daraus folgt co — OÖberflächenspannung. 
2Mo 1 
I —= Doo | 
(II, 90) Pr = pesexp | = =) 


Da p der Dampfdruck des übersättigten Dampfes ist, erkennt man ohne 
weiteres, daß unser Tröpfchen solange nicht wachstumsfähig ist, wie sein 
Dampfdruck über dem des übersättigten Dampfes liegt; d. h. ist r < ry, 
wobei r, durch den Dampfdruck p, definiert ist, bei dem beide Dampf- 
drucke gleich sind, so ist die Wahrscheinlichkeit, daß das Tröpfchen wieder 
aufgelöst wird, größer als die Wahrscheinlichkeit, daß es weiterwächst. 
Die Wahrscheinlichkeit, daß es trotzdem den kritischen Radius erreicht, 
kann man aus dem BoLtzmann-Prinzip zu 


W=exp dk 
: : : : i dQ 
ermitteln. Da sich andererseits die Entropie ® aus ® — Eh m berech- 


net, erhält man für konstante Temperatur 

A 

T >} 

wobei A die Arbeit zur Bildung des Tröpfchens vom Radius r; ist. Das 
negative Vorzeichen steht deshalb, weil bei der Tröpfehenbildung Arbeit 


DD — — 


*) Für das zweidimensionale Problem sind neuerdings auch von RunGE (56) einige 
besonders die Kinetik im Auge behaltenden Rechnungen ausgeführt worden. Die An- 
sätze folgen der gleichen Linie wie die Reihenentwicklungen von VAN DER WAERDEN 


(57) für das Ordnungsproblem. Die Berechnungen gelten aber nur für den geometrisch | 
ungestörten Kristall. 
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geleistet werden muß. A ist also positiv zu nehmen. Damit erhalten wir für 
die Wahrscheinlichkeit der Bildung eines Tröpfehens von Radius r; 


(IT, 91) W=exp— A/KT. 

Für die Anzahl n der pro sec. gebildeten Tröpfehen gewinnt man die 
Beziehung 

(II, 92) n=Kexp— A/kT. 


Die Berechnung, auf die hier nicht eingegangen wird, ergibt, daß K 
gleich der mittleren Anzahl der Zusammenstöße zweier Teilchen pro Sekunde 
ist. Die Arbeit A zur Bildung eines beliebigen Tröpfchens kann man auch 
ausrechnen. Sie beträgt 


4 
(II, 93) sau indie, 


Das erste Glied entspricht den (temperaturabhängigen!) Volumen- und 


‚ Oberflächenanteilen zur freien Energie. Den kritischen Radius kann man 


oa 
aus Gleichung (II, 93) mittels der Beziehung Fr 0 berechnen und erhält 


(II, 94) BR, 

® 

Ganz analog kann man nun für den Kristall verfahren. Hier macht 
BECKER den Gleichung (II, 92) entsprechenden Ansatz 


(II, 95) n = Cexp — Q/kTexp— A/KT. 


An Stelle von K wurde hier © exp — Q/KT eingeführt. Die Notwendigkeit 
ist klar, denn im Kristall kommt für die Diffusion die Schwierigkeit der 
Überwindung der Potentialschwelle für einen Platzwechsel hinzu. Q ist also 
die Aktivierungsenergie für einen Platzwechsel, die aus der Temperatur- 
abhängigkeit der Diffusionskonstanten ermittelt werden kann. A ist wieder 
die Keimbildungsarbeit, die für ein würfelförmiges Entmischungsgebiet der 
Kantenlänge a 
(IL, 96) A=—a?v+6aro ist. 

Setzen wir v und o als bekannt voraus, so berechnet man a, wieder aus 


oA 
—=() an 
03 


4o 
(IL, 97) Al 


v 


co und v können nun z. B. nach dem Wechselwirkungsprinzip der nächsten 
Nachbarn ausgerechnet werden, das wurde von BECKER näherungsweise 
durchgeführt und man gewinnt damit quantitative Vergleichsmöglichkeiten 
mit den experimentellen Ergebnissen. Wir wollen jedoch darauf nicht näher 
eingehen, sondern vielmehr das Ausscheidungsproblem von der kinetischen 
Seite her, insbesondere bei auftretenden Schwierigkeiten in der Gitter- 
geometrie näher beleuchten. 
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Der kritische Radius ist stark temperaturabhängig, das tritt zwar in 
den Gleichungen (II, 94) und (II, 97) nicht deutlich in Erscheinung, aber 
diese Tatsache kann man sich am Kondensationsproblem besonders gut 
klar machen. Der Dampfdruck eines Tröpfchens steigt mit höherer Tem- 
peratur und daher ergibt sich auch sofort ein größerer kritischer Radius. 
Die in Gleichungen (II, 94) und (II, 97) eingehende Größen v und o sind eben 
nicht unabhängig von der Temperatur. Bei sehr tiefen Temperaturen ist der 
kritische Radius sehr klein und in diesem Fall wird die Ausscheidung mehr 
durch die Diffusionsschwierigkeiten als durch den kritischen Radius be- 
hindert. 

Betrachten wir zunächst einmal irgendeinen Vorgang bei konstanter 
Temperatur. Der BECKERr’sche Ansatz hat den Vorzug, daß man auch etwaige 
Gitterverzerrungen, die mit der Entmischung verknüpft sein können, mit 
einarbeiten kann. Allerdings ergibt das Anteile, die nicht allein von der 
Oberfläche des Keimes, sondern auch von dessen Größe abhängig sind. 
Formal kann man den Ansatz (LI, 96) in dem Sinne erweitern, daß man zur 
Öberflächenenergie noch einen allerdings a-abhängigen Teil hinzunimmt, 
und danach den kritischen Radius berechnet. Schreiben wir also versuchs- 
weise 


a? 
A=— +62 [+ ol], 


so erhalten wir nunmehr für die kritische Kantenlänge die Beziehung 
(IL, 98) [w—20(a)]=4[c + (a). 


Gleichung (II, 98) gibt nun guten Aufschluß darüber, was die Hinzu- 
nahme der Gitterspannungen an neuen Ergebnissen bringen kann. Ist die 
Gitterverzerrung nur von der Oberfläche und nicht vom Volumen abhängig, 
so ist 

o’(a) =, 


d.h. an unserem Ergebnis ändert sich nichts Grundsätzliches, die kritische 
Kantenlänge a; wird lediglich zu größeren Werten hin verschoben. In Fig. 8 
wurden die linke und rechte Seite von Gleichung (II, 98) für drei Tempera- 
turen aufgetragen. Wir nehmen hier an, wie es auch von BECKER getan 
wird, daß der Oberflächenanteil zur freien Energie temperaturunabhängig 
ist, das ist natürlich nicht exakt. Die folgenden nur qualitativen Betrachtun- 
gen werden aber auch durch eine Verallgemeinerung in dieser Richtung im 
Prinzip nicht geändert. : 

Für o (a) = const, o’(a) = 0 ist die rechte Seite von (II, 98) eine 
horizontale Gerade, und die linke Seite eine Gerade mit einer v entsprechen- 
den Neigung, die mit sinkender Temperatur immer steiler wird ; der Schnitt- 
punkt beider Kurven entspricht einem Maximum der für die Keimbildung 
zu leistenden Arbeit und ergibt den Wert für die kritische Kantenlänge. 

Ist die Verzerrungsarbeit des Gitters aber stark abhängig vom Volumen 
des Keimes, so werden die Verhältnisse wesentlich komplizierter. Ohne eine 
bestimmte Vorstellung über die Art der Gitterspannungen kann man 
natürlich zunächst nichts Quantitatives aussagen. Aber man kann leicht 
einsehen, daß nicht nur ein Schnittpunkt, sondern mehrere Schnittpunkte 
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der beiden Kurven möglich werden. Daß die Spannungsarbeit bei stark 
abweichender Gittergeometrie des Auszuscheidenden sehr stark anwachsen 
muß, ist klar, denn die Nichtvergleichbarkeit der Gitterabstände wird 
schon bei kleiner Anzahl der Atombesetzung des Keimes sehr schnell er- 
hebliche Spannungen herbeiführen. Das gilt vor allen Dingen für kleine 


y (Sta)+e) 


“pn 0° arg s 
Fig. 8. Kritische Kantenlänge für verschiedene Temperaturen T!> T, > T,,. 
Keime mit nur oberflächenabhängiger Verzerrungsarbeit im Gitter. 


Keime. Für größere ist es schwerer zu übersehen; aber immerhin wird man 
annehmen können, daß bei großen Keimen eine Entspannung durch Umbau 
im Keiminneren auftreten kann, so daß schließlich die Spannungsenergie 
eines Keimes nur noch von der Oberfläche abhängt, vorausgesetzt, daß noch 
kein Bruch mit dem Wirtgitter eingetreten ist. In Fig. 9 wurde nun ein den 


Fig. 9. Angenommener Verlauf von o (a) ; Spannungen im Gitter steigen stark 
an, sind bei genügender Keimgröße aber nur noch oberflächenabhängig. 
v—2 o’ (a) -- - wurde in Abhängigkeit von a für 3 Temperaturen ausgerechnet. 


vorangegangenen Überlegungen entsprechender Ansatz für die Verzerrungs- 
arbeit in Abhängigkeit von der Keimgröße ausgeführt, dem natürlich nur 
ein qualitativer Wert zuzusprechen ist. Desgleichen wurde die Ableitung 
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dieser Funktion eingetragen. Die resultierenden graphischen Darstellungen 
der beiden Seiten von Gleichung (LI, 98) ist in Fig. 10 für drei Temperaturen 
durchgeführt. Man erkennt, daß unter dieser Voraussetzung bei der 
Temperatur T, insgesamt drei Schnittpunkte auftreten, von dem der erste 
und dritte einem Maximum, der zweite aber einem Minimum entspricht. 
Wir haben in diesem hypothetischen Fall zwei kritische Keimgrößen ax, 


und a,,. Es kann nun durchaus der Fall sein, daß die Aktivierungsenergie 
für die Bildung des Keimes ax, so hoch liegt, daß er praktisch nicht auf- 
treten kann. Ist dagegen die Entstehung des Keimes a,, möglich, dann 


gibt es für die Stabilität dieses Keimes eine bestimmte bevorzugte Größe. 
D. h. der Entmischungsvorgang verläuft so, daß gar keine statistische Ver- 
teilung von Keimen verschiedener Größe auftritt, sondern sich ein Gleich- 
gewicht von Keimen von ziemlich gleichmäßiger Größe ausbildet, und damit 
eine relativ gleichmäßige Verzerrung des gesamten Gitters erfolgt. Wird die 
Temperatur erhöht (T, in Fig. 10), dann gibt es schließlich nur einen Schnitt- 
punkt der beiden Kurven, d. h. die gebildeten Keime lösen sich alle wieder 
auf (Rückbildung), vorausgesetzt, daß die Aktivierungsenergie für die Aus- 
bildung der Keime ax, immer noch zu hoch liegt. 


Eine weitere Erhöhung der Temperatur muß schließlich einmal zur 
Ausbildung der großen nunmehr endgültig wachstumsfähigen Keime führen. 
Es ist klar, daß diese rein qualitativen Überlegungen erst durch exakte 
Berechnungen erhärtet werden können. Der hier vorgenommene Ansatz 
ist natürlich schlecht und wurde nur ausgeführt, um Licht auf das grund- 
sätzliche Verhalten des Entmischungsablaufs beim Auftreten starker Gitter- 
spannungen zu werfen. Die Annahme eines würfelförmigen Keimes und der 
Ansatz für die Verzerrungsarbeit sind hier natürlich nur als rein formal zu 
betrachten. In Wirklichkeit wird der Keim wegen der Anisotropie der Span- 
nungen nicht würfelförmig sein. Es ist leicht einzusehen, daß eine solche 
Anisotropie meist zu blättchenförmigen Ausscheidungen führen muß, da 
das Wachstum der Keime bevorzugt in den Richtungen erfolgt, in denen 
die Verzerrungsarbeit am geringsten ist. Ein solches Verhalten wurde in der 
Tat bei den meisten Entmischungen gefunden. 

Es ist klar, daß die drei Schnittpunkte nur unter ganz bestimmten 
Bedingungen auftreten können. Tritt kein Spannungsabbau ein, so sind 
größere Keime bei der Annahme stärker anwachsender Verzerrungsarbeit 
nicht möglich. Das Verhalten wird wesentlich durch die Differenz von 
Spannungsarbeit und Energiegewinn durch die Entmischung bestimmt. Wir 
wollen hier diese an sich einfachen Verhältnisse nicht weiter diskutieren, 
weil das ohne quantitative Ansätze zu spekulativ ist. Aber immerhin besitzen. 
wir z. B. im Falle des Duraluminiums experimentelle Hinweise dafür, daß 
die Verhältnisse beim Entmischungsvorgang in der besprochenen Richtung 
liegen müssen. Hier tritt bei Zimmertemperatur eine Keimbildung auf, die 
sich sehr schnell festfährt und zu einer gleichmäßigen Verspannung des 
Gitters führt, die eine ungewöhnliche Härte des Kristalls bedingt (45). Bei 
höheren Temperaturen tritt nun in der Tat Rückbildung dieser Keime ein, 
wobei gleichzeitig größere Keime wahrscheinlich bereits unter teilweisem 
Spannungsabbau entstehen. | 
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Fig. 10. Graphische Darstellung von a [v—2 0° (a)] — für 3 verschiedene Temperaturen 

und 4 [0 + o (a)] --; die Schnittpunkte der beiden Funktionen ergeben den kritischen 

Radius. Für T, gibt es 3 Schnittpunkte, der erste und dritte entspricht einem Maximum, 
der zweite einem Minimum der Keimbildungsarbeit. T, > T, > T;. 
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Wir wollen damit das Gebiet der Ausscheidung, das heute noch in 
vielen Punkten ungeklärt ist, verlassen und uns einer Theorie zuwenden, 
die eigentlich mehr vom theoretischen Standpunkt von Interesse ist. 


II, 8. Die gruppentheoretische Methode der Fehlordnung. 


Die gruppentheoretische Methode der Fehlordnung hat sich insofern als 
besonders leistungsfähig erwiesen, als durch ihre Anwendung die Lösung 
des zweidimensionalen Fehlordnungsproblems exakt gelungen ist; das ist 
deshalb besonders interessant, als sich hier der Eintritt der Fernordnung 
zwangsläufig aus der Theorie ergibt, und nicht durch einen Ansatz a priori 
eingeführt werden muß. Beim eindimensionalen Modell tritt die Fernordnung 
noch nicht auf, weil diese ja erst durch die höhere Koordination im Netz 
oder Gitter hervorgerufen werden können. Vom Standpunkt der experi- 
mentellen Untersuchungen hat aber das zweidimensionale Fehlordnungs- 
modell noch keine Bedeutung erlangen können, weil zweidimensionale Fehl- 
ordnung noch nicht häufig beobachtet wurde. 


Die gruppentheoretische Methode wurde durch KRAMERS und WANNTER 
(46) für die Behandlung eines Modells des zweidimensionalen Ferromagneti- 
kums durchgeführt. Allerdings gelang diesen Autoren die exakte Lösung noch 
nicht; MONTROLL (47) hat kurz danach die prinzipielle Anwendung für 
dreidimensionale fehlgeordnete Kristalle diskutiert, konnte aber, wie auch 
KRAMERS und WANNIER, nur Näherungslösungen angeben. Weitere Ent- 
wicklung auf diesem Gebiet verdanken wir LASSETTRE und HowE (48), 
MOoNTROLL (49), AsHKkın und LAmB (50) sowie AsHKIN und TELLER (Sl). 
Während sich jedoch alle bisher genannten Autoren mit Näherungslösungen 
zufrieden geben mußten, gelang ONSAGER (52) die exakte Durchführung der 
Berechnung der Zustandssumme für das zweidimensionale Problem. Aller- 
dings ist der dazu notwendige mathematische Apparat außerordentlich 
kompliziert und wir können hier aus diesem Grunde nicht die vollständige 
Durchführung der Lösung bringen, aber es sollen wenigstens die Hauptzüge 
der Theorie dargestellt werden. Ehe wir jedoch dazu übergehen, wollen wir 
für die mathematisch weniger geschulten Leser einige Elemente der Matrizen- 
rechnung, so weit sie hier gebraucht werden, kurz durchsprechen. 


Wir betrachten irgendeine beliebige Matrix A der Ordnung n, deren 
Elemente der i-ten und k-ten Spalte mit A;, bezeichnet werden. Aus- 
geschrieben sieht dann die Matrix also folgendermaßen aus 


| Anı Ana Ann | 


Man kann jeder Matrix Eigenwerte zuordnen, die durch die Gleichung 
(IT, 99) |A—71|=0 


gegeben sind. Die Betragsstriche sollen dabei andeuten, daß die Determinante 
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der resultierenden Matrix gemeint ist. Mit 1 wurde die Einheitsmatrix vom 
Rang n bezeichnet. Diese hat nur in der Hauptdiagonalen von 0 verschiedene 
Elemente der Größe 1. A ist eine skalare Größe; Gleichung (IL, 99) ist all- 
gemein als die sogenannte Säkulargleichung bekannt. Sie repräsentiert 
offenbar eine Gleichung n-ten Grades für die unbekannte Größe . Die 
Lösungen ließen sich sofort angeben, wenn die Matrix A eine Diagonal- 
matrix wäre, also nur auf der Hauptdiagonalen von 0 verschiedene Elemente 
hätte. Wir bezeichnen die n Lösungen der Gleichung (II, 99) mit 2,, A, 


u. ...:.. ‚ An. (1,99) kann aber auch als Bestimmungsgleichung der 
nn-dimensionalen Vektoren p; j =1,......... ‚n) aufgefaßt werden, deren 
Komponenten wir mit pıj, Pay» ---- ---- ‚ Pnj bezeichnen wollen. Da jedes 


4,(5=]1,....,n) der Gleichung (Il, 99) genügt, wird offenbar auch das 
Gleichungssystem 


(II, 100) A,,Ppız + ApP2jt+ --- + Aın Pnj = 4 Pıj» 
AyıPıj + Age Pe; + --- + Azn Pni = A Pa 
AnıPpız  Anspo; 4... 4 Ayn Prnj = A Pnj 


als Bestimmungsgleichung des Vektors p; erfüllt. Die so bestimmten Vek- 
toren nennt man auch die Eigenvektoren der Matrix A. Ein beliebiger 
konstanter Faktor, der das Gleichungssystem nicht ändert, bleibt dabei 


für jedes p; noch offen. In der üblichen Vektorenschreibweise kann man 
‘ für das Gleichungssystem (II, 100) auch schreiben 


(I, 101) Ap=7;p. 


' Diese Schreibweise wird unmittelbar klar, wenn man für den Vektor p; die 


zulässige Matrixstellung 
[Pi 


Ze 

a 

= j | 
WW, 


| Pnj 


' wählt und die Matrixmultiplikationsregel (AB)ıx = 2 A:; B;x anwendet. 


j 
(AB);« ist dabei das Element der i-ten Zeile und der k-ten Spalte des 
Matrixproduktes (AB). Die Elementdarstellung der Gleichung bedeutet 
dann 


(IL, 102) & Ak Pri =APii- 
k 


Das Ergebnis ist also mit dem Gleichungssystem (IL, 100) identisch. 
Wir bilden nun aus den n Vektoren p; eine Matrix derart, daß die i-te Kom- 
ponente des Vektors p; im Kreuzungspunkt der i-ten Zeile und j-ten Spalte 
steht. Wir erhalten damit eine Matrix P und wollen zunächst den Vektor- 
charakter ihrer Spalten nieht im Auge behalten. Da nun aber in Gleichung 


Fortschritte der Mineralogie 1949. Bd. 28, 10 
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(IT, 102) das Element mit den Indizes ij des Matrixproduktes (APP) steht, 
wird also die Gleichung (II, 102) zu 


(AP); = A pi 


oder 
[ 4, 0 
(IT, 103) (AP)=(PD@))mitD() =| °. 
Ki "m 
Rechtsseitige Multiplikation mit P-! ergibt 


(IT, 104) A=(PD () P-). 


Die Umformung (II, 104) setzt allerdings voraus, daß P eine Reziproke 
P-! besitzt; dies ist aber bei linearer Unabhängigkeit der Vektoren p; 
immer erfüllt. Die Elemente A;, der Matrix A lassen sich nun in einfacher 
Weise aus den Elementen (P;, und (P-!),, ermitteln. Es ist 


Ax = (PDWAP")« =zP;(D (AP), 
J 
(IT, 105) (DA)P")x = (Pix, 
Axr=24 P; (Pr. 
J 


Ist die Matrix P unitär, das kann aber immer erreicht werden, wenn A 
hermiteisch oder unitär ist, weil im Gleichungssystem (II, 100) noch ein 
konstanter Faktor offenbleibt, so gilt die Beziehung 


(Pr) =P*; 


(ist die Matrix P reell, so können wir das Sternzeichen für den konjugiert 
komplexen Wert fortlassen). 


Wir erhalten also für eine hermiteische Matrix A 

(II, 106) Ak 2), P;; Pr, =2% Pij P*rj: 
j i 

Nach diesen einleitenden Betrachtungen aus der Theorie der Matrizen 
wollen wir unsnun unserem Fehlordnungsproblem zuwenden. Wir betrachten 
den Fall eines‘ dreidimensionalen einfachen kubischen Gitters, dessen 
1, x 1, x 1; Gitterplätze wir mit einer beliebigen Anzahl von A- und B- 
Atomen zu besetzen haben (l, entspricht der Anzahl der Zellen in der 
Translationsrichtung a,). Wir bauen zu diesem Zweck den Kristall suk- 
zessive aus den 1,-Netzebenen mit 1, x 1,-Gitterplätzen auf. Zur Verein- 
fachung der Ergebnisse fordern wir, daß die Netzebene 1, + 1 mit der 
Netzebene 1 identisch ist (d. h. also der Kristall ist in der vierten Dimen- 
sion in der Weise geschlossen, daß die letzte und die erste Netzebene 
aufeinander liegen). Diese Einschränkung bedeutet keine Beschränkung der 
Allgemeinheit unserer Rechnung, solange 1, sehr groß ist. Aus Gründen der 
Vereinfachung unserer Überlegungen wählen wir für eine A-Besetzung 
irgendeines Gitterplatzes das Symbol + 1 und für eine B-Besetzung das 
Symbol — 1. Da jeder der Gitterplätze mit A oder B besetzt werden kann, 
gibt es insgesamt 2%, für eine einzige Netzebene aber 2: verschiedene 
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Konfigurationen. Für die i-te Konfiguration (der 2") in der Netzebene j 
schreiben wir hier u; (i=]1,...,2,";j=]1,...,1,). Wir wollen uns weiter- 
hin nur auf die Wechselwirkungsenergien der nächsten Nachbarn beschrän- 
ken. Eine Verallgemeinerung in dieser Richtung stößt zwar nicht auf 
grundsätzliche Schwierigkeiten, kompliziert aber die Ansätze sehr. Zur 
Ermittlung der Zustandssumme des Kristalls müssen wir nun die Wechsel- 
wirkungsanteile der einzelnen Nachbarn abzählen. Wir setzen hier 
Vu=YHtHH)=-V\s=V-—)=V 
und Vs=-V-+)=- Vya=-V--+4=-—V\, 
: dann ergibt sich der BorLrzmann-Faktor für den Anbau der Netzebenej + 1 
; an die Netzebene j aus der Abzählung der Bindungen innerhalb der Netz- 
ebene sowie derjenigen Bindungen, die zwischen den Netzebenen j undj + 1 
liegen. Den ersten Anteil bezeichnen wir mit V (ui;), den zweiten mit 
V (wj, Aijzı), wobei beachtet werden muß, daß die Verwendung des 
gleichen Index innerhalb der beiden Netzebenen nicht die gleiche Kon- 
ı figuration in den Netzebenen j und j + 1 bedeuten muß. Der BoLTzMmanN- 
' Faktor ist 
1 


exp — [V(u;) + V(ni, mir). 


| Der Borrzuasx-Faktor für eine Konfiguration des gesamten Kristalls beträgt 
‘ daher 


EN a 
N il 

| (II, 107) exp | —_ - > (Vla;) + Va, Mij+1)) I 
ner 5 


Man erhält die Zustandssumme, wenn man alle Netzebenen unab- 
: hängig alle 2!" Konfigurationen durchlaufen läßt 


| BT, 108)2 = > >> R >> exp - iD (Vu) + V(u De): 


EiTE2i2 Mil, 
% or u 
| er 2 > N > I a b er! Vers) + Ve rij4n) 
Kiı #i2 a), 3-1 L e| 


‚ Schreiben wir statt 
Via) , Vlan+ı) | 


= | 
2 


2 


V (u; ij+1) = v (mi, Mj+n): 

so wird 

. l; 

k Al, 109) IE oe II EXP —V (ui, ij an ı)/KT, 

14142 Aa, Il 

da l, + 1=1 ist. Alle möglichen Werte von v (wjj, 41; +1) sind nun offenbar 
;durch die 21% x 2ll möglichen Konfigurationen der beiden Netzebenen 
‘gegeben. Man kann diesen Werten eine Matrix vom Rang 2! zuordnen. 
. Der Einfachheit wegen lassen wir die Indizes, die sich auf die Netzebenen 
'j, j + 1 beziehen, fort. Mit v;x meinen wir also, daß sich die Netzebene j 


10* 
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in der i-ten Konfiguration und die 'Netzebene j + 1 in der k-ten Kon- 
figuration befindet. Wir bilden also die Matrix 


exp — v,,/KT exp — vj5/KT ... exp — vja11/kT 
exp — v,,/kT exp — vy/kT ... exp — vg51u1,/KT 


(II, 110) —E 


| EXP — Val), Fr EeXP — Yal!, Se + EXP — Voll, gu1,/KT | 


Die Matrix ist aber symmetrisch und reell, denn es ist gleichgültig, 
ob sich die Netzebene j in der Konfiguration i und die Netzebene j + lin 
der Konfiguration k befindet oder ob es gerade umgekehrt ist. Es ist also 
vr = vi; V ist also auch hermiteisch. Damit sind auch ihre 21 Eigen- 
vektoren reell und orthogonal, und wir können zur Berechnung des Elements 
V;x Gleichung (IL, 106) verwenden. Wir erhalten 

olıl, 


(LI, 111) VYir- a Pi» Pk» 


vu 


wobei also }, der v-te Eigenwert der Matrix VW ist. 


Die eben erwähnten Orthogonalitätsbedingungen für die Vektoren p; 
lauten wegen der Unitarität von P 


olıle ij —( füri=s 1; 
(II, 112) (Pi, pj) = Pri Prj = Öj 5 Öj — ] füri —aR 
veill 


Setzen wir nun Gleichung (II, 111) in die Zustandssumme ein, so er- 
halten wir unter Berücksichtigung der vorgenommenen Umindizierung 
% alı ls 
(1, 113) Zee el 


Av Pay,» Prajr,o ® 
2 42 a 


Die Auswertung der Gleichung (II, 113) kann mit Hilfe der Orthogonali- 
tätsbedingungen (II, 112) vorgenommen werden. Um das Verfahren besser 
übersehen zu können, schreiben wir das Glied II explizit auf. Die etwas 


v 
umständlichen Indizes, die sich auf die verschiedenen Netzebenen beziehen, 
vereinfachen wir, müssen aber im Auge behalten, daß die 1. Indizes der p 
sich auf eine Netzebene beziehen, die noch über alle Konfigurationen 
summiert werden müssen. Wir schreiben also 


Hz = (A, Pıı Psı + 4 Pı2 Pag + As Pıs Pas + --- + Aglıl Pıahls Paglıl), 
X (Al Paı Psı + Ag Pa2Ps32 + As Pas Pas +: -- + Agıl, Paelıl, Pzell.), 


X (A Prı Pıı + Ag Pıa Pız + As Pys Pıs + -- + Aglıl, Pielıl, Pıall.)- 
In diesem Produkt betrachten wir zunächst die Glieder 


Al Pı» P»» P»» Ps» - - - Pı» Pıv- 
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Der erste Index gibt die Konfiguration der Netzebenen 1, ...,], an, 
‚und da jedes einzelne Glied des obenstehenden Produktes aber alle mög- 
lichen Konfigurationen unabhängig zu summieren ist, können wir diese 
Summierung schrittweise ausführen. Wir summieren also für die j-te Netz- 


ebene und erhalten 
all, 
A P1» Pav- .. | I Da Bee 
j=o 
Vergleich mit den Orthogonalitätsbedingungen zeigt, daß die Summe den 
Wert 1 ergibt. Das gleiche erhält man für alle Paare mit gleichem 1. Index, 
wenn man über alle Konfigurationen summiert. Also wird das betrachtete 
olls 
Glied gleich Ab. Alle Glieder zusammen ergeben den Beitrag x A". 
vi 
Für alle gemischten Glieder ist es nun leicht nachzuweisen, daß immer 
mindestens ein Paar von Vektorkomponenten Pr, Pr, d. h. mit gleichen 
ersten Indizes aber ungleichen 2. Indizes auftreten muß. Da der erste Index 
die Konfiguration der Netzebene r darstellt, ergibt die Summierung über 


diese Netzebene 
olıl 
x| > Pr»Pru x 


v=1 


und man erkennt aus den Örthogonalitätsbedingungen, daß die Summe, 
also auch das ganze Glied verschwindet, damit haben wir den Wert der 
‚ Zustandssumme errechnet und erhalten 
BR 
(IL, 114) ZEN 
v=1 

; Mit (IL, 114) ist die Berechnung der Zustandssumme auf die Auffindung 
' der Eigenwerte der Matrix V zurückgeführt worden. Gibt es einen größten 
Eigenwert, so kann man, vorausgesetzt, daß 1, eine sehr große Zahl ist, die 
Beiträge der kleinen Eigenwerte vernachlässigen und nur den größten 
; Eigenwert für die Zustandssumme benutzen. Es muß dann natürlich die 
- Bedingung erfüllt sein, daß die Summe der l,-ten Potenzen der Eigenwerte 
| klein im Vergleich zu De ist. In diesem Fall gibt also 

(I, 115) IN 
"Damit haben wir aber noch nichts gewonnen, wenn wir nicht eine Möglich- 
keit finden, den größten Eigenwert der Matrix V zu berechnen. Das scheint 

nun fast eine unlösbare Aufgabe zu sein. Aber gewisse Eigenschaften kann 
man trotzdem aus der Matrixdarstellung entnehmen. Eine kritische Tem- 
peratur muß mit einer Entartung des größten Eigenwerts der Matrix V/ 
verknüpft sein. KRAMERS und WANNIER haben gezeigt, daß man diese 
Tatsache dazu benutzen kann, die kritische Temperatur festzulegen. Sie 
transformieren die Matrix V mit einer Ähnlichkeitstransformation und 
gewinnen damit eine neue Matrix, welche die Eigenschaft besitzt, daß sich 
ihr Eigenwertspektrum gegenüber der ersten Matrix nicht ändert. Gibt es 
eine Entartung des größten Eigenwerts, die bei einer bestimmten Temperatur 
in der ersten Matrix eintritt, so muß die Entartung in der anderen Matrix 
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bei der gleichen Temperatur stattfinden. Diese Überlegungen führten 
KRAMERS und WANNIER, ohne daß ihnen die Lösung des Eigenwertproblems 
selbst gelang, zu der richtigen Festlegung der kritischen Temperatur für das 
zweidimensionale Modell, vorausgesetzt, daß eine kritische Temperatur über- 
haupt existiert. Mit Hilfe der Variationsrechnung können diese Autoren 
bereits eine Näherungslösung des zweidimensionalen Problems geben, die 
recht gut sein dürfte. 

AsHKIN und Lamg (50) konnten zeigen, daß die Entartung des größten 
Eigenwerts nicht — wie KRAMERS und WANNIER annahmen — nur für 
eine bestimmte Temperatur vorhanden ist, sondern in einem mindestens 
in zwei Dimensionen unendlichen Kristall über einen ganzen Temperatur- 
bereich existiert. Und zwar reicht diese Entartung von T = 0 bis zur kri- 
tischen Temperatur T., liegt also in dem Bereich der Fernordnung im Kristall. 
Man kann nachweisen, daß im Fall der Nichtentartung des größten Bigen- 
werts die Besetzung irgendeines Gitterplatzes unabhängig von der Besetzung 
aller Gitterplätze ist, die von dem betrachteten Gitterplatz sehr weit entfernt 
sind. Das ist aber gerade das Charakteristische für fehlende Fernordnung. Ist 
Amax dagegen entartet, so sind diese Wahrscheinlichkeiten nicht mehr un- 
abhängig voneinander. Von MoNTROoLL (49) wurde bewiesen, daß gemäß 
einem Theorem von FROBENIUS, demzufolge der größte Eigenwert einer 
Matrix dann positiv und nicht entartet ist, wenn alle Elemente der Matrix | 
positiv sind, für den endlichen Kristall keine Fernordnung zu erwarten ist, | 
weil ja die in der Matrix V stehenden Glieder alle reell und positiv sind. | 
Erst wenn der Kristall in zwei Dimensionen unendlich wird, kann Amax einer 
Entartung sehr nahe kommen, weil dann einige Elemente praktisch 0 werden | 
können. — Weitere Näherungslösungen kann man auch über die Störungs- 
rechnung der Matritzentheorie gewinnen. Wir brauchen hier nicht darauf | 
einzugehen, da inzwischen von L. ONSAGER (52) die exakte Berechnung des 
größten Eigenwertes ausgeführt wurde. 

ONSAGER konnte seine Berechnung auf einen es: Kristall aus- : 
dehnen, wir wollen seine Methode wenigstens in ihrem Gedankengang er- : 
läutern. Das den Berechnungen zugrunde gelegte Modell ist in Fig. 11 


r => ES 
1% 
\ 
[6) 
eo Bi 
og 
ee 
IN | —__g__0— |,* 
ua Su | | | _x Ix 
rl v=x | zn x 
au x Dr | FE | 
x Sy | | rn |x 
3 Be Kar x 
ER Ri | v2 
2 | ‘o 


Fig. 11. Onsager’s Modell für den zweidimensionalen rhombischen Kristall. 
Anbau der Kette O erfordert Kenntnis der Konfiguration nicht nur dieser Kette, | 
sondern auch der Kette @. 
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wiedergegeben. Wir teilen diesmal das Modell in ringförmig geschlossene 
Ketten ein, die wir, wie im dreidimensionalen Fall schrittweise aneinander- 
fügen. Die ringförmige Schließung der Ketten geschieht hier wiederum nur 
aus Gründen der mathematischen Einfachheit der Darstellung. Ebenso 
wollen wir annehmen, daß die letzte Kette unmittelbar mit der ersten be- 
nachbart ist, auch das steht in voller Analogie zum Lösungsverfahren des 
dreidimensionalen Modells. Die Wechselwirkungsenergien seien hier, wegen 
der verschiedenen Gitterabstände innerhalb der Ketten und senkrecht dazu 
- mit V’ und V bezeichnet. Zwei gleiche Nachbarn innerhalb einer Kette 
haben also die Wechselwirkungsenergie + V’, zwei ungleiche dagegen — V'. 
Für die Nachbarn von Kette zu Kette betragen die entsprechenden Energien 
+ V und — V. Nennen wir noch V/kT—=— H und V’/kT=-—-H/, so sind 
die BoLrzmann Faktoren für ein gleiches Paar exp H bzw. exp H’ und für 
ein ungleiches Paar exp — H’ bzw. exp —H. 


Eine Konfiguration der Kette j, die aus n Gliedern bestehen möge, 
ist also durch Angabe von u; = (1, Us, - - -, in) gegeben. 11, - - -, Un können 
also, wie vorhin definiert, die Werte + 1 annehmen. Wir nehmen nun an, 
daß die relative Wahrscheinlichkeit, die Kette j in der Konfiguration 
8... ‚4n) zu finden, bekannt sei (=P (w,, -:--,4n)). Wir können 
dann berechnen, wie groß die relative Wahrscheinlichkeit ist, die Kettej + 1 
in der Konfiguration (w,, - - -, {'n) zu finden. Offenbar besteht die Beziehung 


(II, 116) 


n n 
,P (un sun) = Plün ---, Mn) opfirz a ee 7 Yu) 
Ne ji=1 j-1 

Der Faktor A muß hier erscheinen, weil es sich nicht um die abso- 
luten, sondern um die relativen Wahrscheinlichkeiten handeln soll. Die 
erste Summe im Exponentialglied zählt dabei die Bindungspaare zwischen 
der j-ten und der vorhergehenden Kette, die zweite Summe dagegen die 
Bindungen innerhalb der neu angebauten Kette. Gleichung (II, 116) ist 
weiter nichts als die zu Gleichungen (II, 11, 12) analoge Formel für das zwei- 
dimensionale Problem. Im eindimensionalen Fall läßt sich die Gleichung 
lösen, während hier wegen der 2" verschiedenen Vektorkomponenten eine 
Lösung ausgeschlossen erscheint. Betrachtet man wie im eindimensionalen 
Modell die P (u, -- -, 4n) als Komponenten eines Vektors, so erkennt man 
sofort, daß die Gleichung (II, 116) eine schr nahe Verwandtschaft zu unserem 
Matrixeigenwertproblem (Gleichungen (II, 100)—(IL, 102)) besitzt. Der 
einzige Unterschied liegt in der Verwendung der verschiedenen u; und u’. 
Wenn wir aber bedenken, daß der physikalische Sachverhalt bei großer 
Kettenzahl für alle Ketten die gleiche ist, fällt auch diese Einschränkung 
und wir haben also wieder die Lösung des Eigenwertproblems, der durch 
Gleichung (II, 116) definierten Matrix V vorzunehmen. 

Wir wollen nun die Matrix explizit aufschreiben. Sie besitzt offenbar 
den Rang 2". Wir vereinfachen zu diesem Zweck die Bezeichnungsweise. 
Die ; und die w'; durchlaufen ja unabhängig voneinander alle möglichen 
Konfigurationen, wir kennzeichnen nun die erste Summe im Exponential- 
glied mit zwei Indizes, wobei der erste die Konfiguration der u’; und der 
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zweite die Konfiguration der u; angibt. Die zweite Summe im Exponential- 
glied braucht, da ja nur eine Konfiguration eingeht, nur einen Index, 
dieser muß aber mit dem ersten der ersten Summe übereinstimmen. Damit 
erhalten wir 


(II, 117) 


[expH 2,] expH' 2, expH 2,expH! 2,...ezpH 2.2 expH7 25 
expH 2,, expH" 2, expH 2, expHl2,...exp$ Hr 2, rer pH 


\\f= 
exp H-2,% expHl I. en en expH 2,",nexpH’ 2,2 
Die Matrix V kann nun aber als Produkt geschrieben werden 
[exp H' 2, 0 ] Ki Dura exp | 
V= 
(0) . 
| exp H’ 2,n exp H 2.n,. .. exp Hera 


Damit haben wir die Beziehung gewonnen 


Man kann also den Anbau in zwei Teiloperationen zerlegen (s. Fig. 11). 
V, entspricht den Wechselwirkungen zwischen den Ketten und VW, den 
Wechselwirkungen innerhalb der neuen Kette. Der Grundgedanke On- 
SAGER’S ist, nicht die Matrizen selbst zur Berechnung heanzuziehen, son- 
dern neue Operatoren einzuführen, deren Eigenschaften leicht zu erkennen 
sind. Mit diesen Operatoren baut ONSAGER eine abstrakte Algebra auf, mit 
deren Hilfe sich die Lösung des Eigenwertproblems sehr viel einfacher finden 
läßt, als mit den unhandlichen Matrizen. 


Gemäß Gleichung (II, 116) müssen wir die Wirkungsweise der Matrix V/ 
auf die Eigenvektoren studieren. Da die Matrix V, zuerst angewendet 
werden muß, wollen wir einen entsprechenden Operator zu bestimmen 
suchen. Wir gehen dazu von Gleichung (II, 116) aus, und zwar greifen wir 
die Konfiguration u; —= w'; heraus. Für dieses Glied hat die erste Summe 
offenbar nur positive Glieder und hat deshalb den Wert n. Ändert sich nun 
in den u; ein Vorzeichen, so bedingt das ein Glied mit negativem Vorzeichen 
in der Summe und Anderung des entsprechenden Vorzeichens im Vektor. 
Wir führen deshalb einen Spiegeloperator c; ein, der die Eigenschaft besitzt, 
in P(u], 44, -..,4n) das Vorzeichen von u; umzukehren. Da außerdem für 
jeden Vorzeichenwechsel an Stelle eines Faktors exp H ein Faktor exp — H 
tritt, ist leicht einzusehen, daß die Wirkungsweise der Matrix V, durch 


(II, 119) V,= II espH+ (exp —H)c; 
j 
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beschrieben wird. Man erkennt aus der Definition der Spiegeloperatoren, 
daß folgende Beziehungen gelten müssen 


(IT, 120) Cc;=c?’=1;C;Ck = cxC;. 


ONSAGER benutzt an Stelle der etwas unhandlichen Gleichung (II, 119) 
folgenden Ausdruck für V| 


(IH, 12) V,= (2 sinh 2H)®2 expH*bmitb = c; 
j-1 
und H*=}$InctghH =artghex&p— 2H 
Man kann Gleichung (II, 121) leicht mit Hilfe der Beziehungen der 


hyperbolischen Funktionen und der Reihendarstellung der Exponential- 
funktion unter Benutzung der Beziehung c;? = 1 verifizieren. 


In analoger Weise kann man nun die Wirkungsweise der MatrixV, 
durch Operatoren ersetzen. Gleichungen (II, 116) und (II, 118) zeigen, daß 


‚ hier nur die Multiplikationen des entsprechenden Terms mit einem Exponen- 


tialfaktor vorzunehmen sind. Definieren wir also einen Operator s;, der die 
Eigenschaft besitzt, aus dem Eigenvektor das j-Argument als Faktor heraus- 
zuziehen, so wird die Wirkungsweise von V, dargestellt durch 


(IL, 122) V; = exp H’a mit ae Sj Sj ]. 


Aus der Definition der Operatoren S; erkennt man sofort die Zusammen- 


hänge 

Sıs Ss; — S; Sk = Sk Sj- 
Damit erhalten wir für Gleichung (II, 116) folgende Beziehung (in Vektor- 
form geschrieben) 


Bern un... m )=(V.p)=(V; V,,p)= (2sinh 2H)j%2 (exp H’a 
x exp H*b, p). 


Zur Auswertung der Gleichung (II, 123) mit Hilfe der dort eingeführten 
Operatoren müssen wir ihre Eigenschaften untersuchen. Alle Operatoren 


Cı> ---, Ca; Sp -- > ---, Sn bilden offenbar eine endliche Gruppe (solange 
n endlich). Die Produkte der einzelnen Elemente sind: 
(I, 124) s?=c/?’=1;$s;$ = 5 $Sj; 06; ck = CK 0; SC; = — CjC;; 


S;Ck = Cxk S;j- 
Diese Beziehungen kann man auch unmittelbar aus den Definitionen ab- 
leiten. Man sieht daran sofort, daß alle Elemente (| =1,...,n)1,$,,C;,$SjcC; 
zusammen mit den nicht unabhängigen Elementen — 1, — Sj; — C,, C; $; 
eine Untergruppe bilden (Quaternion). Die Untergruppe besitzt eine ir- 
reduzible Darstellung, nämlich 
(II, 125) 


D()= Fe: D(s;) = ir h Ik D (c;) = I Disic (eat 


Der Leser kann sich leicht durch Ausmultiplizieren überzeugen, daß mit 
dieser Darstellung die Beziehungen (II, 124) für | =k erfüllt sind. Die Dar- 
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stellung der ganzen Gruppe kann daher mit Hilfe des direkten Produkts 
aus den n gleichen Darstellungen unserer Untergruppen erfolgen. Diese 
Darstellungen sind irreduzibel und der Ordnung 2%. Das gesamte Problem 
ist damit auf eine allgemeine Quaternionenalgebra im 2”-dimensionalen 
Vektorraum zurückgeführt worden. Es wäre nun die Aufgabe, alle Gesetze 
und Beziehungen dieser Algebra gruppentheoretisch abzuleiten und danach 
die Lösung des Eigenwertproblems (II, 123) auszuführen. Der mathematische 
Aufwand für diese gruppentheoretischen Gedanken ist leider nicht ganz 
einfach, so daß ich die nun folgenden rein mathematischen Überlegungen 
aus Raumgründen nicht darzustellen vermag. Da aber die hier verwendete 
Bezeichnungsweise, absichtlich weitgehend an die von ONSAGER angepaßt 


ist, kann der Leser die gruppentheoretischen Ableitungen der ÖNSAGER- 
schen Arbeit ohne besondere Einarbeit weiter lesen. Uns interessieren hier 


lediglich noch die Ergebnisse, die von ONSAGER ausgerechnet wurden. Die 
Größe T. ist durch die Gleichung 
sinh V/KT. sinh V’/kT. = 1 festgelegt. 

Für V = V’ wurde diese Beziehung bereits von KRAMERS uns WANNIER 
gefunden. Der Verlauf der spezifischen Wärme als Funktion der Temperatur 
zeigt für den unendlich großen Kristall zwar eine Unstetigkeitsstelle bei T,, 
aber rechts und links davon ist der Verlauf völlig stetig. Die spezifische 
Wärme wird an dieser Stelle mit In (T — T.) unendlich. Die Energie bleibt 
in jedem Falle eine stetige Funktion der Temperatur. In Fig. 12 ist der von 


Fig. 12. Ergebnis der Onsager’schen Berechnung der spezifischen Wärmen als Funk- 
tion der Temperatur für verschiedene Werte von 

mens ee Er 

vv v210 V 


UÜNSAGER gefundene Verlauf der spezifischen Wärme für verschiedene Ver- 


hältnisse von V’/V wiedergegeben. Das Verhalten interessiert besonders 
bezüglich des Bereichs der kritischen Temperatur. Man erkennt deutlich, 


daß die Fläche unter dem Maximum mit wachsender Anisotropie stark 


absinkt und damit experimentell kaum noch nachweisbar wird und schließ- 
lich ganz in den eindimensionalen Fall übergeht. 
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III. Fehlordnung und Röntgenbeugungsproblem. 


l. Allgemeine Gleichung für dreidimensionale Fehlordnung. 


Das Problem der Röntgenstrahlbeugung folgt an sich der gleichen 
Linie wie die vorausgegangenen Berechnungen, die von thermodynamischer 
Seite her durchgeführt wurden. Wenn wir, wie bei allen bisher behandelten 
Problemen, als Grundlage annehmen, daß die Geometrie des Gitters nicht 
gestört wird, lassen sich die Verhältnisse noch einigermaßen gut übersehen, 
trotzdem führt die praktische Berechnung noch eine nicht unerhebliche 
Komplikation formal mathematischer Natur ein, die aber nicht unlösbar ist. 
Bekanntlich verfährt man ja bei der Berechnung der Röntgenintensität 
eines exakt geordneten Kristalles so, daß man die Streuamplituden der 
Zellen, deren Lage durch Angabe der Translationsvektoren a,, a, 4,3 
gegeben ist, mit ihren Phasen summiert. Die Streuamplitude A ergibt sich 
dann gemäß der Formel 


h, l, ls 
cn j i 8 — 8 
(IL,1) A= >» > > F;;, exp 2 ni (aa + Jaadg + ]303, Zn > 


05-050 

dabei sind 1,, 1,, 1, die Anzahl der sich an die Zelle mit dem Index 0,0, 0 an- 
schließenden Zellen in Richtung der drei Haupttranslationsrichtungen. 
s,, s sind die Einheitsvektoren des einfallenden und gebeugten Röntgen- 
strahls, 4 die Wellenlänge der verwendeten Strahlung. 

Man kann in Gleichung (III, 1) auch die Vektoren des Fourier-Raumes 
einführen, die durch die Orthogonalitätsbedingungen 

0 für v # u 


L) 


(Gr, Bu — 7 für » N 


gegeben sind. Für einen Vektor des Fourierraumes gilt die bekannte Be- 


ziehung 
8 — 89 


hbı+kb, +15; =, 


wobei wir hier die h, k, | ganz allgemein als Vektorenkomponenten definieren 
wollen, die also nicht ganzzahlig zu sein brauchen. Damit wird 


[£ 


- f. , r B—18 nn E k 
exp2 riljıaı 7 JaAdg + J3Q3, —) = exp2riljyh+jek+ jpl). 


‘Wenn wir zur Intensitätsgleichung übergehen, erhalten wir 


(111, 2) 
E 5 FOR. ;. ee N 
J=_ Sa pri re) Eis he). 

Ivjajs J1’l2’,a”, 

Ist F, ;,;, für alle Zellen gleich, so ergibt die Summierung die bekannte 
Intensitätsgleichung für Kristalle. Wir wollen jedoch wissen, wie sich die 
Verhältnisse ändern, wenn F;;,;, statistischen Schwankungen unterliegt, 
in diesem Fall können wir F\ ,,,, nicht vor die Summe ziehen, sondern 


müssen anders verfahren. 
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Wir nehmen deshalb zunächst eine Umformung in der Summierung 
von Gleichung (ILL, 2) vor. Wir setzen 

Y=jhtm W=Shtm, "=, tm 
und erhalten 
(III, 3) 

2 22) E33, Kine exp — 2 ni (m;ıh+m,k-+ m;]). 

Jıjajs mı m, m; 

Wir haben diese Umformung deswegen vollzogen, weil in dieser Form 
unmittelbar das Prinzipielle des Problems zutage tritt. Offenbar hängt bei 
einer statistischen Abhängigkeit der Zellen untereinander die Besetzung 
irgendeiner Zelle davon ab, wie die benachbarte Zelle besetzt war. Das hier 
gestellte Problem sieht also so aus: Gegeben ist die Besetzung irgendeiner 
Zelle, gefragt ist nach der Wahrscheinlichkeit der Besetzung, der um m, 
Schritte in der a,-, m, in der a,- und m, in der a,-Richtung entfernten Zelle. 
Da ja in der Schreibweise der Gleichung (III, 3) die Zelle mit den Indizes 
Jı Ja, Js die Ausgangszelle und die mit den Indizesj, + m,, ja + m,, ja + my 
die Bezugszelle ist, und nach der Besetzungswahrscheinlichkeit gefragt wird, 
kann man offenbar so verfahren, daß man die Mittelwerte 


Er} je» Js Fi m, je+m;, js+m, 
aufsucht, und diese als unabhängig von der Ausgangszelle j,, ja, j; betrachtet. 
Damit läßt sich die Summation über j},, ja, ]; ausführen und man erhält für 
die Intensität 


(ILL, 4) 
hal +11 et max — 
Jen 3 2 (h— 1m) —m;|) (l3— m; )) FF, m,m, 


m, (h—1) m, (a&—1) m; (1.—1) 


xexp—2ri(mn+m,k+m,]). 
Das Hauptproblem ist also, die statistische Abhängigkeit des Mittel- 
wertes der Streuamplitude in Abhängigkeit vom Abstand der beiden Zellen 
zu ermitteln. Erfolgt insbesondere die Besetzung der Zellen statistisch un- 
abhängig, dann wird F F,* m, m, unabhängig von m,, m,, m, mit Ausnahme 
an der Stellem, =m, = m, =(0. Der Mittelwert an der Stelle m, =m, 
— ma — 0 ist offenbar gleich dem Mittelwert der Quadrate aller Möglich- 
keiten der Strukturamplituden; wir schreiben dafür F?. Der Mittelwert für 
m}, m,, m; ist aber gleich dem Quadrat der Mittelwerte. (— F?) Ziehen wir nun 
das Glied mitm, =m, =m; —=0 (=], 1,1, F?) aus der Summe (III, 4) heraus 
und bringen dafür ein Glied 1, 1,1, F? in die Summe hinein, so erhalten wir 
die bekannte Intensitätsformel für den idealen Mischkristall 


sin“ m-l,h sin m Lk sin nzler 


A,5) J=lll; [IP — F2] + F? 


> 


sin zn h” sinn kenne 
d. h. scharfe Interferenzen mit einem diffusen Untergrund, der wie die 
Differenz F2 — F? Intensitätsschwankungen in Abhängigkeit der h,k,l 
unterworfen ist. Die scharfen Interferenzen verschwinden, wenn F? —0. 
Das ist bei einer statistischen Lagenbesetzung gewisser Atomarten immer 
möglich (Lagenfehlordnung). Da das Glied F? unabhängig von der Phase, 
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die durch statistische Besetzungsmöglichkeiten der Zelle hervorgerufen wird, 
ist, F? dagegen nicht, kann F? größeren Schwankungen unterliegen als F?. 
Es ist unmittelbar klar, daß man mit vernünftigem Ansatz für die Abhängig- 
keiten, die man als richtig annimmt, Aussagen über das Verhalten der 
Röntgenintensitäten machen kann, und durch den Vergleich mit den Rönt- 
genaufnahmen feststellt, ob das angenommene Gesetz eine gute Überein- 
stimmung gibt oder nicht. Das ist aber im Prinzip noch unbefriedigend. Die 
zweifellos schönste Lösung wäre die, daß man die vorhandenen statistischen 
Ansätze versucht, in die Theorie der Röntgenbeugung einzubauen; und 
röntgenographisch untersucht, ob die gefundenen Beziehungen mit den ex- 
perimentellen Erfahrungen im Einklang stehen oder nicht. Für den drei- 
dimensionalen Kristall ist das Problem noch nicht befriedigend gelöst. 
Lediglich für den Fall des Systems Au—Cu wurde von A. J. C. WıLson 
ein Versuch unternommen, das Verhalten der Röntgenreflexe mittels einiger 
Modelle für die Übertrukturbildung im Bereich Au Cu, zu erklären. Die 
Berechnungen legen allerdings nicht die statistischen Ergebnisse der im 
theoretischen Teil besprochenen Arbeiten zugrunde. Ein solches Verfahren 
hätte schon deswegen keinen Zweck, weil die Theorien, die sich mit dem 
dreidimensionalen Modell befaßt haben, viel zu grob sind, um ein eini- 
germaßen befriedigendes Bild für das Verhalten der Röntgeninterferenzen 
zu geben. Hier hat lediglich die gruppentheoretische Methode Aussicht 
auf Erfolg. Diese ist aber für den dreidimensionalen Fall nicht ent- 
wickelt und für den zweidimensionalen Fall fehlen geeignete Substanzen, 
an denen der Vergleich von Theorie und Experiment nachgeprüft werden 
könnte. Die bei WrLson ausgeführten Ansätze sind, wie schon gesagt, nicht 
exakt und dürften bestenfalls gute Näherungen zur Erklärung des Nah- 
ordnungsverhaltens der näheren Umgebung eines herausgegriffenen Atoms 
sein. Allerdings kann eine qualitative Erklärung des Röntgenbeugungs- 
bildes durch Kombination zweier Modelle herbeigeführt werden. Die Haupt- 
annahme ist dabei, daß Au—Au-Nachbarn vermieden werden. Mehr als 
die Deutung des qualitativen Verhaltens kann aber von der Näherung nicht 
erwartet werden. 


III, 2. Zweidimensionale Fehlordnung. 


Liegt zweidimensionale Fehlordnung vor, ein Fall, der z. B. bei den 
Kettenstrukturen auftreten kann, dann ergibt sich in der Summierung 
der geordneten Translationsrichtung in Gleichung (III, 4) insofern eine 
Vereinfachung, als diese Summation statistisch unabhängig erfolgen kann. 
Legen wir diese als die a,-Richtung fest, so erhalten wir 


2 

(I.6) J=e— ee] Ba: l, — |mı|) (a — |m;|) F F,m, exp 
sn al u; 

—2rz im, h+mk). 

Bei statistischer Unabhängigkeit der Zellenbesetzung folgt analog zu Glei- 


chung (III, 5) 
in? l Ti, _- snzlh sin z ı,k 
(TIL, 7) —— sin’ all | lı l, ( 9 F?) | FR? IT 1 IMZeTt 9 


sin nz h sn? z k 
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D. h. wir erhalten neben diffusen Gitternetzebenen im reziproken Gitter 
normale scharfe Interferenzen. Verschwindet F? für bestimmte Werte von 
l, so sind die Gitterebenen für alle Werte h, k diffus. Wird dagegen 
FR—F, dann gibt es nur scharfe Interferenzen. Liegt eine reine Lagen- 
fehlordnung in dem Sinne vor, daß die Ketten zueinander zwei gleichwertige 
Lagen einnehmen können, sagen wir (000) und (003), dann wird 


| 1+expril |]? 
D) 4 


r—|E2 
B=|rje 


und wir erhalten 
PR—M=0 für 1= gerade, 
F=0 für 1= ungerade, 


d. h. die Gitternetzebenen im reziproken Gitter mit 1 = gerade haben nur 
scharfe Punkte und diejenigen mit 1 = ungerade sind rein diffus. Dieses 
allgemeine Gesetz ändert sich nicht, wenn wir statistische Abhängigkeiten 
einführen. Die diffusen Gitternetzebenen lösen sich, je größer die statistische 
Nahordnung ist, mehr und mehr in Interferenzpunkte auf. Sind die sta- 
tistischen Abhängigkeiten in den beiden Translationsrichtungen aus Sym- 
metriegründen verschieden, dann werden die diffusen Interferenzpunkte in 
den entsprechenden Richtungen des reziproken Gitters ebenfalls verschieden 
stark diffus sein, also keine Kreissymmetrie aufweisen. Wir wollen jedoch 
diesen zweidimensionalen Fall, der ja gemäß den Arbeiten von ONSAGER (52) 
und AsHkın und LAMmB (50) auch mit statistischer Abhängigkeit quanti- 
tativ gelöst werden kann, nicht näher durchdiskutieren. Die exakte experi- 
mentelle Prüfung der dreidimensionalen und zweidimensionalen Fehl- 
ordnung ist natürlich insofern nicht ganz einfach, weil an sich die Messung 
der Intensität im gesamten reziproken Raum bzw. in einzelnen reziproken 
Gitterebenen von Punkt zu Punkt durchgeführt werden muß. Für diesen 
Zweck würden sich monochromatische WEISSENBERG-Aufnahmen gut eig- 
nen; leider ist bisher für den zweidimensionalen Fall noch kein brauchbares 
Material gefunden worden, mit dem exakte Messungen ausgeführt werden 
könnten. 


III, 3. Eindimensionale Fehlordnung, Berechnung eines 
Spezialfalles. 


Viel günstiger liegen sowohl experimentelle als auch theoretische 
Ergebnisse bei der Behandlung der eindimensionalen Fehlordnung. Hier 
verfügen wir durch Arbeiten von WıLson (14), HENDRICKS und TELLER 
(54) und schließlich in allgemeiner Form durch Arbeiten des Verfassers 
(15) über ein ausreichendes Bild zumindestens für solche Kristalle, in 
denen thermodynamisches Gleichgewicht herrscht. Die letzte Frage ist 
nun vor allen Dingen auch für natürliche Minerale nicht leicht zu beant- 
worten und kann zu schwerwiegenden Konsequenzen in bezug auf das 
Röntgenbeugungsproblem führen. Alle Berechnungen des eindimensionalen 
Modells haben ergeben, daß für diesen Fall, streng genommen, keine Fern- 
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ordnung auftreten kann. In Wirklichkeit tritt sie natürlich dann auf, wenn 
die ungünstigen BoLtzmann- Faktoren eines Fehlers so groß werden, daß der 


! Eintritt eines Fehlers ein unwahrscheinliches Ereignis ist; wenn also ent- 


weder V sehr groß oder T sehr klein ist. Man muß nur darauf achten, daß 
die Temperaturen für die Einnahme des thermodynamischen Gleichgewichts 
nicht zu gering sind. (Aktivierungsenergie für einen Platzwechsel zu hoch im 
Vergleich zur Temperatur.) Die Berechnung komplizierter Fälle des ein- 
dimensionalen Modells erfolgt am besten mit Hilfe von Differenzenglei- 
chungen. Wir wollen dieses Verfahren erläutern. 


Für das eindimensionale Modell vereinfacht sich Gleichung (III, 4) zu 
1-1 
B=R?>(,—|m,|) FFi exp —2rim,| 
Net 


4 sin® zl,h sin zl,k 
mit == 


sin zz h,k sinsck 


' h, k haben also nur für ganzzahlige Werte einen Sinn. 


Das Kernproblem liegt darin, ausgehend von einem beliebigen Glied 
der Kette den Mittelwert des Ausfalls eines um m,-Schritte entfernten 
Gliedes zu ermitteln. Die Lösung dieses Problems ist eine reine Frage der 
Wahrscheinlichkeitsrechnung und ist im Fall von Wechselwirkungen nur 
zwischen nächsten Nachbarn sehr elementar zu lösen. Zur Erläuterung des 
Verfahrens wählen wir den allgemeinen Fall von A- und B-Gliedern unter 
Berücksichtigung der Wechselwirkungen nicht nur der nächsten, sondern 
auch der übernächsten Nachbarn. Wir bauen also unsere Kette auf, indem 
wir Glied für Glied anfügen und die Wahrscheinlichkeiten aus den Wechsel- 
wirkungsenergien berechnen. Da wir auch diejenigen der übernächsten 
Nachbarn mit in die Rechnung einbeziehen wollen, definieren wir neue 
Wechselwirkungsenergien 


Vans, Vase, Vara, Van, Vrar, Vein, VoBA, VeBB; 
wobei die Indizes die zugehörige Dreiergruppe angeben. Wir nehmen dabei 
bewußt davon Abstand, die Wechselwirkungsenergien so zu definieren, daß 
wir Vaaz zerlegen in Vaa + Van + V’an, wobei V’,z die Wechselwir- 


kungsenergie vom 1. zum 3. Glied darstellen soll. Denn es ist in vielen Fällen 
sicher so, daß die Wechselwirkungsenergie der übernächsten Nachbarn nicht 
unabhängig von der Konfiguration der nächsten Nachbarn sein wird. 


Wir wollen nun aus Gründen der Einfachheit der Durchrechnung das 
obige Problem symmetrisch machen; diese Vereinfachung ist aber keine 
grundsätzliche Schwierigkeit, man kann die Lösung auch ganz allgemein 
durchführen. Unsere Symmetriebedingungen sollen lauten: 

1. Vaaa = VoBB; 
2. Varp = Vona, 
3. Vara = Vrar. 
D.h. die Atomarten verhalten sich symmetrisch zueinander, diese Bedingung 
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braucht nicht immer erfüllt zu sein. Dagegen ist Vaap — Vraa für unseren 
Fall von A- und B-Teilchen immer erfüllt. Wir schreiben also | 


ei Var _ __ VBBB _H 
MI Kitas 
VAAEN EV BEA Via 2 VıBB an 
ET ET kT ; 
VpaR un Vara Sf 
kT kT 


Aus diesen Werten lassen sich die Wahrscheinlichkeiten W, yz, daß auf ein | 
x-y-Paar ein z-Atom folgt, bestimmen. Es ist | 


(ILL, 8) Wiss Wers 1 un = as — 7%, 
Warn We Te ee Hy en u: 
Wann = Wrau = + en H, Bi a ar 
Wir = Wen = — Hr re ß. 


l1-+ expH, cosh H2/2 


Diese Beziehungen, die also nunmehr absolute und keine relativen 
Wahrscheinlichkeiten mehr sind, ermöglichen die einfache Berechnung des 
eindimensionalen Falles. Wir bezeichnen mit P,„ die absoluten Wahrschein- 
lichkeiten im Abstand m ein gleiches Atom zu finden, und mit 1 — P„ die 
Wahrscheinlichkeit im Abstand m ein ungleiches Atom zu finden, dann läßt 
sich offenbar folgende Rekursionsformel für die P„ aufstellen 
Pa=Pn-2Pnm-ı® + Pn-: (l ei —rn)) Pn—ı (ls 

+ (l ——n23) (l — Ense) (1 — a) 
oder 
UL) Pr Berner ee 
Die homogene Gleichung zu (LIL, 9) lautet 
(III, 10) 
Dan ap a Ele 20 Era an p—ı) = 
Das ergibt die charakteristische Gleichung 
dt) Brad ae year 


deren Lösungen 


| 


Lo razı 


. 
= ——-— 4 Va — PP? +4(a + B—1) sind. 
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Die allgemeine Lösung der homogenen Gleichung (III, 10) ist also 


(ILL, 13) Pm =A + Bx,n + Cx,m, 
Einsetzen von (III, 13) in die inhomogene Gleichung (ILL, 9) ergibt 
A=}, 


B und C lassen sich aus den Randbedingungen 


P,=1, P,=1 (jenachdem ob ein gleiches oder ungleiches Paar 
P)=0 Ausgangspaar ist) 
berechnen. 


Man erhält so für die Konstanten B und C 


&) Ausgangspaar war AA oder BB (B,, C,) 


Be 2 orrP 
) un 
E19 eV ale Mae 
,=4— 
8 4 4 V ’ 
b) Ausgangspaar war AB bzw. BA (B, + md C,) 
| ir 
Ur: ; — 
(ILL, 15) N 
a: 2+0— Pf 
Ey 


ı Damit ist für diesen Fall das Problem der Auffindung des für die Berechnung 
des röntgenographischen Beugungsproblems wichtigen Mittelwertes FF}, als 
gelöst zu betrachten. 


Ich will an dieser Stelle noch bemerken, daß wir in gleicher Weise auch 
die Zustandssumme des Systems finden, wenn wir nicht, wie wir es hier 
| getan haben, die absoluten, sondern die relativen Wahrscheinlichkeiten, die 
durch die BoLrzmann-Faktoren gegeben sind, benutzen. Man muß dann, da 
die Beziehung P) = 1—P,, für gleiche und ungleiche Atome nicht gilt, 
zwei Differenzengleichungen für gleiche und ungleiche Atome aufstellen. 
‘ Diese lassen sich aber wieder auf eine zurückführen, und das Ergebnis der 
beiden Lösungen addiert, ergibt die Zustandssumme, aus der alle thermo- 
dynamischen Eigenschaften abzuleiten sind. Da das Verfahren ganz analog 
zu der Durchrechnung des eindimensionalen Falles in Kap. II verläuft, 
soll hier nicht darauf eingegangen werden. Der Leser kann an Hand des 
vorstehenden Gedankengangs den Lösungsweg leicht selbst ausführen. 


Wir kommen nun zu unserem Beugungsproblem zurück. Zur Aus- 
wertung der Intensitätsgleichung (11, 5) müssen wir zunächst FF aus- 
rechnen. Da wir zwischen gleichem und ungleichem Ausgangspaar unter- 
schieden hatten, muß zunächst ermittelt werden, mit welcher Wahrschein- 
lichkeit wir beim Durchlaufen des Kristalles ein gleiches oder ungleiches 
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Paar finden werden. Nennen wir die Wahrscheinlichkeiten W, und W,. 
Gemäß (III, 8) sind die Wahrscheinlichkeiten 

1. daß ein gleiches Paar auf ein gleiches folgt a, 

2. daß ein ungleiches Paar auf ein gleiches folgt 1 — x, 

3. daß ein gleiches Paar auf ein ungleiches folgt 1 — ß, 

4. daß ein ungleiches Paar auf ein ungleiches folgt P£. 


Es bestehen die Beziehungen 


W, Il — 8 
nl = 
ai, W\ I 
daraus folgt | 
1—B 1—ı& 
IIL, 16 We «Wa = 
en en) | 
Damit lassen sich die FF#, berechnen. Es ist 
(BETT) 
ng Bo Ike an 
FF = W; u Pas Fi Ar (1 En a) F5 Sr , PragF& RE (Pre | 
Ban ee 5 
+ W, =” N ir (1 — in FE as 5 ern FF Zw (1—Pau) FX n 


F, = Atomfaktor der A-Atome. 
FR = Atomfaktor der B-Atome. 
(A + B-Atome hier im gleichen Mengenverhältnis vorausgesetzt.) 


Einsetzen von (ILI, 13) in (ILI, 17) unter Berücksichtigung, daß zu dem : 
Wert Ps die Konstanten B,, C, und Pu, Bu, ©, gehören, ergibt schließlich : 
für ER", 

(LET, 18), 0% = , I mit! |R miese Ren 
mit K, = W;B, + WıuBu, 
Re ww EMS: 


Setzt man (III, 16) in (ILL, 5) ein, so erhält man schließlich die Inten- : 
sität der gebeugten Röntgenstrahlung 


in2 sel, 1 
(EIT, 19) vn 13 a 
sin? zl 
l; I 1 — x, 122 be; 
+2, —-FBR|LIK i 2 ER S | 
> | x »| \ 21 23,008 12x22 wi a | 


wobei 1, groß angenommen wurde und dadurch einige Glieder bei der Aus-' 
führung der Summation vernachlässigt werden konnten. Wir betrachten! 
nun die vier Grenzfälle der Ordnung 


Il Ba | Du | 
| = 
Bel Bl! 
4.6 = B=0 
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Fall l: 
8=-l,=—-1;) RB=4; R=4ı 
%—=1; $=1 bedeutet aber auf AA folgt A 
auf BB folgt B 
auf AB folgt A 
auf BA folgt B 
d. h. es existieren zwei Ordnungszustände nebeneinander AAAA ........- 
und der äquivalente BBBB ......... Bowie BABeEree 


Gleichung (III, 19) ergibt scharfe Maxima für 1 = ganzzahlig, und 
1 = halbzahlig im Einklang mit der Erwartung der beiden Strukturmöglich- 
keiten. 


Fall 2 
a VERK, IR, 0 
&«—=1;  =0 bedeutet aber auf AA folgt A 
auf BB folst B 
auf AB folgt B 
auf BA folgt A. 
Es existiert also nur ein Ordnungszustand AAA ........ oder der äqui- 
Bee BBB .......... 


Gleichung (III, 19) ergibt demgemäß ein Maximum an der Stelle 
l = ganzzahlig, wie es auch zu erwarten ist. 


wer, 0, = 17 R,- 0R, =} 

@—=0; = 1 bedeutet aber auf AA folgt B 
auf AB folgt A 
auf BB folgt A 
auf BA folst B. 

Es existiert nur der Ordnungszustand ABAB ........ 


Es gibt also Maxima für 1 = ganzzahlig und | = halbzahlig im Fin- 
‚klang mit der Erfahrung. Das Intensitätsverhältnis ist natürlich im Ver- 


| gleich zu Fall (1) verändert. 


| 
Ball2ı ,=-,,=—-1;K, =: R,=+ 


| o=N0; Bel Dean aber auf AA folgt B 

| auf AB folgt B 

auf BB folgt A 

auf BB folgt A. 
Es existiert nur der Ordnungszustand AABBAABB ........ 


Die Maxima liegen an den Stellen 

l = ganzzahlig 

l=+ (mod l) 

l=3 (mod ]). 
Auch das steht im Einklang mit der Erfahrung, weil für die geordnete 
Struktur AABB......... die Strukturamplitude 
S=F,(l1+exp2ril)tesp4rilFfs, (l+exp2nil) 
sein würde. 
Ile 
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Man sieht unmittelbar ein, daß 


®=4|F,+Fs|? für 1 = ganzzahlig 
SE) für 1 = halbzahlig 
S2=2|F,—Fy|? für 1 = # (mod }). 


Wir wollen nun noch die beiden Sonderfälle «a=ß und «a=1-—$ 
betrachten, aus denen einiges über den Verlauf der Röntgenintensitäten 
zu entnehmen ist. 


Für «= £ vereinfachen sich die Gleichungen (III, 12), (ILL, 15) und 
IT I6)Zu 


Die Intensitätsgleichung (III, 19) wird damit 


ET |FıA+FR|? sin? zz 1,1 
2) 45 | 4 sin? sul 
1.1; 1—« 1—« | | 
nr, [22 | | 
+ IF 2 | 8la— Y2a—lcos2rl rn 
oder 
F,+Fr|?’smnll 1 ’ a(l—.a) 

ln a EN 4 

| 4 sineal ' | x = © — (2% —1)co®2n1 


Das erste Glied ergibt unabhängig von h, k scharfe Interferenzen an den ı 
Stellen 1 = ganzzahlig unabhängig von «. | 


Das zweite Glied hat für 2«— 1<0, a<} 


Maxima für 1 = 4 (mod #) proportional zu 


Minima für 1 = 0 (mod #) proportional zu 1 _ 
—ıa 


Für 2«—1> 0 sind die Lagen von Maxima und Minima vertauscht, | 
wie man leicht einsieht. Der Leser wird sich die Gründe für das Verhalten! 
der Röntgenintensität leicht selbst erklären können, wenn man beachtet, 
welche Grenzwerte der Ordnung mit dem Übergang « =0&— 1 verknüpft! 
sind. 
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Für « = 1 — ß ergeben sich folgende Werte für 
Selig =2Rr —- 1; %-=V0); 


3—2x 3—2x 1+ 2x 
B, = 1 — il —.4 
a eo, ia 4 da ıy m rat, 
1+2« 
G=14- 4201. 
Wax W le; 
R,=3; KR = 


Einsetzen in Gleichung (III, 19) ergibt 
|Fı+Fs] sin? se 1,1 
4 sin? ul 


J-R| 


«(1 —.«) 
1—2(2&x —1)cos2n1+(2«— 1) 


Für 2« — 1 > 0 liegen die Maxima des diffusen Gliedes von Gleichung 
(ILL, 19) an den Stellen 1—=0 (mod 1); für 24 — 1 <0 dagegen an den 
Stellen 1— 4 (mod 1); für die scharfen Reflexe hat sich nichts geändert. 

Dieses Ergebnis kann man sich an Hand der Grenzwerte a =0 > a —1 
in analoger Weise deuten. Die durchgerechneten Grenzfälle geben nun nicht 
alle Eigenschaften des hier gewählten Modells wieder, und zwar deshalb, 
‚weil neben diesem Diffus- und Schärferwerden der Maxima noch Ver- 
schiebungen in ihren Lagen auftreten können. Diese Frage wurde vom 
‚Autor (15) aber bereits ausführlich am Beispiel der dichtesten Kugelpak- 
kungen diskutiert. 

Es besteht keine prinzipielle Schwierigkeit, das eindimensionale Modell 
der Fehlordnung auf beliebige Fälle zu erweitern, allerdings bringt die Ein- 
führung weiterer Wechselwirkungen die Komplikation mit sich, daß die 
dabei auftretenden charakteristischen Gleichungen zu unhandlich werden, 
‚um alle Eigenschaften des entworfenen Modells in einfacher Weise über- 
‘sehen zu können. 


+1 |F4 — Fol? 


III, 4. Vergleich mit experimentellen Ergebnissen. 

Die experimentelle Prüfung der hier entworfenen Vorstellung stößt 
aber auf große Schwierigkeiten; dies wurde früher vom Autor auf die 
natürlich strittige Frage des thermodynamischen Gleichgewichts der unter- 
suchten Kristalle zurückgeführt; alle bisher abgeleiteten Formeln gelten 
in Strenge ja nur dann, wenn sich der Kristall im thermodynamischen 
Gleichgewicht befindet. Da wir nun darauf angewiesen sind, das eindimen- 
sionale Modell röntgenographisch am dreidimensionalen Kristall zu prüfen, 
muß vorerst die Frage erörtert werden, inwieweit sich die vom Autor 
durchgeführten Berechnungen auf die eindimensional fehlgeordneten Kri- 
stalle übertragen lassen. 

Für das Eintreten normaler Unordnung sind die Platzwechsel der 
Atome im Kristallgitter die Ursache. Im dreidimensionalen Kristall ist aber 
für das Auftreten eindimensionaler Fehlordnung eine ganze Reihe von 
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Platzwechseln notwendig, nämlich die Umordnung einer gesamten Netz- 
ebene. Die Tatsachen der experimentellen Beobachtung scharfer Gitter- 
stäbe beweist aber, daß ein einmal an einer Stelle einer Netzebene aufgetre- 
tener Fehler sich sehr rasch über die ganze Netzebene ausbreiten muß, da 
sonst keine Schärfe der Interferenzen in zwei Richtungen des reziproken 
Gitters auftreten sollte. Da aber das thermodynamische Gleichgewicht eines 
Kristalls unabhängig vom Weg ist, auf dem es erreicht wird, sollte also das 
bisher ausgeführte Berechnungsverfahren auch für diese Fälle zulässig sein. 

Eine zweite Schwierigkeit erwächst aus der Tatsache, daß bei den hier 
zugrunde gelegten Berechnungen eventuelle Gitterspannungen durch Ände- 
rungen der Gitterkonstanten nicht berücksichtigt wurden. Bei einem in 
dieser Richtung verbesserten Ansatz würde es notwendig, eine zusätzliche 
Arbeit pro Flächeneinheit der Oberfläche der Gebiete verschiedenen Ord- 
nungszustandes einzuführen. Das ist zwar nicht ganz exakt, da an sich die 
eventuell vorhandene Verzerrungsarbeit nicht streng von der Oberfläche, 
sondern auch vom Volumen eines ungeordneten Gebietes abhängt. Die 
statistischen Zusammenhänge der reinen Wahrscheinlichkeiten werden da- 
durch in dem Sinne geändert, daß ein eingetretener Fehler einen gleichen 
„Fehler‘‘ in den benachbarten Schichten wahrscheinlicher macht, als es der 
normalen Fehlerwahrscheinlichkeit entspricht. Dies führt in der statistischen | 
Verteilung der durchschnittlichen ‚Dicken‘ der umgewandelten Gebiete 
zur Bevorzugung besonders der Gebiete größerer Dicke, so daß also in solchen 
Fällen die Berechnungen nickt mehr streng gültig sind. Man kann sich 
aber leicht überlegen, in welcher Richtung sich die Intensitäten der Röntgen- | 
interferenzen beim Vorhandensein solcher Spannungen ändern müssen. 

Für die Berechnung der Fehlerwahrscheinlichkeiten brauchen aber nur 
Verzerrungen innerhalb der Netzebenen berücksichtigt zu werden, da ja 
eventuelle Abstandsveränderungen senkrecht dazu bereits in den Wechsel- . 
wirkungsenergien enthalten sind. 

Wir wollen für die Betrachtungen der Röntgeninterferenzen in verzerrten | 
Gittern die a,-, a,-Achse des Kristalls in die geordnete Netzebene legen; ; 
dann steht der reziproke Vektor b, senkrecht auf a,, a,. 

Gitterkonstantenänderungen in der Unordnungsrichtung (b,) bringen ı 
immer Änderungen in den Phasengliedern unserer Summen in Gleichung : 
(III, 4) mit weil nunmehr nicht mehr der strenge Abstand der Netz- 
ebene eingesetzt werden darf. Diese Tatsache hat zur Folge, daß im rezi- 
proken Gitter kein einheitlicher Identitätsbereich für den Intensitätsverlauf | 
auf den Gitterstäben angegeben werden. Bei eindimensionaler Lagenfehl- : 
ordnung (Fehlordnung mit gleichen Netzebenen, die verschiedene Lagen | 
einnehmen können, z. B. dichteste Kugelpackung) kann man das besonders : 
gut auf dem (O0l)-Gitterstab erkennen, für den ohne Gitterspannungen alle 
Interfer enzen scharf sein müssen. Wenn eine Aufspaltung der Maxima oder ' 
ein Diffuswerden auf diesen Gitterstäben vor allen Dingen bei Interferenzen | 
mit hohen I-Werten auftritt, dann sind solche Störungen da, und man muß | 
in einem solchen Fall sehr vorsichtig bezüglich der Anwendung der bisher ' 
gegebenen Formeln sein. | 

Spannungen innerhalb der Netzebenen sind zwar für den streng geord- | 
neten Kristall mit eindimensionaler Fehlordnung kaum denkbar, da bei! 
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unendlicher Ausdehnung der Netzebenen sie zu empfindlichen Störungen 
im Gitterbau führen müßten. Aber immerhin könnten sie bei nicht beliebig 
in den a,-, a,-Richtungen ausgedehnten Kristallen auftreten. Das bedingt 
natürlich eine gewisse Unschärfe der Gitterstäbe bei großen Glanzwinkeln, 
verbunden mit einem Intensitätsabfall. 

Wir verfügen jedoch über einige experimentelle Fälle, bei denen Gitter- 
verzerrungen keine oder nur eine untergeordnete Rolle spielen. Der eine 
Fall ist das Co, das bei 400° ein Umwandlungsgebiet vom kubischen ins 
hexagonale Co besitzt und von EDWARDS und Lipson (55) experimentell 
und Wırson (14) theoretisch untersucht wurde. Leider wurden nur Pulver- 
aufnahmen ausgewertet, so daß die Untersuchungen nicht ganz als gesichert 
angesehen werden dürfen. Aber immerhin wurden mit einer geringen Schwan- 
kung aus verschiedenen Interferenzen ein einheitlicher Fehlordnungsgrad 
berechnet. Wir haben es hier zwar mit einem Umwandlungsvorgang zu tun, 
der an sich nicht streng in das besprochene Gebiet fällt. Solche Umwand- 
lungen können ja nur eintreten, wenn die Wechselwirkungsenergien tem- 
peraturabhängig sind, was bisher bei den Berechnungen nicht berücksichtigt 
wurde. Für die Fehlerwahrscheinlichkeiten macht das aber nichts aus, da 
ja in ihnen sowieso schon eine Variable, nämlich die Temperatur, enthalten 
ist. Der Zusammenhang zwischen den Fehlerwahrscheinlichkeiten ist dann 
eben nicht allein durch die Temperatur, sondern noch durch die mit T 
variablen Wechselwirkungsenergien gegeben. Ihre Konfigurationsabhängig- 
keit ist in dem Ansatz durch die Einführung verschiedener Wechselwirkungs- 
energien für die möglichen Anordnungen berücksichtigt. Es besteht also 
kein Grund, die Gültigkeit der berechneten Ansätze zu bezweifeln, wenn 
stärkere Gitterspannungen fehlen. Sie sind so allgemein gehalten, daß sie 
bei Einführung genügend weitreichender Wechselwirkungsenergien auch für 
homöopolare und sogar heteropolare Kristalle gelten sollten. 

Die im hiesigen Institut mit synthetischem Zn—S-Material von H. 
MÜLLER*) durchgeführten Untersuchungen, bei denen Temperungen beson- 
ders im Umwandlungsgebiet Zinkblende und Wurtzit bei sehr genauer Ein- 
haltung der Temperaturen von teilweise sehr langer Dauer durchgeführt 
werden, sowie eigene Untersuchungen der verschiedenen SiC-Strukturen*) 
haben ergeben, daß die gewonnenen experimentellen Daten mindestens zu 
einem gewissen Teil mit keiner nur auf dem Wechselwirkungsprinzip der 
Atome beruhenden Theorie zu deuten sind. In diesen Fällen treten sicher 
keine merklichen Spannungen auf, weil sowohl die Gitterstäbe bei großen 
Glanzwinkeln als auch die (OOl)-Interferenzen scharf bleiben. Um das ano- 
male Verhalten klar herauszustellen, seien hier die experimentellen Ergeb- 
nisse kurz beschrieben. 

Die von MÜLLER durchgeführten Untersuchungen am ZnS hoher Rein- 
heit bestätigen im Prinzip die schon bekannten Tatsachen, daß die Zink- 
blende beim Erhitzen bei 1120° in Wurtzit, und Wurtzit seinerseits dicht 
unterhalb 900° wieder in Zinkblende umgewandelt werden. Die Umwand- 
lungshysterese ist wohl auf geringfügige Änderungen der Gitterkonstanten 
zurückzuführen. Genaue Abtastung des Umwandlungsgebietes durch lange 


*) Erscheint demnächst. 
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Temperungen ergaben die Existenz eines sehr schmalen Umwandlungs- 
gebietes von ca. 10—20° Breite bei 900°, in dem plötzlich sehr viele, zum 
Teil völlig neuartige fehlgeordnete Kristallstrukturen des ZnS auftreten. 
Es besteht kein Zweifel darüber, daß bei diesen Versuchen das thermo- 
dynamische Gleichgewicht erreicht werden mußte. Denn bei wesentlich 
tieferen Temperaturen reichten die gleichen Temperzeiten aus, um den 
Wurtzit vollständig in Zinkblende überzuführen. Die neuen Strukturen 
haben nun die Eigenschaft, daß sie ähnlich hohe Perioden aufweisen, wie 
sie auch bei den SiC-Strukturen schon seit langer Zeit bekannt sind. | 

Durchaus analoge Verhältnisse zeigten die eingehenden Untersuchungen 
des Verfassers am SiC. Auch hier erfolgt die Umwandlung analog zum ZnS 
durch reine Lagenumordnung der Si—C-Schichten (Lagenfehlordnung). Hier 
wurden verschiedene ‚Strukturen‘ mit einer Identitätsperiode bis 600 A 
oft in Verbindung mit eindimensionaler Fehlordnung entdeckt. Wenn man 
nun versucht, diese Strukturen nach dem bisher verwendeten Wechsel- 
wirkungsprinzip aufzubauen, wie das vom Verfasser (15) getan wurde, dann 
kommt man hier zur Reichweite der Wechselwirkungen, die völlig sinnlos 
sind (bis zu 200 Ä); besonders beim SiC, wo wir wohl mit Sicherheit an- 
nehmen können, daß die chemische Bindung vorwiegend homöopolar ist, 
und deswegen nicht so hohe Reichweiten in ihren Wechselwirkungen haben 
sollte, kann man keine Begründung für die Ausbildung so langer Perioden 
finden. 

Es hat sich nun gezeigt, daß zur Aufklärung des strukturellen Ver- 
haltens der Aufnahmen fehlgeordnete SiC-Kristalle, fast für jede Aufnahme 
ein neuer mathematischer Ansatz eingeführt werden müßte. Ich habe mich 
deshalb darauf beschränkt, nachzuweisen, daß die Röntgenintensitäten unter 
vereinfachten Maßnahmen verhältnismäßig gut in einem Fall aufgeklärt 
werden können. Der dazu eingeführte Ansatz wird dem Leser zwar zunächst 
ziemlich sinnlos erscheinen, aber es hat sich durch Vergleich mit ca. 150 
weiteren Aufnahmen des SiC gezeigt, daß dieses merkwürdige Verhalten 
für alle Kristalle in gleicher Weise gilt. Zur Auswertung wurde ein Kristall 
herangezogen, der neben starken Reflexen der SiC-6-Schichtstruktur (hkk 
Bez. sh. (15)) auf diffusen Gitterstäben relativ gute, aber doch diffuse 
Reflexe der 15-Schichtstruktur (hkkhk) aufwies. Es wurde nun eine Sta- 
tistik in der Weise aufgebaut, daß Bauelemente hkkhkk und hkkhk ein- 
ander abwechseln, und zwar soll aufhkkhkk mit einer Wahrscheinlichkeit & 
ein gleiches bzw. 1 — « das andere Bauelement folgen. Die entsprechenden 
Wahrscheinlichkeiten für das Bauelement hkkhk wurden mit ß und (1— ß) 
definiert. Die Auswertung der Aufnahme ergab nun eindeutig, daß eine 
Klärung des Intensitätsverlaufs auf den diffusen Gitterstäben dann erzielt 
werden kann, wenn man die Fehlerwahrscheinlichkeiten so ansetzt, daß auf 
gleiche Bauelemente mit wesentlich höherer Wahrscheinlichkeit wiederum 
gleiche Bauelemente folgen, d. h. die Bauelemente müssen trotz ihrer Länge 
„wissen“, zu welcher Kategorie ihre Nachbarn gehören. Das gleiche gilt 
auch für die hochperiodischen Strukturen. Das Wesentliche an der Unter- 
suchung ist aber, daß ein solches Bauprinzip auf kompliziertere Systeme 
anwendbar ist. Man kann feststellen, daß sich verschiedene Teile desselben 
Kristalles strukturell vollkommen verschieden verhalten. Bei MÜLLERS 
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Untersuchungen am ZnS$ ist das besonders deutlich festzustellen gewesen. 
Das Bauprinzip ist also für eine bestimmte Temperatur nicht einheitlich, 
sondern sehr mannigfaltig. Aus diesen Gründen hat der Verfasser schon in 
(15) Bedenken über die kristallchemische Realität dieser Strukturen ge- 
äußert. Diese Bedenken haben sich nun durch die neueren Untersuchungen 
weitgehend bestätigt. 


Ich glaube, daß der Grund für das merkwürdige Verhalten darauf zu- 
rückzuführen ist, daß im Entropieglied der freien Energie bei allen Ab- 
leitungen nicht die Schwingungsentropie des Kristallgitters mit einbezogen 
wurde. Ähnliches wurde schon von Nıx und SHockL£y (11) für das System 
Cu-Au vermutet. Ihre Überlegungen sind allerdings auf diese Fälle nicht 
anwendbar, weil wir es hier einerseits mit einer reinen Störung der Gitter- 
geometrie zu tun haben, und außerdem in den bearbeiteten Fällen eine 
sehr starke Kopplung der Eigenschwingungen eines Atoms durch die homöo- 
polare Bindung vorliegt. Es wäre also notwendig, das Verteilungsgesetz der 
Eigenschwingungen des gesamten Gitters unter Berücksichtigung der Kopp- 
lung der Gitterschwingungen zu ermitteln. Das Problem ist bis heute noch 
nicht exakt für geordnete Kristalle gelöst worden und somit ist die Lösung 
für den fehlgeordneten Kristall ein ziemlich hoffnungsloser Fall. 


Wenn wir jedoch die Hypothese einführen, daß die Besetzungsdichte 
der Energiezustände der Eigenschwingungen stark von der Symmetrie einer 
Anordnung in dem Sinne abhängt, daß die Dichte mit wachsender Sym- 
metrie ansteigt, dann kann man eine qualitative Erklärung der vorliegenden 
Verhältnisse geben. Auch die Einführung der Identität in der fehlgeordneten 
Translationsrichtung bedingt eine höhere symmetrische Anordnung. Die 
Schwingungsentropie des Kristalls verläuft dann in dem Sinne, daß den 
Anordnungen mit hoher Dichte der Energieniveaus ein höherer Wert zu- 
kommt als Anordnungen mit einer geringeren Dichte. Den stabilen Zustand 
erhalten wir wie üblich über das Minimum der freien Energie. Es gilt also 
bei einer bestimmten Temperatur 


(III, 20) | F=E—-T (8, + $,), 
= _ dE En | e=0 
os T= const ds ds ds 


Die Gesamtentropie © wurde also in ihren Konfigurationsanteil ®, und ihren 
Schwingungsanteil ®, zerlegt. Wenn nun die Schwingungsentropie struktur- 
empfindlich ist, dann muß sie sich schon bei geringen Störungen des Gitters 
sehr bemerkbar machen, das würde bedeuten, daß die stärksten Abweichun- 
gen in der Nähe von $ — 1 auftreten. Unter dieser Voraussetzung wurden 
drei mögliche Kurven für die Schwingungsentropie in Fig. 13 und ein ent- 
sprechender Verlauf für die Konfigurationsentropie eingetragen. Man sieht 
daraus, daß die Gesamtentropie (ausgezogene Kurven in Fig. 13) sehr davon 
abhängt, wie groß der Unterschied der Schwingungsentropie im geordneten 
und ungeordneten Zustand ist. Ist ®, im geordneten Zustand wesentlich 
größer, so fällt die Gesamtentropie von S = 1 bis zu einem Minimum und 
steigt dann wieder an. Da in den Fällen ZnS und SiC die Änderung der 
Energie mit dem Ordnungszustand sicher sehr klein ist, können wir den 
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Energieanteil in Gleichung (III, 20) zunächst vernachlässigen. Dann be- 
deutet aber im Fall der abnehmenden Gesamtentropie, daß die freie Energie 
ansteigt, also der geordnete Zustand stabiler ist als der ungeordnete. 


p 


S=A 
Fig. 13. Verhalten der Gesamtentropie bei der Annahme, daß die Schwingungsentropie 
vom Ördnungszustand abhängt. 


Es existiert aber noch ein zweites Maximum, das an der Stelle 8 —0 
liegen muß. Die Höhe dieses zweiten Maximums hängt sehr davon ab, wie 
stark sich Schwingungsentropie und Konfigurationsentropie absolut unter- 
scheiden und wie groß die Differenz von &, und d, fürrS=1undS = 0 ist. 
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Die Vernachlässigung der konfigurationsabhängigen Energie hat nichts zu 
sagen. Wenn wir die Temperatur nur hoch genug wählen können, dann 
wird das Entropieglied in Gleichung (III, 20) en und die Stabilität 
des Zustandes bestimmen. 

Unter diesen Voraussetzungen kann man nun BE Verhalten der ZnS- 
und SiC-Strukturen verstehen. 

1. Diese Erscheinungen treten bei hohen Temperaturen auf, in einem 
Gebiet also, in dem das Entropieglied der freien Energie sehr bedeutend ist. 

2. Die Ausbildung gewisser Ordnungszustände bei hohen Temperaturen 
trotz gleicher oder annähernd gleicher Energie E der Strukturen ist auf die 
Abnahme der Gesamtentropie mit steigender Unordnung zurückzuführen. 

3. Die Frage, welche Strukturen mit größter Häufigkeit auftreten, ist 
durch die Symmetrie bedingt. Langperiodische sind also seltener als kurz- 
periodische, hochsymmetrische häufiger als unsymmetrische Strukturen. 
Damit würde sich also auch zwangsläufig erklären, warum gerade bei den 
hochperiodischen Strukturen des SiC die mit der höheren rhomboedrischen 
Symmetrie am häufigsten sind. 

4. Auch die Existenz vieler, in ihrer Struktur verschiedener Gebiete wird 
damit klar, da der Ordnungszustand durch Wahrscheinlichkeitsgesetze be- 
stimmt wird. Die Neigung zur Identität, die ja ebenfalls eine höhere Sym- 
metrie der Anordnung darstellt, ist dabei das einzige Bauprinzip. Führt man 
die in der vorangegangenen Überschlagsrechnung vernachlässigte Energie 
als Funktion des Ordnungszustandes wieder ein, so läßt sich einsehen, daß 
diese nur in dem Sinne wirken können, daß solche Ordnungszustände, die 
eine geringere Energie besitzen, bevorzugt werden. Das bedingt aber die 
Verwandtschaft der in den Fällen ZnS und SiC gefundenen neuen Strukturen. 
Auch die Tatsache, daß im Fall ZnS und SiC zum Teil die gleichen hoch- 
periodischen Strukturen auftreten, wird damit zwangsläufig erklärt, weil 
die Schwingungsentropie der praktisch gleichen Strukturverhältnisse einen 
sehr verwandten Verlauf besitzen muß. 

Schon die experimentelle Untersuchung des eindimensionalen Modells 
und der Vergleich mit den theoretischen Überlegungen haben gezeigt, daß 
alles das, was die Theorie vom heutigen Standpunkt zur Klärung des Fehl- 
ordnungsproblems im Kristall beitragen kann, im Grunde genommen nur 
an Spezialfälle geknüpft ist. Das qualitative Verhalten der thermodynami- 
schen Eigenschaften wird zwar richtig wiedergegeben, die quantitative Über- 
einstimmung ist dagegen wenig befriedigend. Trotzdem darf nicht verkannt 
werden, daß mit Hilfe der Theorie Re. Verständnis für das Verhalten der 
Kristalle erhalten wurde. Für die kristallstrukturellen Eigenschaften wurden 
die vorhandenen Ansätze noch nicht so weit entwickelt, daß man mit ihnen 
an die exakte Deutung der Röntgenaufnahmen dreidimensional fehlgeord- 
neter Kristalle herangehen kann. Aber immerhin hat die Untersuchung des 
eindimensionalen Modells wichtige Fingerzeige gegeben, in welcher Richtung 
die Theorien erweitert werden müssen. Für Fälle, in denen Gitterverzer- 
rungen auftreten, sind wir nach wie vor darauf angewiesen, den Fehl- 
ordnungsvorgang unter voller Ausnutzung allen experimentellen Materials 
zu erforschen. Die Theorie hat zwar schon viel geleistet, aber es ist wenig 
im Vergleich zu dem, was noch erreicht werden muß. 
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IV. Zusammenfassung 


Im ersten Teil dieser Arbeit wurden die gesamten theoretischen Ansätze, 
die von der thermodynamischen Seite zur Lösung des Fehlordnungsproblems 
gemacht wurden, diskutiert. Allen Theorien ist gemeinsam, daß sie nur auf 
solche Fälle angewendet werden können, bei denen keine Störungen der 
Gittergeometrie auftreten. 

Der BraaG-WiırLıam’sche Ansatz hat die große Schwäche, daß er den 
Kristall von vornherein in richtige und falsche Gitterplätze einteilt, ohne 
sich darum zu kümmern, ob diese Einteilung physikalisch sinnvoll ist oder 
nicht. Diese Unvollkommenheit wird in der BETHE’schen Theorie, welche 
die Kristallenergie aus den Wechselwirkungsenergien errechnet, vermieden. 
Dadurch kann in vielen Punkten eine Verbesserung erzielt werden, aller- 
dings ist die Durchführung der Rechnung nur mit Hilfe von Näherungs- 
methoden erreicht worden. 

KIRKwoonp’s Semi-Invariantenmethode beruht auf Reihenentwick- 
lungen der thermodynamischen Funktionen. Die abgeleiteten Formeln sind 
allgemein gültig; für die Ausrechnung der Mittelwerte der Potenzen der 
Energie muß aber auf ein Näherungsmodell zurückgegriffen werden. Der, 
große Vorteil liegt in ihrer Verallgemeinerungsfähigkeit und Exaktheit, wenn 
die Reihenentwicklung genügend weit durchgeführt wird. Die Berechnung 
hoher Glieder ist aber eine mühevolle Arbeit. 

Die machtvollste Theorie ist wohl für das dreidimensionale Fehlord- 
nungsmodell die quasichemische Methode; sie verdankt das vornehmlich 
ihrer mathematischen Einfachheit. Die Ableitung der Formeln beruht auf 
einer Gleichgewichtsbetrachtung der einzelnen Bindungspaare, in Analogie 
zum Massenwirkungsgesetz der Chemie. Sie ist in ihrer elementaren An- 
wendung in der Güte mit der BETHE’schen ersten Näherung gleichwertig. 
Die Genauigkeit kann aber noch verbessert werden. 

Die Methode der Cluster-Integrale geht auf die statistisch-thermo- 
dynamische Behandlung realer Gase zurück und zählt rein kombinatorisch 
die Anhäufungsmöglichkeit von Partikeln in Gruppen zu 1, 2, 3 Teilchen 
usw. ab. Die Methode kann aber auch auf Kristalle übertragen werden. In 
diesem Fall werden die Integrale durch Summen ersetzt, die wesentlich 
einfacher auswertbar sind, weil die Anzahl der Konfigurationsmöglichkeiten 
beschränkt wird. j 

R. BECKER hat einen Ansatz ausgeführt, der das Problem mehr von 
der kinetischen Seite beleuchtet. Durch einfache Betrachtungen von Ober- 
flächen- und Volumenenergien in Abhängigkeit von der Temperatur können 
wichtige Aussagen über das Problem der Keimbildung bei Entmischungen 
gemacht werden. 

Die vom thermodynamischen Standpunkt vollkommenste Methode ist 
zweifellos die gruppentheoretische Methode. Diese berechnet die Zustands- 
summe aus den Eigenwerten einer bestimmten Matrix. Das Matrix-Eigen- 
wertproblem ist zwar für praktisch unendlich große Matrizen schwer zu lösen. 
ONSAGER ist jedoch die exakte Durchführung der Lösung für ein zweidimen- 
sionales Modell gelungen. Das ist ein besonderer Gewinn, weil man damit 
einen guten Maßstab für die Güte der übrigen Näherungen besitzt. 


{ 
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Im zweiten Teil wird das Fehlordnungsproblem von der röntgeno- 
graphischen Seite kurz beleuchtet. Die Deutung der Beugungsbilder drei- 
dimensional fehlgeordneter Kristalle steckt praktisch noch in den Kinder- 
schuhen, hier existieren nur einige Ansätze, bei denen man jedoch selbst 
der qualitativen Deutung noch mit großer Vorsicht begegnen muß. Nur für 
zweidimensional und eindimensional fehlgeordnete Kristalle kann das Beu- 
gungsproblem befriedigend für den Fall des thermodynamischen Gleich- 
gewichts gelöst werden. Für eindimensional fehlgeordnete Kristalle wird an 
der Durchführung des bisher noch nicht berechneten Falles von A-, B- 
Statistik mit Wechselwirkungen bis zu den übernächsten Nachbarn gezeigt, 
wie man die Wechselwirkungsenergien und die Temperatur zur Errechnung 
der röntgenographisch allein interessierenden absoluten Fehlerwahrschein- 
lichkeiten heranziehen kann. Damit wurde die Brücke zu den thermo- 
dynamischen Berechnungen geschlagen. Die Prüfung der Theorie mit dem 
Experiment stößt aber schon im eindimensionalen Fall besonders bei hohen 
Temperaturen auf ernsthafte Widersprüche, die darauf hinweisen, daß die 
hier durchgesprochenen Berechnungen nicht immer anwendbar sind. Für 
ZnS und SiC wurde wahrscheinlich gemacht, daß das in diesen Fällen auf- 
tretende anomale Verhalten durch die Vernachlässigung der Schwingungs- 
entropie bei der Berechnung der freien Energie hervorgerufen wird. Nimmt 
man an, daß die Dichte der Energieniveaus der Schwingungszustände im 
geordneten Zustand wesentlich höher liegt als im ungeordneten, so kann 
man eine qualitative Deutung aller anomalen Eigenschaften in diesen beiden 
Fällen herbeiführen. 


Marburg, Mineralogisches Institut der Universität. 
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Versucht man die Phosphate (Arsenate, Vanadate) auf kristallchemi- 
scher Grundlage systematisch zu ordnen, so stößt man auf beträchtliche 
Schwierigkeiten (vgl. 96). Diese Tatsache dürfte vielleicht überraschen, da 
die Isotypiebeziehungen zwischen Phosphaten und Silikaten eine analoge 
kristallchemische Klassifikation nahelegen könnten. Diese Schwierigkeit 
mag an der Doppelrolle liegen, die das P5+-Ion als Zentralkation spielt. Im 
Sinne der Definition gehören die Phosphate zu den anisodesmischen Gitter- 
typen (20). Die enge kristallstrukturelle Verwandtschaft zu den Silikaten 
jedoch offenbart eine deutliche mesodesmische Tendenz. Tatsächlich liegen 
Verknüpfungen von PO,-Tetraedern vor, wie z.B. im Mg,[P,0,]. 

Einleitend sind demnach die beiden Fragen zu klären: Wie werden 
1. die [4]-Koordination und 2. die Doppelrolle des PO,-Tetraeders als aniso- 
-desmischer und mesodesmischer Typus kristallchemisch verständlich ? 


Wir werden versuchen, die erste Frage zunächst auf energetischem 
Wege zu lösen. Tatsächlich führt schon eine einfache Betrachtung zu dem 
Verständnis, daß PO,-Komplexe wahrscheinlich sind (nach Koss, vgl. (1) 
und (72)). Q sei die Energie, die frei wird, wenn das Komplexion AX, aus 
einem Molekül (AX,„) und einem oder mehreren einfachen Ionen (X'-) ge- 
bildet wird. Ist ferner Ey, die Energie, die der Bildung eines Moleküls aus 
seinen Ionen entspricht, und Ex die Energie, die bei der Bildung des Kom- 
plexes aus den konstituierenden Ionen frei wird, dann gilt die Beziehung: 
Q = Ex — Ey. Rein ionogene Bindung vorausgesetzt, können die Ex und 
Ey in einfacher Weise berechnet werden, wobei nur COoULoMB’sche Anzie- 
hung und Abstoßung berücksichtigt zu werden braucht. Wird @ ein Mini- 
mum für einen Komplex AX,, so ist dieser als wahrscheinlichster zu er- 
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warten. Führt man die Rechnung für zentrale A-Ionen mit den Ladungen 
1, 2, 3, 4, 5, 6 und 7 und für einwertige und zweiwertige X-Ionen aus, so 
ergibt sich für die zu erwartenden Koordinationen folgende Tabelle: 


Ladungen des zentralen A-Ions 1 2 3 4 5 6 7 
KZ für einwert. X 2 4 5 6 8 8 8 
KZ für zweiwert. X l 2 3 4 4 5 6 


Nies (72) hat diese Betrachtungsweise noch weiter vertieft: Sei A 
der elektropositive, X der elektronegative Partner eines binären Kom- 
plexes mit den entsprechenden Wertigkeiten w, und wx, so nennt NIGGLI 
den Ausdruck w‚/wx die Bindungszahl. Bei heteropolarer Bindung hat 
diese Größe die Bedeutung, daß im Mittel ein A-Ion wı/wx X-Ionen direkt 
zu binden vermag. Betrachten wir wieder die Inselradikale AX,, und es sei 
die Bindungszahl für sie w4/wx = m, die Koordinationszahl n, so ist in 
Bruchteilen der Wertigkeit des Liganden jeder Ligand zu m/n abgesättigt. 
Nun wird die Annahme postuliert, daß einkernigen Komplexionen norma- 
lerweise nur dann eine größere Beständigkeit zukommen kann, wenn die 
Liganden mindestens halb an das Zentralion gebunden sind, d. h. es muß 
seinn/m > #. Demnach ist n = 2m die normalerweise zu erwartende maxi- 
male KZ für stabile einkernige Radikale. Damit ergibt sich folgende Tabelle: 


Wertigkeit des Zentralions il 2 3 4 5 6 U 8 
Bindungszahl b. einwert. X 1 2 4 5 6 7 8 
maximal zu erwart. KZ 2 4 6 3, 1022 21416 
Bindungszahl b. zweiwert. X ı | ee a ae 
maximal zu erwartende KZ 1 2 3 4 5 6 7 8 


Man erkennt sofort, daß beide Tabellen beträchtliche Widersprüche zur 
Erfahrung ergeben. So ist z. B. nach der ersten Tabelle für 6wertige Ka- 
tionen [5]-Koordination zu erwarten, während dies etwa für S$+, Cr6+ usw. 
nicht zutrifft. Diese Ausnahmen können auf zwei Gründe zurückgeführt 
werden: Bei der Rechnung ist die Voraussetzung benutzt worden, daß sich 
die Moleküle durch Vereinigung freier Ionen bilden. In Lösungen müßte die 
Hydratationsenergie berücksichtigt werden. Dadurch wird offenbar das 
Energieminimum nach Komplexen mit kleinerer Ionenzahl verschoben, 
d. h. die KZ nimmt kleinere Werte an. 

Der zweite Grund ist in der Raumbeanspruchung gegeben: Der Be- 
reich einer KZ wird bekanntlich durch das Radienverhältnis begrenzt. Des- 
wegen treten statt der Komplexe NO, CO, SO, 
tatsächlich folgende auf NICH I0 

Dagegen verhalten sich SiO,*" und PO,?- durchaus, wie sich auf Grund 
der 1. Tabelle erwarten läßt. 

Die geometrische Existenzbedingung der [4]-Koordination lautet: 


0,225 S r,/rx < 0,414 für tetraedrische und 0,414 < ry/rx < 0,732 für | 


quadratische Konfiguration. Letztere wollen wir aus unserer Betrachtung 


ausschließen. Folgende Tabelle gibt Ionenradien und Radienguotienten (be- 
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zogen auf O°-) für eine Reihe von Kationen A, die kristallchemisch als Kerne 
von [AX,]-Tetraedern auftreten können (Radien nach PAuLıng (75a)). 


A r,mÄ ra/to A r,inÄ rufto 
Be?* 0,31 0,222 ser 0,29 0,207 
AI*+ 0,50 0,357 Sest 0,42 0,300 
Si4t 0,41 0,292 Cr 0,52 0,372 
Get 0,53 0,378 Mos* 0,62 0,443 
Irtt 0,66 0,472 (leer 0,26 0,186 
Os 0,67 0,478 Bros 0,39 0,278 
we 0,68 0,486 kit 0,50 0,358 
Pas 0,34 0,243 Mn’+ 0,46 0,328 
Asst 0,47 0,336 
Vor 0,59 0,422 


Aus der Tabelle ergibt sich eindeutig, daß für P tetraedrische Koordina- 
tion zu erwarten ist. Lediglich liegt das Vanadium mit seinem Quotienten 
bereits im Übergangsbereich, immerhin noch so nahe der geometrischen 
Grenze, daß seine Zugehörigkeit zum Tetraedertypus verständlich bleibt. 
Neben Ge und Si fallen die übrigen vierwertigen Kationen heraus, insbeson- 
dere sollte bei ihnen keine tetraedrische Koordination zu erwarten sein. Sie 
zeigen daneben auch keinerlei Tendenz zur mesodesmischen Vernetzung, 
so daß eine vom heteropolaren Typus stark abweichende Bindung anzu- 
| nehmen ist. Hier erklärt also schon der beträchtliche Größenunterschied 

zum Sit+-Ion das völlig andersartige Verhalten. Dagegen zeigen einerseits 
ı das Al®+, andererseits das P*+* usw. sowohl hinsichtlich der absoluten als 
‚ auch relativen Dimensionen eine weitgehende Verwandtschaft zum Sit. 
' Offenbar spielt, wenn wir von der Gruppe Ir—W absehen, auch die 
‚ elektrostatische Ladung eine Rolle. Tatsächlich können wir annehmen, daß 
‘die Bindungsverhältnisse weitgehend verwandt sind, wogegen zum W*+ 
‘oder auch zum S$+ bemerkenswertere Unterschiede vorliegen. Auch die 
' Elektronenstrukturen der entsprechenden Atome, die auf verwandtschaft- 
‚liche Tendenzen im homöopolaren Bindungsanteil hinweisen können, deu- 
"ten die erwähnten Beziehungen an. Wenn wir nur die Valenzelektronen 


i betrachten, so haben wir für Si das Schema 3s?3p?, dem das Schema für Ge 


‘ mit 4s?4p? entspricht. Aluminium mit 3s?3p!, Phosphor mit 3s?3p®? und 
‘ Arsen 4s?4p? zeigen gleichartige Konfigurationen. Dagegen fallen Wolfram 
_ mit 5d?6s?, Osmium mit 5d66s? und Molybdän mit 4d?5s! heraus. Wider- 
-spruchsvoll bleibt in diesem Zusammenhange vielleicht nur das Vanadium, 
das eine äußere Elektronenkonfiguration mit dem Schema 3d?4s? besitzt. 

Wenn wir nach dem Vorbild von PauLine (75a) die Partikelabstände 
für die wichtigsten tetraedrischen Gruppen zusammen betrachten, so wird 


| das oben Gesagte noch klarer: 


| Ww—-o Se) PO SO) (le) 
in WO4- in SIO4 POT in SO2- in 010,7 
ra+troinÄ 2,08 1,81 1,74 1,69 1,65 
beob.d.in Ä 1,73 1,60 1,55 1,51 (1,48) 
A 0,35 021 0,19 0,18 0,17 
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Wesentlich ist hier wiederum, daß SiO, und PO, unmittelbar neben- 
einander stehen. Deutlich weicht dagegen WO, ab, während wir zwischen 
PO,, SO, und C10, offenbar einen mehr allmählichen Übergang vorfinden. 
Wir verstehen aus diesem Verhalten ohne weiteres die enge kristallchemi- 
sche Beziehung zwischen Phosphaten und Sulfaten. 


Daß andererseits bei SO, keine mesodesmische Tendenz zu erwarten a 
ergibt sich aus den Werten für die elektrostatische Bindungsstärke p = 
in folgender Tabelle: 
5% BeO, AIO, 
p=2.. 050 0,75 


Sio, 
1 


PO, so, c10, 
1,25 1,50 1,75 


Auch aus dieser Tabelle ergibt sich sehr gut der Zusammenhang zwi- 
schen Al — Si—P! 

Die Verkürzung der Atomabstände gegen die berechneten Radiensum- 
men weist nach PAuLInG (75a) auf Anteile einer kovalenten Doppelbindung 
hin. Danach können bei Berücksichtigung von Doppelbindungen folgende 
Strukturformeln aufgestellt werden: 


RO _— Om m 
ae le ae Rh 
ER 


Jede dieser Doppelbindungen soll einer Resonanz zwischen den drei 
Strukturen A::Ö:, A+ :Ö:- und A++0-- entsprechen. Allerdings beinhal- 
ten die angegebenen Strukturformeln nicht die konventionellen Wertig- 
keiten der Zentralatome. So ist z. B. das Chloratom nicht als siebenwertig 
angesetzt. Um auf die Konvention Rücksicht zu nehmen, schlägt PAULING 
folgende Formulierungen vor: 


eu re oO 
| | | | 
Co "0—P=0 ”—8=-0 —-a=-0| | 
| | | | 
0- 0: 0 oO 


Dabei könnten die entsprechenden Resonanzstrukturen in Klammern ı 
angefügt werden, falls eine präzisere Beschreibung der Gruppierungen not- 
wendig erscheint, z. B.: 


(Ol WS (On 

| | 
20 SZ -o0—s7 =0 0=-St=-0 usw 

| | 

Oz Ol On 


Im allgemeinen wird ein derartiger Komplex Struktur und Dimensio- - 
nen auch in verschiedenen Kristallgittern beibehalten, da die Bindungs- 
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kräfte innerhalb eines Komplexes gewöhnlich stabiler sind als zwischen 
Komplex und umgebendem Kation, zumal wenn diese geringe Ladungen 
aufweisen. PAULING weist in diesem Zusammenhange auf eine bemerkens- 
werte Ausnahme hin. Im Ag,PO,-Gitter beträgt der P—O-Abstand 1,61+ 
0,03 Ä im Gegensatz zu 1,55 A in anderen Phosphaten (HELMHoLzZ, L.: 
J. Chem. Phys. 4. 1936. 316). Man kann sich diese Erscheinung etwa fol- 
gendermaßen erklären: Im Ag,PO, ist jedes Sauerstoffatom tetraedrisch 
von einem P- und 3 Ag-Atomen umgeben. Die Bindungen der Silberatome 
weisen kovalente Tendenzen auf. Da das Sauerstoffatom nur 4 kovalente 
Bindungen herstellen kann, muß der kovalente Charakter der Ag—O-Bin- 
dungen z.T. den Doppelbindungscharakter der P—O-Bindung verhindern. 
Konsequenterweise führt dies zu einer Ausweitung des PO,-Komplexes. 


A. Anisodesmische Analogien 
Krokoit-Typus 


VON GLISZCZYNSKI (23) hat die Isotypiebeziehung zwischen Krokoit, 
Monazit und Xenotim röntgenographisch genauer untersucht. Er kommt 
zu dem Ergebnis, daß eine weit größere kristallographische Verwandtschaft 
zwischen Xenotim und Monazit als zwischen Xenotim (oder auch Monazit) 
und Krokoit anzunehmen ist. Nach VON GLISZCZYNSKI kann vermutet wer- 
den, daß gewisse Monazite magmatischen Ursprungs zur Zeit ihrer Bildung 
eine deformierte Xenotimstruktur besitzen, die erst allmählich in die eigent- 
liche Krokoitstruktur übergeht. Dagegen weisen die Turnerite von. vorn- 
herein Krokoitstruktur auf. Die Beziehungen der Strukturen ergeben sich 
auf Grund folgender Zusammenstellung, wobei die Zirkonstruktur noch 
miterwähnt sei: 


Ei bo Co ß Rer. 
Zirkon Zr SiO, 6,58 (6,58) 5,93 Di) — I4/amd 
CaCrO, — - — Zirkonstruktur 
Xenotim Y PO, 6,88 (6,88) 6,03 Di, — I4/amd 
Turnerit Ce PO, 6,782 6,993 6,455 103° 40’ en — P2,/e 
Monazit Ce PO, 6,78 6,99 6,45 104° @, — P2/n 
BiPO, 6,74 6,95 6,41 104 0, b2,n 


Krokoit Pb CrO, 7,1068 7,410 6,771 102037 On—P2n 


Bei geeigneter Kombination der Formelschreibweisen von MACHA- 
TSCHKI (60) und Strunz (95) erhalten wir für Zirkon den Formeltypus 
X[3[z0,] und für den Krokoittypus X! [z0O,]. Bei kleinerem Radius für 
das X-Ion wird der Zirkontypus, bei größerem Radius der Krokoittypus 
gefunden. 

Nach ZEMANN (113) besteht Isotypie zwischen Monazit und BiPO,. 
Die Gitterkonstanten sind für das Wismutphosphat etwas kleiner anzu- 


setzen, wobei allerdings die maximale Abweichung 1% kaum überschreiten 
12* 
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dürfte. Auch die Parameter für BiPO, sind ähnlich wie beim Monazit. Die 
Bi-Atome sind im BiPO, von 10 Sauerstoffatomen unregelmäßig umgeben. 
Die Bi—O-Abstände liegen zwischen 2,40 und 2,98 Ä. Die Phosphoratome 
sind fast regelmäßig von 4 Sauerstoffatomen in Abständen zwischen 1,48 
und 1,51 Ä umgeben. Nach Mooney (71) hat die tetragonale Modifikation 
von BiAsO, Scheelit-Struktur mit [8]-Koordination um Bi. 


Palmierit-Typus 


_ Formeltyp: X,[z0,].- 
Nach Strunz (93) sind Whitlockit und Palmierit isotyp, nicht iso- 
morph. Folgende Zusammenstellung der Gitterdaten dürfte diese Beziehung 
belegen: 


Palmierit PbK;[SO,], Whitlockit Ca,[PO,], 
2029,48 10,32 = 2 x 5,16 
er 20,61 36,9 =2 x 18,45 
c0208,161 3,976 
MR ZA) 
ame 11.56 13,65 = 2x 6,825 
— 420 28 44° 06° 
Da (?) ei) 


Die Daten für PbK,[SO,], sind an synthetischem Material gewonnen 
worden. FRONDEL (21), der Whitlockit röntgenographisch untersucht hat, 
hat darauf hingewiesen, daß strukturell keine Verwandtschaft zum Apatit 
vorliegt. Immerhin liegt eine gewisse Ähnlichkeit in den Gitterdimensionen 
vor: | 


Whitlockit Fluorapatit 
Raumsruppe we... BR e C 6,/m 
a0 (hex.Koord. 2.2 20.521032 9,37 
EN 0 36,9 6,88 (36,9/6) 


Alunit-Typus 


Die Kristallarten sind ausgesprochen pseudokubisch. Als Formeltypus 
kann folgender angegeben werden: 


WIs+sY,[S[(OH), | (20,).]- 


__ Nach GoOSSNER (24) kann wegen des Ersatzes von SO, durch PO, (OH) 
beim Hamlinitgitter die Raumgruppe D3, nicht mehr vorliegen, sondern 
nur eine Raumgruppe ohne Symmetriezentrum. Als möglich wird C}, als 
Untergruppe von D;, angegeben. Wegen des pyroelektrischen Verhaltens 
folgert HENDRICKS (32) auch für Alunit die Raumgruppe C}, — R3m. Nach 
LEmmon (52) können im einzelnen drei Reihen unterschieden werden: | 
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Mineralname W(X) Y z ar Conkcolaoı cis a 
Alunit K Al Ss 6,96 | 17,35 | 1,246 | 1,252 [89010 
2 Natrojarosit Na Fe S 7,18 | 16,30 | 1,135 1,104 93 50 
‘3 Jarosit K Fe Ss 7,20 | 17,00 | 1,181 | 1,249 |89 15 
3 Ammoniojarosit NH, Fe Ss 7,20 | 17,00 | 1,181 = z0525 
= Argentojarosit Ag Fe Ss 7,22 | 16,40 | 1,136 | 1,106 |93 42 
Plumbojarosit Pb Fe S 7,16 | 16,90 | 1,180 — —— 
synthetisch H,0 Fe | Ss 7,16 | 16,90. L180 1 = — 
= Woodhouseit Ca Al Ss? | SB ee = —_. 
& Svanbergit Sr Al SP BIN) | — — 
= Hinsdalit Pb Al SP 1,2677 |89 40 
S Corkit Pb Bol SP 1,1842 |91 16 
S Beudantit Pb Fe SAs 1,1842 '91 16 
Pseudowarvellit Ca AL) SB = 
© Deltait Ca AH) P 6,98 | 16,10 | 1,154 | — == 
3 Hamlinit Sr AUF) EBEN TE 890, u SAT S8 13 
2 Gorceixit Ba AH), P — 
3 Dussertit Ba |Fe(H)| As | | — 
5 Plumbogummit | Pb AH) P — 
Floreneit Ce Al 12 1,1901 |88 56 
Stiepelmannit |Ce,Y,Yb| Al | P | 6,75 | 16,52 | 1,224 | 1,272 |ss 32 


Nach LARSEN (51) zeigen Pulveraufnahmen von Pseudowavellit und 
Deltait sehr ähnliche Gitterdimensionen. Unterschiede wurden lediglich in 
der chemischen Zusammensetzung der beiden Minerale gefunden. Deswegen 
vermutet LARSEN, daß beide Kristallarten eine isomorphe Mischkristall- 
reihe bilden, wobei Ca gegen Al mit OH auszutauschen ist. Nach STRUNZ (95) 
stehen Alunit- und Beudantitreihe im gleichen Symmetrieverhältnis wie 
Kalkspat und Aragonit. Demnach handelt es sich bei den Kristallarten der 
Beudantitreihe um ‚„Doppelsalze‘“ (nicht um Mischkristalle!). Zwischen 
der Alunit- und Hamlinitreihe besteht nach STRUunz (95) das Verhältnis der 
Homöotypie. 

Gipstypus 

Formeltypus: X'®[z0,]: 2 H,O. Nach Strunz (92) sind Weinschenkit 
und Gips isotyp. Bei der folgenden Zusammenstellung ist die veränderte 
Aufstellung von Gips zu beachten. Bei der gewöhnlichen Aufstellung ist die 
Elementarzelle allseits flächenzentriert. Nach der Transformation [101]>aneu 
und [001] — Cneu erhält man eine einfachere basiszentrierte Zelle (a0 bo’ 
co ß). Außerdem wird Z = 4 statt 2 = 8. 


Gips Weinschenkit 
CaSO, :-2H,0 (Y, Er)PO, :2H,0 
8508,51 a0 = 6,48 
bo 1515 bo, = 15,12 


6,28 Co = 6,28 
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Gips Weinschenkit 
RE 5 DE ER I a0:b5 205 = 0,4298717.10,415 
ß’ = 1279 24° ß = 129° 24° 
C—A2/a CL, —A2/a(?) 
Zu rt Z=4 


Ferner besteht Isotypie zwischen Gips und Brushit bzw. Pharmakolith 
bzw. Ardealit. Dabei tritt an Stelle des Caleiums CaH in das Gitter ein. Die 
Beziehung zwischen Gips, Brushit und Ardealit wurde von HALLA (30) er- 
kannt. Den Hinweis auf die Verwandtschaft zwischen Gips und Pharma- 
kolith verdanken wir GoOSSNER (24). Ardealit ist als ‚Doppelsalz‘‘, nicht als 
Mischkristall aufzufassen. 


Gips Brushit Pharmakolith Ardealit 
Ca[SO,]:2H,0 CaH[PO,] :2H,0 CaH[As0,]-2H,0 CaH[PO,S0,] - 4H;! 
0 10,7 10,47 10,97 10,47 
b, = 15,15 15,15 15,40 14,64 
an = 08 6,28 6,29 6,28 
bo :C0 = 0,69°71°20,415 ebenso 0,712 :1 : 0,409 0,715 :1 : 0,429 
B2— 9811587 9515 96 36 979 
Zu 8 8 8 4 


B. Übergangsgruppe und Isomorphie 


Eine Reihe von Phosphaten mit der Tendenz zur Bildung anisodesmi- 
scher Inselstrukturen weisen eine ausgesprochene Neigung zu Isomorphien 
auf, wobei P insbesondere durch Si und/oder S diadoch vertreten wird. 
Solche Beziehungen scheinen offenbar dort besonders ausgeprägt zu sein, 
wo Strukturen mit intermediärem Charakter (zwischen anisodesmisch und 
mesodesmisch stehend) aufgebaut werden. Hierher gehören zunächst ein- 
mal Alkali—Erdkali-Phosphate vom Formeltypus X,[z0,] bzw. W,[z0,]. 
Nach KLEMENT und STECKENREITER (45) kommt die Verbindung NaCaf[PO,] 
in zwei Modifikationen vor, wobei die -(Tief-)Modifikation rhombisch, die 
a-(Hoch-)Modifikation hexagonal kristallisiert. KLEMENT und STECKEN- 
REITER haben bei diesem Phosphat die Ersatzmöglichkeit Ca—Sr, Ba und 
P—Si, S untersucht. Dabei zeigte es sich, daß unter Beibehaltung der 
Struktur von NaCa[PO,] leicht P durch Si + S ersetzt werden kann. Be- 
vorzugt bildet sich bei diesen Mischverbindungen die Hochtemperaturform. 
Folgende Mischungen wurden dargestellt: NaCa[PO,], NaCa,[Si,PO;>], 
Na,Ca,[P;SO,,], NazCa,[SiP,O,5], Na;Ca; [SiPSO,,], Na,Ca [PS,0, 3]. 

Im K—COa-Phosphat ist nach KLEMENT und STECKENREITER der Er- 
satz des P durch Si-+-S nicht in allen Fällen ohne Strukturänderung 
durchführbar. 

Mineralogisch von Bedeutung sind zwei Phosphatgruppen, die weit- 
gehende Diadochie in der angegebenen Richtung zu erkennen geben. Es 
handelt sich um die Wagnerit- und Apatitgruppe. Als bemerkenswert ist 
hervorzuheben, daß in beiden Fällen die Phosphate und deren Strukturen | 
zuerst bekannt waren und erst später die entsprechenden (P, Si, S)-Misch- 
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kristalle auf synthetischem Wege gefunden wurden. Die Zusammensetzung 
der beiden Gruppen kann kristallchemisch wie folgt formuliert werden: 
Wagnerit Y,[A| zO,] 
Apatit X,9I[A| (z0,)3]. 
In Erweiterung der MACHATSCHKT’schen Symbolik steht hier für F, Cl, 


) OH, 0, CO, usw. A. 


KLEMENT und ZUREDA (46) geben als Grenze des Kationenradius für 
das Auftreten der Wagnerit- bzw. Apatitstruktur 1.0 Ä an. Für rx > 10Ä 
findet sich der Apatittypus (z. B. Cd, Ca, Sr, Pb, Ba), dagegen fürry <1.0Ä 
der Wagnerittypus (z. B. Mn, Cu, Zn, Mg). Mit Wagnerit isomorph sind 
Triplit, Triploidit und Sarkinit (Mn,[OH|AsO,]). KLEMENT und GEMBRUCH 
(44) haben die zu Wagnerit isomorphen Mischkristalle LiMg,[F, IPSO,] und 
NaMsg;[F, | PSO,] gefunden. 

Noch wesentlich umfangreicher sind die Isomorphiebeziehungen inner- 
halb der Apatitgruppe. Nach McConxerL (64) bestehen folgende Substi- 
tutionsmöglichkeiten: X — Ca, Mg, Mn, Sr, K,Naund C (\);z=PP, As, V, 
Si, Sund © (!) und schließlich A = F, Cl, OH, Ö (nicht aber CO,!). Nach 
zahlreichen Untersuchungen, insbesondere von McConnELL und Mitarbei- 
tern (63—68) und KLEMENT und Mitarbeitern (40—43) sind Mischkristalle 
mit weitgehend diadocher Vertretung von P durch Si+S möglich. Als End- 
glied der (P, Si, S)-Mischkristallreihe kann der von MCConnNELL (63) zuerst 
untersuchte Ellestadit Ca,,[F's | (SiSO,),] aufgefaßt werden. Zwischenglieder 
im (P, Si, S)-Mischbereich sind Wilkeit Ca,[(F, O) | (PO,, SiO,, SO,);] und 
Britholith (Na, Ce, Ca),[(F, OH) | (SiO,, PO,),]. Die Isomorphie des letzt- 
genannten Minerals zu Apatit wurde von MACHATSCHKI und HÄGELE (29) 
röntgenographisch nachgewiesen. Der Ersatz von P°* durch Si?* und der von 
Ca?+* durch Ce?+ bedingt eine leichte Aufweitung des Gitters gegenüber 
Apatit. In folgender Tabelle sind natürliche und synthetische Glieder der 
Apatitreihe mit den entsprechenden Gitterdaten zusammengestellt: 


Mineralname X A zZ &o Co Co/&o 
ANsehne Ca ICLOH IB, 9,36 | 6,88 | 0,732 
Spaneh a a oe Ca F SiP,S 9,45 | 6,96 | 0,736 
Spann Ca OH SiP,S 9,44 | 6,96 | 0,737 
Dilketemaee un. Ca F,O Si,P,S | 9,48 | 6,91 0,729 
Billestadita an ar. Ca F Si,S 9,54 | 6,99 | 0,732 
Hydroellestadit. ... . . (a OH Si,S 9,54 | 6,99 0,732 
BritholcHe Na,Ce,Ca|)| F,OH IS, 9,61 | 7.022°.07730 
Krancolthes mer Ca F, OHO 12510. 9,34 | 6,89 | 0,738 
EynthWwe en Na,Ca, F S 9,49 | 6,87 | 0,724 
Sale, Se NaCa, F P,S 9,52 } 6,90 | 0,725 
Salat. au ee Ca F As — _ (0,714) 
Byromorphiu a mr Pb Cl ir 9,95 | 7,32 | 0,736 
Nimetestti arts. 0m Pb (el As 10,36 | 7,52 0,726 
Nanadıntme er: Pb [61 V 10.4728 7.432 0/7110 


1 morphologisch 
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Nach Strunz (Naturwiss. 30. 1942. 65) stellt das Mineral Steenstrupin 
ein Silikat vom Formeltypus Apatit dar, bei dem ein Anteil von 0—20 % Si 
durch P vertreten sein kann. Nach Mitteilungen von MACHATSCHKTI (62) 
jedoch ist Steenstrupin nicht mit Apatit verwandt, und seine Struktur 
ist unbekannt. 

Heftig erörtert wurde die Frage der Karbonatapatite (3, 5—7, 65—68). 
Neuerdings hat TH. GEIGER (22) eine wesentliche Untersuchung zu diesem 
Problem veröffentlicht. Aus den GEIGER’schen Resultaten ist zu entneh- 
men, daß letzten Endes immer noch nichts Sicheres über die Konstitution 
der Karbonatapatite gesagt werden kann. Insbesondere wird darauf hin- 
gewiesen, daß die Existenz der eigentlichen Karbonatapatite im Sinne von 
McConNELL und GRUNER (68) keineswegs als gesichert angenommen wer- 
den darf. Im Gegenteil sprechen verschiedene experimentelle Befundegegen 
einen gittermäßigen Einbau des Kohlenstoffs, und gewisse Eigenschaften 
der Karbonatapatite weisen auf die Einlagerung von submikroskopischen 
Caleitteilchen in Fluor- und Hydroxylapatiten hin. Vor allem aber ist nach 
GEIGER die Möglichkeit der anomalen Mischkristallbildung zwischen Apatit 
und CaCO, in Betracht zu ziehen, d. h. es würde gegebenenfalls eine ge- 
setzmäßige Einlagerung einer CaCO,-Modifikation in Apatit bei verschie- 
denen Dispersitätsgraden vorliegen. Im folgenden geben wir die GEIGER- 
sche Übersicht über den Stand des Problems wieder: 


Karbonat-F-Apatit Karbonat-OH-Apatit 


Sehr geringe Unter- 
schiede zwischen den 


Eher für einen Einbau Pulverdiagrammen 

von C ins Apatitgitter keine Anhaltspunkte 
: Veränderung des Karb. 

sprechen: 


Apatitdiagr. in dasj. von 
F-Apatit beim Erhitzen 


Eigenschaften, die keine Erniedrigung der Licht- und Erhöhung der Doppel- 

Entscheidung zulassen: brechung gegenüber F- und OH-Apatit. Dichte- 
erniedrigung gegenüber F-Apatit. Chemische Zu- 
sammensetzung 

Eher für eine Einlagerung Sehr große Ähnlichkeit | Identität der Pulver- 

von CaCO, im Apatit zwischen den Pulver- diagr. von OH- und 

sprechen: diagr. von F- und Karb.- | Karb.-OH-Apatit 
F-Apatit | 


Bildung von freiem CaO bei der thermischen Behand- 
lung unter gleichzeitigem CO,-Verlust 


Dinn und KLEmEnT (19) kamen bei ihren Apatituntersuchungen auch 
zu solchen Typen, bei denen die Zusammensetzung nicht mehr der Ideal- 
formel X,[A | (20,),] entspricht. Diese Tatsache überrascht nicht, wenn man 
bedenkt, daß das Apatitgitter doch offenbar sehr stabil ist. Insbesondere 
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können in diesem Gitter Leerstellen auftreten oder auch überschüssige Ionen 
auf Zwischengitterplätzen vagabundieren. Folgende Mischungen mit Ca- 
Überschuß bzw. Ca-Mangel wurden hinsichtlich der Möglichkeiten des Auf- 
baus eines Apatitgitters untersucht: 

Ca1o.5 [F, | SiP;O,,] 

Ca105 [F, | Sl] 

Ca 9.5 [F, | SiP,38,0;,] 

Die synthetisch gewonnenen Produkte kristallisieren hexagonal und 
ihre Pulverdiagramme entsprechen durchaus dem des Apatits. 

In diesem Zusammenhange ist vielleicht erwähnenswert, daß nach 
PECoRA und FAHEY (77) beim Lazulith-Skorzalith (Fe, Mg)Al,[OH | PO,] 
offenbar keine (Si, P)-Mischkristalle vorliegen. Der beobachtete SiO,-Gehalt 
wird auf Verunreinigungen zurückgeführt. 

Das System Ca0—P,0,—SiO,. Im Zusammenhang mit derMisch- 
kristallbildung mit Diadochie von P/Si erscheint es wichtig, wenigstens kurz 
auf das System Ca0—P,0,—SiO, einzugehen, das vor allem von TRÖMEL 
(47, 101—103) und BREDIG (8—10) intensiv untersucht worden ist. Nach 
BREDIG (9) existieren in dem genannten System wenigstens 3 Kristallpha- 
sen, von denen keine in den binären Systemen Ca0O—Si0, oder Ca0—P,0, 
oder gar SiO,—P,0, bekannt waren. Diese sind: 5CaO -P,O, SiO, (Silieocar- 
notit), 4+CaO-P,O,-SiO, und 7Ca0-P,0,-2SiO,. BREDIG konnte plausibel 
dartun, daß die letztgenannte Kristallphase sehr wahrscheinlich nicht eine 
ternäre Verbindung, sondern eine feste Lösung von Ca;[PO,], in Ca,SiO, 
darstellt. Das Calciumorthosilikat selbst besitzt eine einfache hexagonale 
Kristallstruktur, die dem Typus der Hochtemperaturformen der Alkali- 
metallsulfate entspricht. Nach SCHLEEDE (81) besitzen Cas,(SiO,)}.(PO,)s 
und Ca,,K,(SiO,)], eine mit $-K,SO, übereinstimmende Struktur. Nach 
TRÖMEL und Mitarbeitern (103) scheint ganz allgemein zu gelten, daß sich 
das Orthosilikat 3 CaO-P,O, mit Oxyden wie SiO,, FeO, MnO usw. im flüs- 
sigen Zustand nur beschränkt mischt. Das ist deswegen besonders bemer- 
kenswert, weil die «-Modifikation von 3Ca0 -P,O, bis zu 5% SiO, unter 
Bildung recht stabiler Mischkristalle aufnehmen kann. Weiterhin bilden sich 
aus Schmelzen, die aus Calciumorthophosphat und Calciumorthosilikat zu- 
sammengesetzt sind, in weitem Umfange Mischkristalle. 


C. Mesodesmische Analogien 
Olivin-Typus 

Formeltyp: Y,[z0,]. 

Als erstes Phosphat mit Olivinstruktur wurde von ZAMBONINI und 
Marossı (112) der Lithiophilit Li(Mn, Fe)[PO,] untersucht. Dann folgte die 
Arbeit von ZAMBONINI und Lavzs (111) über das Lithiumorthophosphat 
Li,PO,. Anschließend berichteten GoSSNER und STRUNZ (27) über den Tri- 
phylin Li(Fe, Mn)PO,. Zur Triphylingruppe gehört außerdem der Natro- 
philit Na(Mn, Fe)PO,. T#ıLo (99) hat aus Mischungen von YNH,PO,: H,O 
und Li,CO, durch Erhitzen auf 800° C eine Reihe von Verbindungen der 
Zusammensetzung YLiPO, synthetisch hergestellt. Dabei zeigte sich, daß 
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die Ni-, Mg- und Co-Verbindung völlig dem Olivin entspricht. Die Zn-Ver- 
bindung zeigt nach THıLo ein ähnliches Diagramm. Mn- und Cu-Verbindung 
entsprechen dem Lithiophilit, während die Cd-Verbindung zwischen Tri- 
phylin und Monticellit liegt. Schließlich entspricht das CaLi[PO,] dem 
Monticellit. 

In folgender Zusammenstellung sei die besprochene Isotypiebeziehung 
an Hand der Gitterdaten erläutert: 


30 bo Co BEIDE:TC 
Olivin 6,00 4,77 10,26 1,256 :1 :2,151 
(Mg, Fe),[SiO,] 
Triphylin 6,00 4,67 10,34 1,285 : 1 : 2,214 
Li(Fe, Mn)[PO,] 
Lithiophilit — — — 1,290 :1 :2,211 
Li(Mn, Fe)[PO,] 
Lithiumphosphat 6,07 4,86 10,26 1,2495: 1e7 2511 
Li,[PO,] 
Natrophilit 6,32 4,97 I ea 


Na(Mn, Fe)[PO,] 

Es ist auffallend, daß das Li,PO, ein Y-Kation mehr hat als der Olivin. 
Auf Grund des strukturellen Aufbaus des Olivins kann man sich die Be- 
ziehung einfach klarmachen: Im Olivin bilden bekanntlich die Sauerstoff- 
ionen eine schwach deformierte hexagonal dichteste Kugelpackung. Die Te- 
traederlücken werden z. T. von Si-Ionen, die Oktaederlücken z. T. von den 
(Mg, Fe)-Ionen besetzt. Aus rein geometrischen Gründen ist stets die Anzahl 
der oktaedrischen Koordinationszentren einer dichtesten Kugelpackung 
gleich der Anzahl der vorhandenen Kugeln (in diesem Falle also der Anzahl 
der Sauerstoffionen). Olivin besitzt sonach pro Formeleinheit des Olivins 
vier Oktaederlücken, von denen jeweils zwei von den Kationen besetzt und 
zwei unbesetzt sind. Es ist also ohne stärkere Deformation der Sauerstoff- 
packung möglich, in den oktaedrischen Lücken alle drei Li+-Ionen derart 
unterzubringen, daß in jeder Lücke höchstens ein Lit sitzt. 

Andererseits sind nicht einmal die zwei Li-Ionen bzw. diadoche Kat- 
ionen zum Aufbau einer Olivinstruktur notwendig. BJOERLING und WeESsT- 
GREN (4) haben gezeigt, daß dem Triphylin unter Oxydation von Fe?* und 
Mn?°+ das Lithium ohne Änderung des Kristallgitters entzogen werden kann, 
und zwar bis zur Erreichung des Endproduktes Heterosit (Fe?+, Mn3+)PO,, 
das somit nunmehr die halbe Zahl Kationen zweiter Art enthält. 


Zu dieser Gruppe Triphylin isomorpher, Olivin isotyper Phosphate ge- 
hören: 


a bo ec ab 3e 
Ferrisicklerit 5,94 4,79 10,09 1.240 21 : 2,106 
Li<,(Fe?*, Mn?*)[PO,] 
Sicklerit 6,00 4,77 10,22 1,258 : 1 22,143 
Li<,(Mn?*, Fe2*)[PO,] 
Heterosit 5,82 4,76 9,68 15,223.217227034 


(Fe3*, Mn®*)[PO,] 
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Zwischen den Li-Mineralen mit Olivinstruktur, die hier erwähnt wur- 
den, besteht praktisch vollständige Mischbarkeit. Sie sind also unterein- 
ander isomorph und zum Olivin isotyp. 

U. W. sind Isotypiebeziehungen zur Phenakitgruppe bis jetzt noch 
nicht gefunden und mitgeteilt worden. 

Auch zum Verständnis der Verwandtschaft zwischen Strengit, Skorodit 
und Norbergit (Strunz (Zs. Krist. 99. 1938. 513)) ist die Möglichkeit des 
Ersatzes von 1 Fe durch 2 Li im Olivingitter von Bedeutung. Folgende 
Gegenüberstellung soll die strukturelle Beziehung verdeutlichen: 


Strengit Skorodit Norbergit 
Fe[PO,] : H,O Fe[AsO,] : H,O Mg;[(OH, F) | SiO,] 
a0 — 10,06 10,26 Co — 10,20 
bu, = 9,85 9,98 2b, = 9,40 
co —= 8,65 8,88 = 872 
Dy,—Pceab D. be DD. _ P men 

ZzE = 8 Zz=8 z= 4 


Nach STRUNZ sind im Skorodit die H,O und OÖ in einer wenig defor- 
mierten, hexagonal dichtesten Kugelpackung (analog dem Norbergit) an- 
geordnet. Im Norbergit sind die Oktaederlücken zur Hälfte, im Skorodit und 
Strengit dagegen nur zu 4 von Kationen besetzt. 


Andalusit-Typus 


Formeltyp: YEIYIS[A | zO,], mit A = O, OH. Nach Strunz (88) 
' besteht Isotypie zwischen Andalusit und Libethenit/Adamin (vgl. folgende 
' Übersicht!). 


Andalusit Libethenit Olivenit Adamin 
AL,[O | SiO,] Cw[OH | PO,] Cw[OH| AsO,] Zn,[OH | AsO,] 
201,76 8,43 8,62 8,31 
bo = 1,90 8,08 8,20 8,51 
Co = 5,56 5,90 5,94 6,06 
a, :bo : co = 0,982: :1:0,704  1,04:1:0,730 1,05::1:0,724 0,974 : 1 : 0,712 
zZ = 4 4 4 4 


12 
Raumgruppe D;, — P nnm. 


Die Struktur des Olivenits wurde von HERITSCH (34, 35) bestimmt. 
Allerdings wird die Beziehung zur Andalusitstruktur erst erreicht, wenn 
man bei Libethenit und Olivenit gegenüber der ursprünglichen Aufstellung 
a, und b, vertauscht. 

Nach RıcHmonxD (79) zeigt der Tarbuttit Zn,[OH | PO,] kristallchemi- 
sche Beziehungen zur Libethenitgruppe, die allerdings komplizierterer 
Art sind. 
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Granat-Typus. 


Formeltyp: X,$®!Y, [9 [z0,];. Isotyp mit Granat sind nach BUBECK 
und MACHATSCHKI (14) Berzeliit, nach McConneELL (70) und nach THıLo 
(100) die synthetischen Phosphate Na,Al,[PO,], und Ca,NaMg,[PO,];. Die 
synthetischen Phosphate wurden durch Zusammenschmelzen bei etwa 500°C 
hergestellt. Die Verbindung Na,Al,[PO,], entspricht strukturell dem Kalk- 
Eisen-Tongranat (Grossular) und die Verbindung Ca,NaMg;,[PO,], dem 
Kalk-Eisengranat (Andradit). Über die Beziehungen der hierher gehören- 
den Minerale gibt folgende Übersicht Auskunft: 


Grossular Andradit Griphit 
Ca,(Al, Fe),[SiO;]; Ca;Al,[SiO;]; (Na, Al, Ca, Fe),Mn,[(OH), | (PO,)3»;] 
a0 = 11,83—11,91 12,02 12,26 
Mg-Berzeliit Mn-Berzeliit 
(Ca,Na)(Mg, Mn),[AsO;], 
a0 #12,35 12,46 
Zi 5 Raumgruppe: on — la3d 


Griphit kann weitgehend zum Plazolith Ca;Al, (OH), | (SiO,),] parallel 
gestellt werden. Plazolith zeigt nach Untersuchungen von PABsr (73) große 
Ähnlichkeit mit Grossular. Die Plazolithstruktur kann als defektes Granat- 
gitter aufgefaßt werden, wobei Si statistisch auf 3 der 24fachen Punkt- 
lage des idealen Granatgitters verteilt ist, während 0° und OH- eine 96fache 
Punktlage im Verhältnis 2 : 1 besetzen. Im Griphitgitter dürften entspre- 
chend 4 der z-Plätze unbesetzt bleiben, wobei gleichzeitig 4 der O0?” durch 
OH- ersetzt werden. 


Es ist sicher, daß zwischen den Phosphaten mit Granatstruktur und den 
entsprechenden Silikaten keine ausgesprochene Mischbarkeit besteht. Im- 
merhin weisen die Analysen der Granate in einigen seltenen Fällen auf sehr 
geringe P,O,-Gehalte hin. 


Zirkon-Typus 


Formeltyp: X[®![zO,]. Vgl. auch Krokoittypus (8.179). Folgende Über- 
sicht stellt die hier in Frage kommenden Verbindungen zusammen: 


3 bo Co/&o Z Rer. 

Zirkon Zr[SiO,] . - . - 6,58 5,93 091 4 DW —IA4/amd 
Xenotim Y[PO,] . . . 6,88 6,03 0,877 4 ” z 
Y[PO,] !synth. . . . 6,862 6,174 0,90 
YASO Er 6,390 6,269 0,91 
IVO 7,126° 6,179 ° 0,867 4, 70, Aland 

1 STRADA & SCHWENDIMANN (84) oder 

® Broc# (13) Dr — l4/amd 


Nach Strunz (95) ist Xenotim nur in geringem Umfange isomorph mit 


Zirkon. Xenotim enthält in der Natur gelegentlich 5 und mehr %, SiO,, 
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wobei Sit* diadoch P®+ vertritt und der Ladungsausgleich durch Substitu- 
tion von Th*+ und Zr*+ (an Stelle von Y3*) hergestellt wird. 

Die geringe Mischbarkeit zwischen Zirkon und Xenotim ist bemerkens- 
wert, da beide Minerale häufig auf der gleichen Lagerstätte gemeinsam und 
sogar in regelmäßigen Verwachsungen vorkommen. Im übrigen scheint das 
Zirkongitter eine geringe Stabilität zu besitzen, worauf die Isotropisierungs- 
erscheinungen hinweisen. MACHATSCHKI (58) hat versucht, diese Frage auf 
die unregelmäßige Koordination um das Zr-Ion und die relativ großen Ab- 
stände zurückzuführen. 
| Homöotypie besteht offenbar zwischen Zirkon und KH,PO,, dessen 
© Struktur von West (107) eingehend bestimmt wurde (H 2,-Typ). Die PO,- 
Gruppen entsprechen den SiO,-Tetraedern, während die K-Ionen 8 Sauer- 
stoffionen in nahezu gleichen Entfernungen (zu zwei tetragonalen Disphenoi- 
den angeordnet) als Nachbarn besitzen. Die Gitterdimensionen können aus 
folgender Zusammenstellung entnommen werden: 


Zirkon Zr[SiO,] KH;[PO,] 
a0 = 6,58 7,43 
6 = 5,9 6,97 
Cola = 0,901 0,938 
Di, — I4/amd DD. 1494 


Titanit-Typus 

Formeltyp: XIYIS[A | zO,]. 

Auf die Isotypie zwischen Titanit und Tilasit wurde zuerst von STRUNZ 
(89) hingewiesen. Die Struktur von Durangit bestimmte KokKoros (49), 
der gleichzeitig auf die Verwandtschaft dieses Minerals mit Titanit auf- 
merksam machte. Untersuchungen an synthetischem Durangit wurden von 
MACHATSCHKI (59) durchgeführt. Eine Übersicht über die Gitterdaten der 
drei genannten Kristallarten vermittelt folgende Zusammenstellung: 


Tilasit Titanit Durangit 

Formel: CaMg[F | AsO,] CaTi[O | SiO,] NaAl[F | AsO,] 

&o 7,56 7,43 7,30 

bo 8,95 8,70 8,46 

Co 6,66 6,55 6,53 
20:D,.265 0,839 : 1 : 0,750 0,854 :1 : 0,753 0,856 : 1 : 0,772 

ß 120° 59 119 43 119 22 

va 386 367 352 

Z 4 4 4 


Raumgruppe: en — C2/Ja. 


Es besteht nach Strunz Isotypie zwischen Titanit und Tilasit, Titanit 
und Durangit und vermutlich auch zwischen Titanit und Cryphiolith (89a). 
Isomorph sind wahrscheinlich Tilasit, Cryphiolith und Durangit. RıcH- 
MOND (79) faßt im Gegensatz zu Strunz den Cryphiolith als Ca-reichen 
Wagnerit auf und betont, daß die Resultate von STRUNZ nicht mit seinen 
Beobachtungen übereinstimmen. 
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Eine weitere kristallchemische Verwandtschaft besteht zwischen den 
Mineralen der Tilasit- und denen der Descloizitgruppe. a* und b* sind für 
die Glieder der Tilasit- und Descloizitgruppe etwa gleich, während c* der 
Tilasitgruppe [301] in der Descloizitgruppe entspricht. Nach STRUNZ (9%) 
entsteht der Descloizittypus aus dem Tilasittypus vielleicht durch eine Art 
innere Verzwillingung nach (001). Im folgenden werden die Daten von 
Adelit und Tilasit, die im Verhältnis der Dimorphie zueinander stehen, ge- 


genübergestellt: 
Tilasit Adelit 
CaMg[OH | AsO,] CaMg[OH | AsO,] 
= 1,56 7,43 
DoE= 8:95 8,85 
= 6,66 
c.sin (180 — ß) = 5,71 co = 5,88 
a,:bo::Co = 0,839 :1 : 0,750 0,840 : 1 : 0,665 

0212027597 90° 

On c> ? (rhombisch) 


Thortveitit-Typus 


Formeltyp: Y;l® [z,0,]. Nach MACHATSCHKI (57) besteht Isotypie zwi- 
schen Thortveitit und Magnesiumpyrophosphat Mg,P,0, (vgl. folgende Zu- 
sammenstellung). Beim Thortveitit-Typus sind je zwei SiO,-Tetraeder zu 
einer Doppelgruppe verknüpft, beim Phosphat entsprechend der Struktur- 
formel: 


SS a 
0—P-0-P—0 


(6) “ N Ö 
Thortveitit Sc,[Si,0,] Ms;[P,0-] 
20656 13,28 — 2 x 6,64 
bes 28:58 8,36 
Co = 474 9,06 = 2 x 4,53 
ß= 103° 08 102112 
03 —C2 ? 
2h m 


MACHATSCHKI weist auf die Verdoppelung der a- und c-Achse beim 
Pyrophosphat gegenüber dem Thortveitit hin. Diese ist durch schwache 
Zwischenschichtlinien angedeutet. Bemerkenswert ist weiter, daß die Härte 
beim Übergang vom Silikat (H = 6—7) zum Phosphat (H — 4) im gleichen 
Maße sinkt wie beim Übergang vom Grossular (H — 7) zu dem mit ihm 
isotypen Arsenat Berzeliit (H — 4—5). Demnach ist für die Härte nicht die 
Bindungsintensität innerhalb eines Komplexes, sondern die Bindungsstärke 
zwischen Komplex und umgebendem Kation entscheidend. 

Weitere Phosphate mit [P,0,]-Gruppen sind von Levı und PEYRONEL 
(53) und PFYRONEL (78) untersucht worden. Es handelt sich um die isomor- 


phen Verbindungen Zr[P;0,], Si[P;0,], Ti[P;0;], Hf[P,O,], Sn[P;0,] und 
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U[P,0,]. Als Raumgruppe wurde Tn® — P a3 festgestellt. Isotype Silikate 
sind nicht bekannt (das Kation müßte 6wertig sein!). Der P—O-Abstand 
beträgt in diesen Strukturen 1,55 Ä, wie in den Inselphosphaten mit [PO,]- 
Tetraedern auch (vgl. z. B. Werrs (106) S. 418). Überraschend ist, daß 
Verbindungen mit so verschiedenartig sich vertretenden Kationen (Zr, Si, Ti!) 
eine ausgesprochen enge kristallchemische Verwandtschaft zeigen sollen. 


Von den Ringstrukturen, wie wir sie bei den Silikaten kennen, haben 
wir eine Analogie bei den Phosphaten für Dreier- und Viererringe. Nach Tuıto 
(120-124) ist Trimetaphosphat Na,[P,O,] isotyp mit Wollastonit Ca,[Si,0,]. 
Das Aluminiummetaphosphat — gewöhnlich als Al(PO,), geschrieben — be- 
sitzt nach PAULING und SHERMAN (76) Gruppen zu je [PO,]-Tetraedern als 
Baueinheiten, etwa nach folgendem Schema: 


(6) 16) 


0— P—-0—-P—O 


16) Ö 


Der P—O-Abstand innerhalb eines solchen Ringes beträgt 1,51 Ä. Die 
nicht zu zwei benachbarten [PO,]-Tetraedern gehörenden Sauerstoffionen 
umgeben die Al-Ionen oktaedrisch im Abstand 1,83 Ä. Der Gruppe kommt 
die Formel [P,O,5]*" zu, so daß dem Aluminiummetaphosphat die richtigere 
Formel Al,[P,O,s] zuzuschreiben ist. Nach Tuıro (120) sind auch beim 
Na,[P,O,,] Viererringe vorhanden. 

Zwar sind auch bei den Silikaten Viererringe bekannt, so z. B. beim 
Neptunit Na,FeTi[Si,0,,] oder beim Kainosit Ca,(Ce, Y)[CO, | Si,032]- 
1—2 H,0. Es scheint jedoch bis jetzt kein zu Al,[P,O,5] isotypes Silikat fest- 
gestellt worden zu sein. Eine gewisse Verwandtschaft könnte zum Zunyit 
Al,,[AlO, | (OH, F),sCl | Si,O,,] bestehen, wie folgende Gegenüberstellung 
andeutet: 


| | 
Ve) 

| 

Ö Ö 


A1,P,O Zunyit 

80 = 13,63 a0 — 13,82 
Te — 143 Ds .#43m 
Z= 16 Z=4 


Nach Pavrine (75) liegen im Zunyitgitter allerdings keine [Si,O75]- 
Viererringe, sondern Gruppen von je 5[SiO,]-Tetraedern vor. 


Unendliche Ketten der Zusammensetzung [PO,]) sind nach TrıLoz.B. 
beim MAnpreır'schen Salz (NaPO,), in Analogie zu den Pyroxenen vor- 
handen. Es ist weiter zu erwarten, daß solche Kettenstrukturen auch in Salzen 
wie Ca[PO,], vorliegen (WELLSs (106)). U. a. vermutet auch HENDRICKS (33) 
bei den strukturell noch nicht bestimmten Polyphosphaten Ca,P,O,, und 
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CaP,O,, Strukturen, die in ihrem Aufbau den polymeren [Si,0,ı Ketten 
der Amphibole entsprechen. 


Nach den Arbeiten von Strunz (87, 89) war bekanntgeworden, daß 
Epidot und Ardennit isotyp sind. Nun wurde bisher angenommen, daß Epi- 
dot ein Inselsilikat darstellt, obwohl gerade beim Epidot prismatische bis 
säulige Habitustypen vorherrschen. Tatsächlich kommt neuerdings Iro (39) 
auf Grund röntgenographischer Untersuchungen zu dem durchaus plau- 
siblen Ergebnis, daß beim Epidot ein Bandsilikat von ganz neuem Typus 
vorliegt. Die einfachen SiO,-Ketten sind beim Epidot etwas anders angeordnet 
als bei den Pyroxenen. Ferner sind beim Epidot je vier [SiO,]-Ketten mit- 
einander verknüpft, so daß [Si,O,]-Bänder nach folgendem Schema ge- 
bildet werden: 


Ö 10) (0) (6) 
Be (0) N 0) N DE 
Er 
0) \ Ö \ Ö ı Ö i (0) 
a 
0-520-4-0-8-0-8-0 
ee, 

(0) Y 10) \ (0) i (6) N Ö 
ON RE 


Dabei ist je eines von 4 Siliciumatomen der Bänder statistisch durch Al 
ersetzt. Die Bänder wiederum sind durch Ca-, Al- bzw. Fe- und Al-Ionen 
miteinander verknüpft: Ca in [8]-, Al(Fe) in [6]- und Al in [4]-Koordination. 
Letzteres wird weder durch Fe noch durch Si vertreten. Damit kann man 
dem Epidot folgende Formel zuschreiben: 


Ca,(Al, Fe)Al[0,0H | AlSi,O,] 
oder allgemein: 
X,[S)Y [9] A1#I[0,0H | 2,03]. 
Entsprechend ist die Formel für den mit Epidot homöotypen Ardennit 
Mn,Mn,Al,[(0,0H), | Al,Si,(P, As, V)O,s] - 2H;0. 


Nach Strunz (89) besteht keine Diadochie zwischen Si und (P, As,V), 
sondern es liegt ein ‚„Doppelsalz‘“ mit einem bestimmten stöchiometrischen 
Verhältnis Si: (P, As,V) vor. Nach Strunz ist ferner die Ardennitstruktur 
als durch innere polysynthetische Verzwillingung nach (001) aus Epidot ent- 
standen zu denken. Einen ähnlichen Aufbau zeigt auch Zoisit, wobei aller- 


dings die Verzwillingung nach (100) erfolgt. Der Strukturvergleich ergibt 
sich aus folgender Zusammenstellung: 


Epidot Ardennit 
80 = 8,96 8,72 
= 38 5,83 
eo 210,20 

2c, c08 (A—90%) — 18,43 18556 
BE 11507247 Br= 908 
VE, Ze): 
0 — P2,/m D!} — Pnmm 


Nach MACHATSCHKI (54) liegt dagegen beim Nagatelit eine Misch- 
kristallbildung vor, wobei sich Si und P diadoch vertreten. Nach den Unter- 
suchungen von Iro müßte dem Nagatelit die Formel 

(Ca, Ce)(Al, Fe),[O,0H | Al(Si,P);O,] 
zugeschrieben werden. Nagatelit wäre isomorph mit Orthit. 

Offen ist noch die Frage des Erikits, bei dem ebenfalls Si + P im Gitter 
eingebaut sind. STRUNZ (95, S. 188) stellt das Mineral zum Ardennit, wäh- 
rend MACHATSCHKI (61) seine kristallchemische Verwandtschaft zum Apatit 
hervorhebt. 


Zweidimensionale Vernetzungen sind bisher bei den Phosphaten u.W. 
noch nicht gefunden worden. Anklänge an Netzsilikate finden wir wohl 
beim Autunit Ca[UO, | (PO,),] : nH,O. Die Strukturschichten werden von 
[PO,]-Tetraedern und von [UO,]-Oktaedern gebildet. Die großen Kationen 
und Wassermoleküle besetzen die weiten Hohlräume, die zwischen den 
Netzen des Autunitgitters vorhanden sind. Untersuchungen über die kri- 
stallchemischen Zusammenhänge solcher Netzstrukturen, wie es das Autu- 
nitgitter darstellt und zu denen bei den Silikaten offenbar noch kein Kor- 
relat besteht, dürften recht interessante Resultate erzielen. 


Trimerit-Typus. 

Trimerit ist nach Untersuchungen von GOSSNER und BESSLEIN (25) 
isotyp mit Beryllonit. Als Formeltypus können wir X [z,0,] angeben. Die 
Gitterkonstanten der beiden Mineralarten stimmen recht gut miteinander 
überein (vgl. folgende Zusammenstellung). 
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Trimerit Beryllonit 
CaMn;[BeSiO,], Na[BePO,] 
a Ol 225%28,06 8,13 
I, = 100) 7,76 
CE 21,805 82,x.13,93 14,17 
Be2 bare, 2 1,060571°31,833 Bo boten . 1.04827212: 1,826 
0902092 90° 
Z7= 16 2. =.12 
Cm —P2y/e (2) Cy—P2le 


Fortschritte der Mineralogie 1949. Bd. 28. 13 
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Beide Minerale besitzen ausgesprochen pseudohexagonale Symmetrie 
mit (010) als pseudohexagonale Basis. Die Spaltbarkeit verläuft in beiden 
Fällen parallel (010). Die genaueren Zusammenhänge zwischen den beiden 
Strukturen sind noch untersuchungsbedürftig. 


Datolith-Typus. 


Nach Strunz (Zs. Krist. 93. 1936. 146) besteht Isotypie zwischen 
Datolith und Herderit. Die Strukturformel für diesen Gittertypus dürfte 
lauten: X[A | z,0,]. 


Datolith Herderit 
Formel Ca[OH | BSiO,] Ca[(OH, F) | BePO,] 
2 9,64 9,80 
bo 7,62 7,68 
© 4,82 4,80 
DEAD ORICH 1,266 : 1 : 0,633 1,276 :1 : 0,625 
ß 90° 09° 90° 06° 
\ 352 314 
Z 4 4 
Rer. Ch —P 2ıle 0) —P2ıle 


Zu Datolith isotyp bis isomorph sind die beiden Minerale Homilit 
Ca,Fe[O | BSiO,], und Gadolinit Y,Fe[O | BeSiO,]>- 

Beim Datolithgitter liegen dreidimensionale Vernetzungen von[BO,P- 
und [SiO,]--Tetraedern bzw. [BeO,]°-- und [SiO,]-- bzw. [PO,P--Tetra- 
edern vor. 

Nach RıcHMmoND (79) besteht kristallchemische Verwandtschaft zwi- 
schen Amblygonit (Li,Na)[(F,OH) | AlPO,] und Herderit. Die Beziehung 
ist allerdings komplexer Art. 

Auch die Beziehung zwischen Trimerit- und Datolith-Typus bedarf 
noch weiterer Untersuchungen. 


SiO,-Modifikationen. 


Zwischen verschiedenen Phosphaten und den SiO,-Modifikationensind 
weitgehende Isotypiebeziehungen gefunden worden. Hierher gehören: 
I. BPO, und BAsO,-Isotypie mit Hoch-Cristobalit nach SCHULZE (82), 
II. FePO,-Isotypie mit Tiefquarz nach CAGLiort (16). 

III. AlPO,- und AlAsO,-Isotypie mit Tief-Cristobalit nach STRADA 
(83), ebenso nach CAGLIorTi (16), 
AIPO,-Isotypie mit Tief-Quarz nach HUTTENLOCHER (38), 
Berlinit-Isotypie mit Tief-Quarz nach STRUNZ (91), 
AlAsO,-Isotypie mit Tief-Quarz nach MACHATSCHKI (56), 


AIPO,-Isotypie mit allen Modifikationen von SiO, nach WINK- 


HAUS (109). 


I. Die Angaben von SCHULZE (82) können nur als Beschreibung einer 
dem Tief-Cristobalit ähnlichen Struktur angesehen werden. ae selbst 
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weist darauf hin, daß das BO,-Tetraeder, das im Wechsel mit dem PO,-Te- 
traeder die Lagen der SiO,-Tetraeder in der Cristobalitstruktur vertritt, 
sehr stark deformiert ist. Auch die Gitterkonstanten von BPO, stimmen 
nicht besonders gut mit denen des Tief-Cristobalit überein: 


ao Co Co/&o 
BIRNEN 4,332 + 0,006 6,640 + 0,008 1,533 + 0,014 
Tief-Cristobalit SiO, 4,96 6,92 1,395 


II. CacLıoti (16) weist Isotypie zwischen FePO, und AlAsO, nach. Das 
Arsenat zeigt nach MACHATSCHKI (56) eine Quarzstruktur mit verdoppelter 
c-Achse. 


III. STRADA (83) fand für die Struktur des AIPO, die tetragonale Raum- 
gruppe S,— I und bestimmte die Gitterkonstante a, — 4,85+ 0,005 Ä 
und c,/a. zu 1,362 mit zwei Formeleinheiten in der Elementarzelle. Gleich- 
zeitig untersuchte STRADA auch AlAsO, und stellte völlige Isomorphie zwi- 
schen den beiden Verbindungen fest. 

CAGLioTI deutete die Struktur von AlPO, als pseudotetragonal, wahr- 
scheinlich rhombisch. Er fand in anderer Aufstellung als STRADA a, —= 7,135, 
co = 6,850 Ä mit vier Formeleinheiten in der Elementarzelle. 

Im Gegensatz hierzu gibt HUTTENLOCHER (38) für AlPO, Quarzstruktur 
mit verdoppelter c-Achse an. Er bestimmte die Gitterkonstanten zu: 


ao Co Co/ao 
AIPO, ren ar a en RA: 4,93 SIE 0,01 Dr 5,47 SIT 0,01 IR TI 
OQusraSI Ober: 4,903 5,393 1,092 


Raumgruppe: D% bis D. 


Nach diesen Ergebnissen darf die Übereinstimmung im Kristallbau von 
Aluminiumphosphat und Tiefquarz als erwiesen angesehen werden. Dabei 
ist allerdings zu beachten, daß der Ersatz der Si-Atome des Quarzes je zur 
Hälfte durch Al®+ bzw. P5+ eine Verdoppelung der Elementarzelle in Rich- 
tung der c-Achse bedingt. 

Die Untersuchungen von STRUNZ (91) an natürlichem AIPO, (Berlinit) 
stimmen nahezu völlig mit den HUTTENLOCHER’schen Resultaten überein. 
Strunz fand a, = 4,92, & = 10,91 = 2 x 5,455 Ä und «Ja. = 2,27 
—2x 1,135. Z = 3 und Raumgruppe: D; — 03,2. STRUNZ konnte die 
Frage, ob Mischkristalle zwischen Berlinit und Tiefquarz möglich sind, nicht 
entscheiden. j 

In einem Vortragsbericht von BUERGER (15) wird darauf hingewiesen, 
daß für AIPO, drei SiO,-Modifikationen (Quarz, Tridymit, Cristobalit) be- 
kannt sind. Leider fehlen dort weitere Angaben. 

In einer ausführlichen Untersuchung von WINKHAUS (108, 109) (vgl. 
auch TRÖMEL und WınkHaAus (104)) wird festgestellt, daß AlPO, in allen 
sechs Modifikationen des SiO, (Tiefquarz, Hochquarz, Tieftridymit, Hoch- 
tridymit, Tieferistobalit, Hocheristobalit) auftreten kann. Die DEBYE- 
SCHERRER-Aufnahmen der verschiedenen Strukturen von AlPO, und von 
SiO, stimmen bei Zimmertemperatur gut miteinander überein. 
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Je nach den bei der Herstellung gewählten Temperaturen erhält man 
das AIPO, bis zu 800° C mit Quarz-, bis 1100° C mit Tridymit- und darüber 
mit Cristobalitstruktur. Erhitzt man das bei 500—600°C gebildete AIPO, 
mit Quarzstruktur über 800° C, so bildet sich die Tridymitstruktur, die sich 
bei 1100° C ebenso ohne Verzögerung in die Cristobalitstruktur umwandelt. 
Die Umwandlungsgeschwindigkeiten des AlPO, sind merkbar größer als die 
des SiO,. Bei einigen Herstellungsverfahren wurden AlPO,-Produkte er- 
halten, die sich vom amorphen Zustand nach thermischer Behandlung so- 
fort in die Tridymit- oder Cristobalitstruktur umwandelten. 

Daß auch die verschiedenen AlPO,-Modifikationen die gleichen glatt 
ablaufenden «--Umwandlungen zeigen, die beim SiO, bekannt sind, konnte 
ebenfalls bei Temperaturen bis 1200°C festgestellt werden. Die Umwand- 
lungstemperaturen sind beim AIPO, im allgemeinen etwas niedriger als beim 
SiO,. Beim kubischen «-AlPO,-Cristobalit wurde die Gitterkonstante bei 
250°C zus. = 7,11 +0,01 Ä ermittelt. Sie ist innerhalb der Fehlergrenze 


der des «-SiO,-Cristobalits gleich (7,12 Ä nach NIEUVENKAMP, Zs. Krist. 96. 


1937. 454). 

Auch von AlPO,-Tridymit wurden Hochtemperatur-Aufnahmen meh- 
rerer Modifikationen parallel zum SiO,-Tridymit gemacht, die aber wenig 
ausgeprägte Unterschiede zeigen. In diesem Zusammenhange sind noch 
nähere Untersuchungen erforderlich. 

Erst nach Fertigstellung unserer Untersuchungen ist uns eine ausführ- 
liche Arbeit von BEck (114) bekanntgeworden, die über röntgenographische 
und thermische Versuche mit den AlPO,-Modifikationen berichtet. Grund- 
sätzlich finden die Ergebnisse von WInKHAUS (109) ihre Bestätigung. Auch 
Beck fand, daß die Beziehungen zu SiO, ziemlich eng sind und daß die Um- 
wandlungen viel rascher erfolgen als bei den entsprechenden SiO,-Modi- 
fikationen. Außerdem sind die Umwandlungstemperaturen (mit Ausnahme 
bei Berlinit) niedriger als bei SiO,. BEcK stellt folgendes Umwandlungs- 
schema auf: 

AIPO, 


815° +49 i 1025°+50° >1600 
Berlinit <——> Tridymit-Form <- — Cristobalit-F. <—> AlPO,-Schmelze 
5860+20 93043 130° (?) 2100-+5° 
Px<———>a B<— <> 0 Bx——& 


Und in Parallele dazu: 
S18:0, 
70° 1470° 1710° 
Quarz <- — Tridymit <— —- Cristobalit <— 
5780 117 163° 2202700 
Pa Pe > Pa 


—- SiO,-Schmelze 


Weiter hat GRUNER (117) eine Arbeit veröffentlicht, in der über orientiertes 
Weiterwachsen von AIPO, auf Quarz berichtet wird. Bei 150—250° traten 


schneeflockenförmige Kristalle von AIPO, in Tridymitform auf. Bei 300 bis 
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480° wurden AIPO,-Kristalle mit Cristobalittypus beobachtet, während 
sich nach GRUNER Berlinit im Temperaturbereich von 850—1000° in die 
Cristobalitform umwandelt. 


Schließlich konnte die Frage der Mischkristallbildung zwischen AlPO, 
und SiO, auch durch die Untersuchungen von WINKHAUS noch nicht ge- 
klärt werden. 


Eine Notiz von BRILL und DE BRETTEVILLE (12) beschäftigt sich mit 
dem Bindungscharakter von AlPO,. Die Untersuchungen führten zu dem 
Ergebnis, daß im AlPO,-Gitter größere heteropolare Anteile vorliegen. 


Abschließend sei ausdrücklich betont, daß dreidimensionale Vernetzun- 
gen von PO,-Tetraedern allein bisher offenbar noch nicht festgestellt 
worden sind. 
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